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ABSTRACT 


This report contains a brief description of the Stromberg Datagraphixs 4060 
(SD-4060) Orbit Comparison Plot (OCPLT4) Program, along with user informa- 
tion and a source program listing. This program was developed by Computer 
Sciences Corportation under Task Assignment 096 to supersede the SC-4020 
OCPLT4 Program, which was developed in early 1970. The object program is 
currently on tape number 564M, and filed under Program Number 498 at GSFC 

Program library. 

In addition to correcting several errors that existed in the original program, 
this program incorporates the following new features. 

• For any satellite whose observations are processed by the Definitive 
Orbit Determination System (DODS), the orbital uncertainty esti- 
mates (OUE) can be obtained via appropriate card input with no major 
modification to the program. 

• All satellite-related information (e. g. , plotter scales, cutoff limits, 
plotting frequencies) is user controlled via card input. 

• Not all components of OUE must be obtained. The user has the 
option of obtaining only the radial component if there is no need for 
the other two components. 

• The altitude and time graph formats are controlled by the user and 
are not stored for specific satellites. 
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SECTION 1 - INTRODUCTION 


The purpose of the SD-4060 OCPLT4 Program is to generate an instruction tape 
for the Stromberg Datagraphixs 4060 (SD-4060) plotter (see Reference 1). The 
resulting graphs display component differences between two satellite position 
vectors within an overlapping time period. These differences are called the 
orbital uncertainty estimates (OUE). The following set of three orthogonal vec- 
tor components is plotted: 

• The radial component 

• The component normal to the radial component in the orbital plane 

• The component normal to the radial component and normal to the 
orbital plane 

Each component is plotted on an individual graph. 

The components can be plotted on either a linear or a logarithmic ordinate scale, 
and against an altitude or a time abscissa scale. The choice of abscissa scale 
is a function of the satellite's altitude. Low-altitude satellites with several rev- 
olutions per day are usually plotted against a time scale, and the results are 
called time graphs, whereas high-altitude satellites which complete only one 
revolution in several days are usually plotted on an altitude scale, and the re- 
sults are called altitude graphs. When requested, these graphs also display the 
time distribution plots of observations used in obtaining the converged elements 
that provide the overlapping ephemerides. These graphs provide the experi- 
menter with the OUE that can be used for analyzing definitive orbit results 
(see Reference 2). 

OCPLT4 provides OUE graphs for any satellite whose observations are proc- 
essed by the Definitive Orbit Determination System (DODS) on the IBM Sys- 
tem/360. 
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Inputs to the SD-4060 OCPLT4 Program include the vector compare (VC) tapes, 
which are generated by DODS Ephemeris Comparison Subsystem; and a working- 
observations -file tape, which is generated by the DODS Differential Correction 
(DC) Subsystem. 

Output from the SD-4060 OCPLT4 Program consists of a printout detailing what 
was accomplished by the run, and an instructions tape for the SD-4060 plotter 
to plot the OUEs. Usually, the SD-4060 plotter will provide 16-mm microfilm 
frames, one frame for each OUE graph, although 35-mm can be requested. 

Hard copies can be obtained from either film format upon request. 

The SD-4060 OCPLT4 Program has been compiled under FORTRAN IV, level H, 
optimization level 2, on the Goddard Space Flight Center (GSFC) IBM System/360 
Model 95, operating under OS using Release 19. 6. No changes are necessary 
to run this program on the M&DO IBM System/360 Model 75. 

The remaining sections of this user's guide present detailed information on pro- 
gram input (with sample deck setup), program output (including error messages), 
sample plotter output graphs, and operating information (with timing estimates). 
Also presented are the programming approach utilized, brief descriptions of 
subroutines, and a source program compilation listing. 
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SECTION 2 - PROGRAM INPUT 


2.1 USER OPTIONS 

All satellite-related variables are user controlled in this version of OCPLT4. 
Variables include satellite name, ID number, and date of run, all of which 
appear on the plots. Grid labeling and grid spacing are also user controlled to 
provide the flexibility required to process a wide variety of satellites. Other 
user inputs are the upper and lower cutoff limits for graphs. These inputs 
allow the user to control the overall appearance of the plots. 

The user controls the following in a single job submission: the type of abscissa 
(altitude or time 1 ); the type of ordinate scale (linear or logarithmic); whether 
or not observation data distribution will be plotted; and whether the radial com- 
ponent only, or all three OUE components, will be plotted. 

2. 2 TAPE INPUT 

OCPLT4 requires at least two input tapes. The first, the VC tape, is generated 
by DODS using function 1 of the COMPARE verb (see Reference 3). This is a 
nine-track EBCDIC tape which is loaded on any 2400 series tape drive. It con- 
tains the Orbit Comparison Report (see Reference 3). This report is obtained 
by comparing two overlapping ephemerides (satellite-position time histories). 
Both ephemerides must be generated at equally spaced and corresponding time 
points in the overlap region. The differences between the two satellite position 
vectors at each point in time are expressed as differences between three orthog- 
onal components of the vectors. The Orbit Comparison Report consists of a 
tabulation of the two ephemerides, the three component differences (which are 
the OUEs), and the total vector difference as a function of time. Several Orbit 
Comparison Reports (also called VC Reports) could be written onto a single file 


1 Either plot, or both, may be generated from a single job submission. 
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on a tape, and several files could be written onto a single tape, but they must 
appear in ascending time order. OCPLT4 will process up to 24 VC files in a 
single job submission. 

The second tape is the working-observations -file tape. This is a nine-track 
binary tape which is likewise loaded on any 2400 series tape drive. It contains 
the working observations data (see Reference 4), as generated by DODS using 
the SETDC verb (see Reference 3). It should contain the observations from a 
time period which extends by at least three hours on both sides of the period 
covered by all of the VC files to be plotted. These observations could be in con- 
centric or nonconcentric time order. This is determined by the relationship 
between the epoch of elements and the start time used in creating the working- 
observations-file tape. When the epoch date precedes, or is equal to, the start 
time of data, observations will be in ascending time order (nonconcentric). If 
epoch is between start and end time, observations will be in concentric order. 

In case no data distribution plots are requested, a tape must Still be mounted; 
it may be a dummy tape. When using the SETDC verb for this purpose, the 
standard DODS Job Control Language (JCL) should be overridden so that the 
working-observations-file data are output on tape instead of disk. 

2. 3 CARD INPUT 

At least 15 data cards are required for each OCPLT4 run. These cards must 
appear in the data deck in the order indicated by card number (Card 1, Card 2, 
etc.). The format for each card is defined on the following pages. 
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CARD 1 


Format 

Column 

Internal 

Variable 

Name 

Description 

A1 

1 

CON 

Indicates whether data on working ob- 
servations file tape is concentric 
edited: 

= C, concentric edited 
f C, not concentric edited 

IX 

2 


Blank 

A8 

3-10 

SNAME 

Satellite name (e. g. , SSS-1) (left 
justified) 

IX 

11 


Blank 

15 

12-16 

ISAT 

Satellite identification no. (e. g. , 
71961) 

IX 

17 


Blank 

16 

18-23 

IRUN 

Computer run date in YYMMDD for- 
mat (e. g. , 720912) 

IX 

24 


Blank 

11 

25 

LOG 

Indicates type of scale on Y-axis of 
graph: 

= 0, linear scale 
= 1, log scale 

IX 

26 


Blank 

11 

27 

MANY 

Controls labeling interval for the hours 


scale (X-axis) on the data distribution 
plot when altitude graphs are desired: 

= 0, label every hour 
= 1, label every 4 hours; this prevents 
overcrowding of the hours labels 
and as a rule should be used when 
there are more than 2 days be- 
tween apogee and perigee 
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CARD 2 


Format 

Column 

Internal 

Variable 

Name 

Description 

F10.0 

1-10 

XL 

Lower limit of X-coordinate on altitude 
graph, thousands of km 

F10. 0 

11-20 

XR 

Upper limit of X-coordinate on altitude 
graph, thousands of km 

F10.0 

21-30 

YB1 1 

Lower limit of Y-coordinate on altitude 
or time graph, radial component (km) 

F10.0 

31-40 

YB2 1 

Same as above except for in-plane 
component 

F10.0 

41-50 

YB3 1 

Same as above except for normal-to- 
plane component 

F10.0 

51-60 

YT1 1 

Upper limit of Y-coordinate on altitude 
or time graph, radial component (km) 

F10.0 

61-70 

YT2 1 

Same as above except for in-plane 
component 

F10.0 

71-80 

YT3 1 

Same as above except for normal-to- 


plane component 


!when the log mode is being used, 


these limits must be integer powers of 10. 
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CARD 3 


Format 

F10.0 

F10. 0 

F3.1 


NOTE: 


Internal 

Variable 

Column Name Description 


1-10 XGRID 


11-20 XLABEL 


21-23 FMTX 


Length of interval for drawing grid 
along the X-axis, altitude option only 
(thousands of km) 

Length of interval for labeling grid 
along the X-axis, altitude option only 
(thousands of km) 

Format for labeling grid along the 
X-axis, altitude option only. FMTX 
is of the form W. D, where W is the 
maximum number of characters in a 
label, including decimal point but not 
the sign; and D is the number of places 
to be displayed to the right of the deci- 
mal. If the X-axis were to be labeled 
from 0. to 140. , FMTX would be 4. 0. 


This card must be included, but should be left blank when using the time 
option only. 
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CARD 4 


Format 

Column 

Internal 

Variable 

Name 

Description 

F10.0 

1-10 

YGRID1 1 

Length of interval for drawing grid 
along the Y-axis of altitude or time 
graph, radial component (km) 

F10.0 

11-20 

YGRID2 1 

Same as above except for in-plane 
component 

F10.0 

21-30 

YGRID3 1 

Same as above except for normal-to- 
plane component 

F10.0 

31-40 

YLAB1 1 

Length of interval for labeling grid 
along the Y-axis of altitude or time 
graph, radial component (km) 

F10.0 

41-50 

YLAB2 1 

Same as above except for in-plane 
component 

F10.0 

51-60 

YLAB3 1 

Same as above except for normal-to- 
plane component 

F3.1 

61-63 

FMTY1 

Format for labeling grid along the 
Y-axis of altitude or time graph, 
radial component (see FMTX on 
card 3) 

IX 

64 


Blank 

F3. 1 

65-67 

FMTY2 

Same as above except for in-plane 
component 

IX 

68 


Blank 

F3.1 

69-71 

FMTY3 

Same as above except for normal-to- 
plane component 

IX 

72 


Blank 


■’'These fields should be left blank when using the log mode, since the log mode 
provides its own grid generation and labeling for the Y-axis. 
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CARD 5 


Format 

Column 

Internal 

Variable 

Name 

Description 

F10.0 

1-10 

ERRL01 

Lower cutoff limit. If the radial com- 
ponent is less than ERRL01, the com- 
ponent is set equal to ERRL01 and 
plotted. ERRL01 is in km. 

F10.0 

11-20 

ERRL02 

Same as above except for in-plane 
component 

F10.0 

21-30 

ERRL03 

Same as above except for normal-to- 
plane component 

F10.0 

31-40 

ERRHI1 

Upper cutoff limit. If the radial com- 
ponent is greater than ERRHI1, the 
component is divided by 10 before 
plotting, and an appropriate message 
is displayed on the plotted output. 
ERRHI1 is in km. 




If, after dividing by 10, the value of 
ERRHIl is still exceeded, data are 
plotted outside the graph (user should 
then increase the scale accordingly 
and resubmit this run). 

F10.0 

41-50 

ERRHI2 

Same as above except for in-plane 
component 

F10.0 

51-60 

ERRHI3 

Same as above except for normal -to- 


plane component 
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CARD 6 


Format 

Column 

Internal 

Variable 

Name 

Description 

16 

1-6 

IDAT 

YYMMDD of start time of period to be 
plotted 

IX 

7 


Blank 

16 

8-13 

IDAT1 

YYMMDD of end time 
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CARD 7 


Internal 

Variable 

Format Column Name Description — 

j 4 4_4 IH Hours and minutes of start time 

(HHMM), where 

HH = hour-of-day 

MM = minute -of-hour 

(Cannot precede start time on first 

VC report to be plotted) 
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CARD 8 


Internal 

Variable 

Format Column Name Description 

A4 1-4 TIMEY"'’ Indicates type of graph to be plotted. 

If = TIME, only time graphs are 
plotted; if left blank, both altitude 
and time graphs are plotted. 


^■For time graphs only— user must specify TIMEY = TIME and NSS6 = 0 or 
blank. For altitude graphs only — user must leave TIMEY blank and specify 
NSS6 = 1 (see Card 9). 
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CARD 9 


Format 

Column 

Internal 

Variable 

Name 

Description 

11 

1 

NSS1 

Dummy, leave blank 

11 

2 

NSS2 

Dummy, leave blank 

11 

3 

NSS3 

Dummy, leave blank 

11 

4 

NSS4 

Data distribution flag: 

= 1, eliminate data distribution part 
of graphs 

= 0, do not eliminate data distribution 
part of graphs 

11 

5 

NSS5 

Debug printout flag: 

= 1, suppress debug printout 
= 0, do not suppress debug printout 

11 

6 

NSS6 1 

Graph flag: 

= 1, suppress time graphs 
= 0, generate both altitude and time 
graphs 


1 For time graphs only — user must specify TIMEY = TIME and NSS6 = 0 or 
blank. For altitude graphs only — user must leave TIMEY blank and specify 
NSS6 = 1. 


2-11 



CARD 10 


Format 

Column 

Internal 

Variable 

Name 

Description 

11 

1 

NCOMP 

Indicates number of range difference 
vector components to be plotted: 

= 1, only the radial component is 
plotted 

= 3 , all three components are plotted 

IX 

2 


Blank 

F3.0 

3-5 

TFREQ 

Plotting interval for time graphs (min- 
utes) (Equals the frequency of selecting 
points from VC report, must be inte- 
gral multiples of T3DIFF x 60) (See 
Card 11) 

IX 

6 


Blank 

F8.0 

7-14 

APOGEE 1 

Satellite apogee (to nearest km) 

IX 

15 


Blank 

F8.0 

16-23 

PERIGE 1 

Satellite perigee (to nearest km) 


^Used to determine the plotting interval for altitude graphs. 



CARD 11 


Format 

Column 

Internal 

Variable 

Name 

Description 

F4.0 

1-4 

T3DIFF 

Time between comparison points, in 
seconds (available from VC output) 
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CARD 12 


Format 

Column 

Internal 

Variable 

Name 

A6 

1-6 

TAPE 1 

IX 

7 


12 

8-9 

IFILE 2 


Description 

VC tape number 
Blank 

Number of VC reports on this tape 


^One card per tape must be specified for each VC tape number for any com- 
bination of tapes and files on tape, up to 24 files. There may be more than 
one file per tape. 

2 There may be more than one VC report per file and more than one file per 
tape. IFILE = total number of VC reports on the specified tape. 
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CARD 13 


Internal 

Variable 


Format 

Column 

Name 

Description 

None 



Blank card; delimiter for comparison 
tapes 
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CARD 14 


Format 

Column 

Internal 

Variable 

Name 

Description 

A6 

1-6 

TAPE 

Working-observations -file tape number 

IX 

7 


Blank 

12 

8-9 

IFILE 

= 1 


(Only one file will be processed per 
single submission. The time span of 
data must extend on both sides of the 
total VC reports time span. ) 
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CARD 15 


Internal 

Variable 

Format Column Name Description 

None Blank card; indicates end of card 

input. 

NOTE: There can be more than 15 data cards because card 12 may be repeated 
up to 24 times. 

Appendix A provides a sample deck input. 
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2.4 FILES 


The SD-4060 0CPLT4 Program uses only one file, a temporary disk data set, 
FT22F001. This data set contains time-sorted information from the working- 
observations-file tape and is used in plotting the data distribution portion of 
altitude or time graphs. There is a record for each observation. These rec- 
ords have the following format. 


Internal 

Variable 

Format Name Description 


IX 

A4 


3X 

16 

IX 

14 

IX 

12 


NOSTOP 


ITIME8 

ITIME9 

ITYPE 


Blank 

Indicates end of data. In the last record, 
NOSTOP is blank. In all other records, 
NOSTOP is ABCD. 

Blank 

YYMMDD 

Blank 

Hour-of-day and minute -of- hour (HHMM) 
Blank 

Type of observation: 

= 1, R range data 

= 2, H ) mini track direction cosines 

= 3, m I data 

= 9, R range-rate data 

= 17, RAO-X) radio antenna observation 

= 18, RAO-Yl angles data 
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SECTION 3 - PROGRAM OUTPUT 


3.1 TAPE OUTPUT 

The program's output is a seven-track binary instruction tape (data set 
SC4060ZZ), which is used as input to the SD-4060 plotter. The format of this 
tape is described in Reference 5. 

3.2 SYSTEM PRINTER OUTPUT 

This section presents information on normal printer output and on error mes- 
sage output. 

3.2.1 Normal Printer Output 

As processing is initiated, the program prints out some of the input variables 
to enable the user to spot check possible input errors along with the start and 
end times of the first VC report to be plotted. As processing proceeds, the 
first task of the program is to rearrange the concentric sorted observations 
from the working-observation-file tape in ascending time order, when neces- 
sary. The time and type of the rearranged observations are printed out. Each 
rearranged time and type is preceded by the letters "ABCD. " 

When both altitude and time graphs are requested along with the data distribu- 
tion plots, as in the sample output (see Appendix B), the program will first plot 
altitude graphs for each of the three OUE components from the first VC report, 
with the associated data distribution information and then the time graphs with 
the associated data distribution graphs. To indicate that the program has fin- 
ished reading the VC report, a flag "AT 1004” is printed. The backspacing of 
this VC report, needed when both altitude and time graphs are requested, is 
shown by A3COMP=7 until the VC report is backed to the first data point. 

ITGPH indicates that the time plots will be plotted next. This sequence is re- 
peated until all VC reports have been processed. 
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The main portion of the printout is concerned with the data distribution portion 
of the graphs. Because the data distribution is identical for all three OUE com- 
ponents, they appear in triplicate. When altitude graphs are plotted, the portion 
of orbit being plotted is indicated by apogee-to-perigee (A-P) or perigee-to- 
apogee (P-A) pass. The time span between A-P or P-A is indicated by the 
PERIOD PLOTTED, and the YYMMDD HHMM of the start and end times of the 
period, and also by the integer hour difference between the start and end times. 
For time graphs, this period is fixed at 24 hours. 

The type and quantity of data available from the working-observation-file tape 
during the period being plotted is also indicated. The number denoted in the 
message ". . . PASSES PLOTTED xx" refers to the number of minutes contain- 
ing one or more observations from a single station. Thus, if within 1 minute, 
one or more observations were obtained from one station; the number of passes 
is increased by one and one asterisk is plotted in the data distribution plot at a 
location corresponding to the hour and minute of the observation. 

Refer to Appendix B for a more detailed description of normal printer output. 

3.2.2 Debug Output 

As a further aid to the user, debug printout will be displayed if column 5 on 
data card 9 is 0 or is left blank. 

This printout, which supplements the normal printout, contains several flags 
to help identify where in the program the computation takes place, the values 
of several key variables, the computed location in plotter units of the first point 
to be plotted, the hours for altitude plots, and the geocentric distances for the 
time plots and their coordinates on the respective graphs. 

3.2.3 Error Message Output 

If the start year-month-day of the current VC report is greater than the end 
year-month-day of the previous report, the message TIME SPAN INCORRECT 
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ON THIS VC REPORT will be printed, along with the start and end year-month- 
day in question. Finally, the message OCPLT4 WILL PROCEED TO NEXT VC 
REPORT TO SEARCH FOR CORRECT TIME SPAN will be printed. 

If the time period to be plotted extends beyond the end time of the time-sorted 
working observations file information on the temporary disk data set, the mes- 
sage REQUESTED TIME SPAN TO BE PLOTTED EXCEEDS OBSERVATION 
TIME will be printed, and program execution will terminate. 

The Integrated Graphics Software (IGS) System is a subroutine library used by 
OCPLT4 to generate an instruction tape for the SD-4060 plotter. In the event 
that OCPLT4 gives an illegal command to the IGS System (such as a command 
to plot a number off scale), an appropriate error message from the IGS System 
will be printed. A complete list of these error messages can be found in 
Table 3-4 of Reference 1, and is reproduced verbatim in Appendix C. 

3.3 GRAPHIC OUTPUT 

The final products of this program are graphs which display the OUE for an 
orbit, along with the data distribution information. This section describes the 
two types of graphs (altitude and time) generated by the OCPLT4 Program. 
Appendix B illustrates a complete set of altitude and time graphs. 

3.3.1 Altitude Graphs 

Altitude graphs are usually requested when a satellite's orbital period is greater 
than 24 hours. An altitude graph presents the OUE components as a function of 
geocentric distance, and also includes a separate data distribution plot. 

Six altitude graphs are normally generated for each orbital period. The ab- 
scissa of each graph represents the satellite's radial distance from the center 
of the earth in 1000-km units. The first three graphs (Figures B-2 through B-4) 
are plotted for one-half an orbit, from apogee to perigee; and the other three 
graphs (Figures B-5 through B-7) are plotted for the remaining half of the orbit, 
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from perigee to apogee. The ordinates of the three graphs for each half orbit 
are the three components of the range difference vector: the radial component, 
the component in the orbital plane normal to the radial component, and the com- 
ponent normal to the orbital plane. Two altitude graphs will be generated when 
user specifies radial component only. 

At the bottom of each graph is a separate plot, which is produced at the user's 
request. The observations that are available from the working-observations- 
file tape are represented on this plot versus universal time (UT) (see Fig- 
ures 3-1 and 3-2). The time span of this plot corresponds to half of the orbital 
period. Asterisks represent the data distribution for three sets of observation 
types: radio antenna observation (RAO) X and Y angles; Minitrack direction 
cosines (l and m); and range and range rate (R and R). The asterisks become 
darker as more observations are available at a given time, as from several 
stations (see Figure 3-3). When no observations are available on the tape, 
during the time interval of a plot, or if a blank observation tape is mounted, 
the message NO DATA FOR THIS PERIOD will appear in place of the asterisk? 
(see Figure 3-2). 

The following user input information appears in the title of each graph: 

Run date (e.g. , 720912) 

Satellite name (e.g., SSS-1) 

Satellite ID (e.g., 71961) 

The grid spacing and coordinate labeling are user controlled. 

The numbers along the altitude OUE curve indicate UT in hours of day along 
the trajectory. These hour numbers start with the hour of apogee and end with 
the hour of perigee for the apogee-perigee graphs, and are in reversed order 
for the perigee-apogee graphs. As would be expected, these times are gener- 
ally not equally spaced. 
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MISSION AND TRAJECTORY ANALYSIS DIVISION, OODDARD SPACE FLIGHT CENTER RON DATE 720912 

ORBITAL UNCERTAINTY ESTIMATE FOR SSS-1 < 7 '«i > 
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Figure 3-1. Altitude Graph of Perigee-Apogee Radial Component 
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Figure 3-2. Altitude Graph of Apogee-Perigee 
Component Normal to Orbital Plane 
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Figure 3-3. Time Graph Component in Orbital Plane Normal to Radial 
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A similar time span is printed in the data distribution plot; however, the time 
spacing in this plot is uniform. The units here are also UT, and they corre- 
spond to the hours of the day for the data distribution. In addition, the obser- 
vation dates appear on the data distribution plots. A date is printed for every 
computed day within the trajectory's time span (see Appendix B, Figure B-l). 

3.3.2 Time Graphs 

Time graphs are usually requested for satellites with short orbital periods (two 
or more revolutions per day). The time graph is basically similar to the alti- 
tude graph, with the following exceptions: 

• Three graphs are normally generated, one for each component of the 
OUE, for each 24-hour period, starting at midnight UT. 

• The abscissa of each graph is divided into hours of day UT. 

• The numbers which appear along the OUE curves indicate radial 
distance from the center of the earth along the trajectory, in 1000-km 
units. 

The remainder of the graph is similar to the altitude graph. 

The data distribution plot in the time graph presents the same types of observa- 
tional data as the altitude graph. The plot corresponds to the 24-hour period 
covered in the OUE portion of the graph. 
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SECTION 4 - OPERATING INFORMATION AND SAMPLE JCL SETUP 

4. 1 OPERATING INFORMATION 

This section describes the minimum system configuration for the OCPLT4 
Program and gives timing estimates for program execution. 

4.1.1 System Configuration 

For the IBM System/360 Model 95 or Model 75, the minimum system configura- 
tion required to support the SD-4060 OCPLT4 Program consists of the follow- 
ing: 

• Three nine-track tape drives. 

• One seven-track tape drive. 

• Direct access space for an intermediate file. 

• Standard system input and output files. 

• The system data set for the SD-4060 named SYS2. SC4060 or 
SYS2. SD4060. 

• An SD-4060 plotter. 

4. 1.2 Timing 

A reasonable IBM System/360-95 timing estimate for OCPLT4 to process and 
plot a period of 1 month of data for 90 time graphs using a program load module 
is as follows: 

CPU = 3 minutes 
I/O = 15 minutes 

No timing estimate is needed for the SD-4060 plotter; however, turnaround is 
usually a few days. 
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4.2 JCL REQUIREMENTS 


Figure 4-1 shows the Job Control Language (JCL) required to execute OCPLT4 
using the program load module. 

Data set SC4060ZZ is the seven-track output instructions tape used for input to 
the SD-4060 plotter. Data set FT20 is allocated to the nine-track VC tapes. 
There can be as many as 24 of these VC files or tapes. Each file requires an 
FT20 card. Data set FT23F001 is a nine-track working-observations-file tape. 
For detailed information on these tapes, see Subsections 2.2 and 3.1. Data 
set FT22F001 is a required intermediate disk file, described in Subsection 2.3. 
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//USER JOB CARD 


//EXEC LOADER, REGI ON=390K, PARM= 'SIZE-390K' 

//GO. SYSLIB DD DSN = S YS 2 . S C4060 , D I S P = SH R ^ -ccd-vyvyv i 

//GO.SYSLIN DD DSN=OBJSET, UNI T=2400-9 , DI SP-COLD, PASS) , VOL SER , 

// DCB = (RECFM=FB,LRECL=80,BLKSIZE=3200), LABE L=( 1, BLP) 

//GO. FT06F001 DD DCB = B LKS I ZE=l4l , S PACE = C CY L, C 5 , D ) 

//GO. SC4060ZZ DD DSN=BURKE3UNIT=7TRACK, 

// DCB=(DEN=1,TRTCH=C,RECFM=F,BLKSIZE=1024) , 

// LABEL=Cl, BLP), DISP=CNEW, PASS), VOL=SER-BLA'NK 3 

//GO. FT20F001 DD UN I T= 2400-9 , VOL = SE R- XXXXX, LABEL- ( 1 , B LP), 

// DCB=CLRECL=133;RECFM=FBA;BLKSIZE=3325)/DISP-COLD;KEEP) ;DSN-PHI1 


/ /GO. FT22F001 DD UN I T = D I SK , DI S P= C NEW , PAS S ) , DCB- C B LKS I ZE- 22 , RE C FM F) , 

// DSN=SB3,SPACE = (CYL,C4,2)) „v, vvv s 

//GO FT23F001 DD UN I T = 2400-9 ,LABEL=( , B LP) , VOLUME- SE R- XXXXX , 

// DCB=CRECFM=VBS;LRECL=104;BLKSIZE=3436);DISP=COLD;PASS) 

//GO . DAT A5 DD :: 


NOTES: 


1. OCPLT4 system tape number 

2. Data set. name for output tape (user-specified) 

3. Output tape number (assigned by computer operator or user-specified) 

4. VC tape number (user-specified). There may be as many as 24 FT20F0xx 
cards in a single OCPLT4 run. These cards must be in ascending time 
order, one for each VC file. See Appendix A for samples showing how 
data set names and labels change for succeeding files. 

5. Working-observations-file tape number (user-specified) 


Figure 4-1. JCL Setup for Executing the OCPLT4 Program 
(SD-4060 Version) 
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SECTION 5 - PROGRAMMING METHOD AND SUBROUTINE DESCRIPTIONS 


5. 1 PROGRAMMING METHOD 

The first function performed by OCPLT4 is that of reading data cards and 
initializing variables for control of titling, grid generating, and grid labeling. 
Values read from input cards are carried into the grid drawing subroutine, 
TITLES. This information remains constant during execution of the entire pro- 
gram. Before the working-observations-file tape is processed for the data 
distribution portion of the plots, the tape is first time-sorted and rewritten on 
disk. This must be done in case the data on the working-observations-file tape 
was sorted concentrically. 

During execution of the program for altitude graphs (see Figures 3-1 and 3-2), 
each P-A and A-P period is determined. Component values from the VC report 
are selected for plotting when the radial distance has changed by at least AR km 
for the previous value, where AR equals the quantity (A-P)/100. 

Once a period has been completed, the subroutine DATAPT is called to plot 
data distribution within the time span of the period. The subroutine TIMTCK is 
also called to develop a time scale along the component curve to allow corre- 
lation between time and altitude. Time values are plotted on altitude graphs 
at 5-hour intervals for altitudes above a radial distance of 100, 000 km, and 
at 1-hour intervals for altitudes below a radial distance of 100,000 km. 

When time graphs are to be generated (see Figure 3-3), the three range differ- 
ence vector components are plotted against time (one day per graph). Sub- 
routine ALTCK is called to develop an altitude scale along the curve, for corre- 
lation with time. Altitude values to the nearest kilometer are plotted at 1-hour 
intervals along the time curves . 
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5.2 SUBROUTINE DESCRIPTIONS 


This section lists the subroutines the OCPLT4 Program uses from the SD-4060 
subroutine library, and describes the main routine (MAIN) and the calling 
sequences for the nine subroutines of the OCPLT4 source program. A listing 
of the source program appears in Appendix D. 

5.2.1 SD-4060 Subroutines Used 

The OCPLT4 Program uses the generalized subroutines for the SD-4060 (see 
Reference 1) to generate all plots. These subroutines do the plotting, generate 
the grids, and label the graphs and grids. The following is a list of the SD-4060 
subroutines used by OCPLT4: 

PAGEG XNORMZ NUMBRG EXITG 

LEGNDG YNORMZ LABELG 

GRIDG OBJCTG SETSMG 

SUBJEG MODESG LINESG 

5.2.2 OCPLT4 Source Program Subroutines 
5. 2. 2.1 MAIN Routine 

MAIN contains all the logic that controls the various options available to the 
user, computes all the coordinates for generating the OUE graphs, and also 
serves as the executive routine for all other subroutines. The following steps 
are the primary divisions of the MAIN routine: 

1. After initializing constants and flags that identify the options re- 
quested by the user, reading input cards and checking requests for 
consistency, and printing several messages to the user, the pro- 
gram will proceed, if no inconsistencies exist; otherwise the job is 
terminated. 
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2 . 


MAIN will rearrange the observations from the working observation 
file in time ascending order, if necessary. 

3. MAIN will call on A5READ to read the first (next) VC report and 
will select and restore the OUE values along with the corresponding 
time and range to be used in generating the OUE graphs. 

4. If altitude graphs are requested, MAIN will determine whether an 
A-P or P-A segment should be plotted next. Then the values of OUE 
components are checked by MAIN to ensure that they are within the 
requested limits. If a value is below the requested lower limit, it 
will be set to the lower limit and plotted. If it is above the upper 
limit, the value is divided by 10 and checked again. Should the new 
value exceed the upper limit, a message will be printed to that 
effect (see page 2-7). This process is continued until the entire 
graph for each OUE component is constructed from the information 
on one VC report. Similar activites take place when time plots are 
requested. 

5. After each OUE graph is constructed, the corresponding data dis- 
tribution plot is developed, if requested. 

6. All the plotting information and instructions to generate the OUE 
graphs for each VC report are stored on the output tape. When 
one VC report is finished, the next report is read and processing 
starts with step 3. This is repeated until all VC reports have been 
processed. Then the program terminates. 

5 . 2 . 2 . 2 Subroutine DATA PT 

This subroutine computes coordinates for and plots the data distribution between 
the times bounding each graph. 
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The calling sequence for subroutine DATAPT is: 

CALL DATAPT (ITME1, ITME2, ITME3, ITME4, XIX, INDTE) 

Argument Description 

ITME 1 Start YYMMDD 

ITME2 Start HHMM 

ITME 3 End YYMMDD 

End HHMM 

Location on page of left limit of data distribution graph com- 
puted in internal units used by the plotting routine 

INDTE Indicates whether it is an altitude or a time graph, and 

whether or not this pass-through requires reading of data 
tape or plotting of previously determined points: 

= 0, read and store data distribution points to be plotted for 
the altitude graph 

= 1, read and plot the stored points on the altitude graph 
= 3, same as 1, but for time graphs 
= 4, same as 0, but for time graphs 

5. 2. 2. 3 Subroutine TIMTCK 

TIMTCK plots the hours along the OUE curves for the altitude graphs. 

The calling sequence for subroutine TIMTCK is: 

CALL TIMTCK (JK, JNDTE) 

Argument Description 

Indicates number of values to be plotted 

Indicates type of component to be plotted: 

= 0, radial component 

= 1, component in orbital plane normal to the radial compo- 
nent 

= 2, component normal to the orbital plane 


JK 

JNDTE 


ITME4 

XIX 
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5. 2. 2. 4 Subroutine ALT CK 


ALTCK plots the satellite's geocentric distance along the OUE curves for the 
time graphs. 

The calling sequence for subroutine ALTCK is: 

CALL ALTCK (KJ, JNDTE) 

The ALTCK argument description is the same as for TIMTCK (\/ith KJ replacing 
JK). 


5 . 2 . 2 . 5 Subroutine TIT LE S 
TITLES plots and labels the graphs. 

The calling sequence for subroutine TITLES is: 

CALL TITLES (MTYPE, MSKIP) 

Argument Description . 

MTYPE Indicates the component to be plotted: 

= 1, radial component 

= 2, component in the orbital plane normal to the radial com- 
ponent 

= 3, component normal to the orbital plane 

MSKIP Indicates the part of the graph to be plotted or that cards are 

to be read: 

= 0, plot altitude graph from apogee to perigee 
= 1, plot time graph 

= 2, plot altitude graph from perigee to apogee 
= 5, read data cards 
= 6, label titles above graphs 

5. 2. 2. 6 Subroutine TAPES 

This subroutine reads and stores all VC tape numbers and the working- 
observations-file tape number, as well as the number of VC reports on each 
VC tape. This subroutine terminates program execution when all the VC re- 
ports on all input VC tapes have been processed. 
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The calling sequence for subroutine TAPES IS: 

CALL TAPES (IBLAP) 

Argument Description 

IBLAP Indicates whether input data cards 12 and 13 are to be read or 

whether to process the next VC report: 

= 0, read the next VC report on file; if none is available, 

read next file on tape; if none, request the next tape to be 
mounted on tape drive; if none, proceed to terminate job. 
= 5, read and store all tape numbers from cards 12 and 13. 

5. 2. 2. 7 Subroutine BSFTAP 

BSFTAP backspaces the current VC report on file to the beginning of that re- 
port, if necessary (i. e. , when both the altitude and the time graphs are to be 
plotted and the program has finished the altitude graph, the report is back- 
spacing to do the time graph). 

The calling sequence for subroutine BSFTAP is: 

CALL BSFTAP (NF) 

Argument Description 

NF FORTRAN file number of VC report to be backspaced 

5. 2. 2. 8 Subroutine TCONVO 

TCONVO, which was incorporated from DODS, converts times from DODS units 
to calendar units. 

The calling sequence for subroutine TCONVO is: 

CALL TCONVO (TIMDUT, IOUTIM, SEC) 

Argument Description 

TIMDUT Number of DODS units of time (DUT) from 0 h September 18, 

1957, to the calendar time 
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Argument 


Description 


IOUTIM The array containing the year, month, day, hour, and minute 

of calendar time 

SEC Seconds of minutes of calendar time (less than 1 minute) 

5. 2. 2. 9 Subroutine A 5READ 

A5READ reads data from the VC report and converts the components into a form 
useful for the main program. Conversion is done by separating the decimal and 
exponential portions of the components and of the range. 

The calling sequence for subroutine A5READ is: 


CALL A5READ (ITIME5, ITIME6, RADI, IEXP1, RAD2, IEXP2, RAD3, 
IEXP3, RANI, IEXP4, 13EOF) 


Argument 

Description 

ITIME5 

YYMMDD 

ITIME6 

HHMM 

RADI 

Decimal portion of the radial component iO < RADI < l) 

IEXP1 

Exponent associated with RADI 

RAD2 

Decimal portion of the component in the orbital plane normal to 
the radial component (0 < RAD2 < 1) 

IEXP2 

Exponent associated with RAD2 

RAD3 

Decimal portion of the component normal to the orbital plane 
(0 < RAD 3 < 1) 

IEXP3 

Exponent associated with RAD3 

RANI 

Decimal portion of the reference range vector (0 < RAN l < l) 

IEXP4 

Exponent associated with RANI 

I3EOF 

End-of-file indicator: 
= 1, end-of-file 
= 0, not end-of-file 


5.2.2.10 Subroutine B5READ 

B5READ reads UT from the working-observations-file tape and converts this 
time to calendar time. 
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The calling sequence for subroutine B5READ is: 


CALL B5READ (I3YMD, I3HM, I3TYP) 


Argument 


Description 


I3YMD 

I3HM 

I3TYP 


YYMMDD of observation 


HHMM of observation 

Type of observation: 

= 1, R range data 

^ ! minitrack direction cosines data 
= 3, ml 


= 9, R 
- 17, X} 
= 18, y! 


range-rate data 
RAO angle data 
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APPENDIX A - SAMPLE INPUT DECK SETUP 


The following list of cards is a sample input deck, including the JCL cards. 
The OUEs to be plotted are time graphs for the SSS-1 satellite for the time pe- 
riod August 11, 1972, to September 10, 1972. See Section 2.3 and Figure 4-1 
for a description of card images. This sample input was not used to obtain the 
sample output (Appendix B), but is presented to show the changes on FT20F0xx 
cards for multiple VC report tapes. 


// Z8NJBSSS JOB (G 10141841 E , P , G00080 , 005005 ) , 95 .QQG ,MSGLE VEL = I 1,1 ) 

// EXEC LOADER, R EG I0N= 390K, PARM= • S I Z E = 390K • 

//GO.SYSLIB 00 DSN=SYS2.SC4060,0!SP=SHR 

Z/G0.SYSL IN 00 0SN=0B J SET ,UN I T= 2400-9 , 0 1 SP= ( OLD , P ASS ) , V0L=SE R=1 24 1 M, 
// 0CB=(RECFM=FB,LRECL=80,BLK SIZE =3200) ,LABEL=(1 ,BLP) 

// GO • FT 06 F00 1 UD 0CB= BLKS IZ E= 141 , SPAC E = f C YL , ( 5 , 1 ) ) 

//GO.SC4060ZZ 00 DSN=BURKE»UNIT=7TRACK, 

// DCB=<0EN=l,TRTCH=C,RECFM=F,BLKSIZE=1024) , 

// L ABEL = ( 1, BLP), 0ISP=(NEW, PASS), VOL =SER=BLANK 
//G0.FT20F001 00 UN I T =290 0- 9 , V0L= SER = 2924P , L ABEL = ( l , BL P ) , 

// 0CB= (LRECL = 133, RECFM=F BA, BLKS IZE= 3325) ,D1 SP = ( OLD, KEEP) , 0 SN=PH 1 1 
// GO • FT 20 F002 00 UN I T = 2 AO 0-9 , V0L = REF = * . F T20F 00 1 , L ABE L= ( 2 , BLP ) , 

// DCB=(LRECL=133,RECFM=FBA,BLKSUE=3325) ,D! SP=( OLD , KE E P ) ,DSN=PH12 
/ / GO . FT 20 F00 3 00 UN I T = 2 AO 0-9 , V0L = REF = * . F T2 0F00 1 , L ABE L= ( 3 , BLP ) , 

// OCB= ( LRECL = 13 3,RECFM=FBA,BLKSI ZE = 3 325) , 0 1 SP = ( OLD , KE EP ) ,DSN=PHI3 
/ / GO . F T 20 F004 DO UNI T = 2 AO 0- 9, V0L = REF = * . F T20F 00 1 , LABE L= ( 4 , BL P ) , 

// DCB= (LRECL=133,RECFM=FBA,BLKSUE=3325) ,01 SP=1 OLD, KEEP) ,DSN=PHI4 
//GU.FT20F005 00 UN I T = 2 AO 0- 9 , V0L = RE F = * . FT 20F001 , L ABE L= ( 5 , BLP > , 

// 0CB=(LRECL=133,RECFM=FBA,BLKSUE=3325) ,01 SP=( OLD , KE E P ) ,DSN=PHI5 
//G0.FT20F006 OD UN I T = 2 400-9, V0L = REF = * . F T 20F 00 1 , L ABE L= ( 6 , BL P ) , 

// 0CB=(LRECL=133,RECFM=FBA,8LKSIZE=3325> ,01 SP=( OLD , KE EP ) ,DSN=PHI6 
/ / GO . FT 20 F007 00 UN I T = 2 AOO- 9, V0L = REF = * . F T20F 00 1 , L ABE L= ( 7 , BL P ) , 


A-l 


// DCB= <LRECL = 133,RECFm = FBA,BLKSI ZE = 3325) , Q I SP =( OLD , KE E P ) ,DSN=PHI7 
//GU.FT20F008 00 UN I T = 2400-9 , VOL = SER = 3 3976H , LABE L = ( 1 , BL P ) , 

// OCB = (LRECL=133,RECFM=F8A,BLKSI ZE=3325> , D 1 SP= < OLD , KE EP ) ,DSN=PHI8 
//GO.FT20F009 DO UN l T = 2 AO 0-9 , VOL = RE F = » . F T 20F 008, L ABE L= ( 2 , BLP > , 

// OCB=(LRECL=133,R£CFM=FBA,BLKSIZE=3325) , 0 I SP = < OL 0 , KEE P ) ,DSN=PHI9 
//G0.FT20F010 00 UN 1 T = 2 AO 0-9, VOL = R E F = * . F T 20F 008, L ABE L= ( 3 , BLP ) , 

// DCB= ( LRECL=133,RECFM=FBA,BLKSI ZE=3325) ,01 SP = ( ULO.KEEP) ,DSN = PHI 10 
//GO.FT20F011 DO UN 1 T = 2 AO 0- 9, VOL = RE F = * . F T20F008, L ABE L= ( 4 , BL P > , 

// DC8=(LR£CL=133,RECFM=FBA»8LKSIZE=3325),DISP=(0LD»KEEP) ,DSN=PHIU 
//G0.FT20F012 00 UN I T = 2 AO 0- 9 , VOL=RE F = * . FT 20F008 , L ABE L= ( 5 , BL P ) , 

// 0CB= <LRECL=133,R£CFM=F8A,8LKSIZE=3325> ,DISP = < OLD, KEEP) .DSN-PHI12 
//GO.FT20F013 00 UN I T = 2400-9 , VOL = RE F = * .F T 20F 008, L ABE L= ( 6 , BLP ) , 

// DCB= <LRECL = 13 3,RECFM = FBA,BLKSi ZE=3325) ,DI SP = ( OLO.KEEP) ,DSN=PHI 13 
//CO.FT20F014 00 UN I T = 240 0- 9, VOL = RE F = ♦ . F T 20F 008, L ABE L = ( 7 ,BLP ) , 

// DC8=(LRECL=133,R£CFM=FBA,BLKSIZE=3325),DISP=t OLO.KEEP) ,DSN=PH114 
// GO . FT 22 FOO 1 OD UN I T = D I SK , 0 1 SP = ( NE W , P ASS > ,DCB =( BLK S I Z E=22 , RE CFM=F ) , 
// 0SN=£83,SPACE=(CYL,[4,2)) 

/ / GO . FT 2 3F00 1 00 UN I T= 2400- 9, LABE L = I 1 , BLP ) , VOL =SER=2684P , 

// DCB=<KECFM=VBS,LRECL=104,BLKSIZE=3436) , D I SP= ( OLD , PASS ) 

//G0.0ATA5 DO * 


CARD 1 

c sss-: 

L 71961 

720919 

CARD 2 

O 

• 

o 

* 

o 

• 

o 

. 1 

CARD 3 

2.0 

10.0 

A .0 

CARD 4 




CARD 5 

.52 

.52 

.52 

CARD 6 

7208 1 1 

720910 


CARD 7 

0000 



CARD 8 

TIME 



CARD 9 

1 



CARD 10 

3 3.0 

33200. 6700. 

CARD 11 

60. 



CARD 12a 

292AP 

1A 


CARD 12b 
CARD 13 

33976H 

1A 


CARD 14 
CARD 15 

268AP 

01 



/* 
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APPENDIX B - SAMPLE OUTPUT 


The SD-4060 OCPLT4 Program output consists of two parts: the IBM System/ 
360 printer output, and the SD-4060 plotter output. 

The printer output provides the user with a description of the accomplished 
processing and reflects user input information, type of plot requested (time, 
altitude, or both), period to be plotted, plotting interval on graphs, input tape 
numbers, types and number of observations plotted in the data distribution 
box, and additional messages when appropriate or as requested by the debug 
option. 

Figure B-l is a sample printout for the SSS-1 satellite. Usually, only time 
plots are required for this satellite; however, in this run both time and alti- 
tude plots were requested. The definitions given below are numbered to 
correspond to the entries on Figure B-l. 

1. Displays part of the input parameters. (See description of input 
card images. ) 

2. Total period to be plotted — July 12, 1972, to July 30, 1972. 

3. Time and altitude plots are requested. 

4. Apogee and perigee heights as input by user. 

5. All three components of the OUEs are plotted (NCOMP). Interval 
between points on time graphs (TFREQ). Interval between points 
on altitude graphs (RFREQ). 

6. Interval between comparisons on VC report. 

7. VC tape number (contains 14 VC reports). 

8. Working-observations-file tape number (with one file). 

9. Start and end time of first VC report that is being processed (both 
at 0 hours). 
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10. Indicates that the data on working-observations-file tape from 
720710 0534 to 720710 2047 was rearranged into time ascending 
order. The data from 720710 2048 to 720812 1838 was in proper 
time order. The total time span of the observations on tape 2814H 
is from 720710 0534 to 720812 1838. 

11. The first portion of altitude graph to be plotted is an apogee-to- 
perigee pass. 

12. The period plotted for this portion of the graph is from 720712 0001 
to 720712 0349. The start time of the above period corresponds 

to the time of the first data point on the first VC report. However, 
it does not always correspond to an apogee or perigee point. 
Subsequent start times do correspond to proper labels. 

13. These dates correspond to dates in item 12. The first date 
appears on the left corner of the data distribution box. The second 
date, if different, appears as the last date in the plot (see Figures 
B-2 and B-8). 

14. The difference between the last and the first hour printed in the 
data distribution plots. 

15. Type and quantity of data plotted in the data distribution plot. 

This reflects the contents of the observations in the working- 
observations-file tape (tape number 2814H). However, "passes" 
refers to the number of asterisks in the data distribution plot. Only 
one observation per minute is plotted and counted in this number. 

16. Items 12 through 15 are repeated for each component of the OUE, 
if all three components are plotted. 

17. Items 12 through 16 are repeated for all A-P and P-A passes 
until the entire VC report is plotted for altitude plots. 
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18. Indicates completion of a VC report. 

19. Indicates the back spacing of the VC report to the beginning of 
the report when both altitude and time graphs are requested. No 
backspacing is needed when only one option is requested. 

20. Indicates that time graphs are being prepared by the program. 

21. IDATE is the start date; ITIME6 is the start hour of the data on 
the time graphs. The data distribution box information, similar 
to these described for the altitude graphs (items 12 through 16) 
is repeated. This time the span of the graphs and plots is 24 
hours. 

22. When an entire VC report has been processed for both options, if 
needed, the next VC report is called in and the start and end time 
of the VC report is printed as is the information from 9 through 
21. This is repeated until all VC reports on all VC tapes are 
processed. 

Because SSS-1 completes approximately three orbits in a 24-hour period, 
there are six sets of altitude plots, two for each orbit (see Figures B-2 
through B-7) , and one set of time plots (see Figures B-8 through B-10). 

Each set consists of the three OUE components. 

The plotter output provides the graphic display of the OUE component (see 
Section 3. 3). Figure B-2 through B-10 show a typical set of altitude and time 
graphs. The altitude graphs consist of two sets: apogee-perigee and perigee- 
apogee. 
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• tt-9 
72C710 < C £ 

720710 -tit 
720710 1242 
-■72 0710 124 2- 
720710 1242 
-TKTtO 1 24 2 
72C71C 124fc 
-72C7iO--t24 1 
72C710 1247 


2 

-3- 

2 


-3- 


2 

3 

2 
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Figure B-l. Sample Printout for the SSS-1 Satellite (1 of 2) 
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Figure B-l. 


Sample Printout for the SSS-1 Satellite (2 of 2) 
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xx -jzwzotxoo r» — 


RUN DATK 720912 
( 71961 > 


MISSION AND TRAJECTORY ANALYSIS DIVISION, GODDARD SPACE FLIGHT CENTER 

ORBITAL UNCERTAINTY ESTIMATE FOR SSS-1 



DATA DISTRIBUTION 


RADIAL DISTANCE FROM CENTER OF EARTH * 1000 KILOMETERS 



Figure B-2. Altitude Graph for SSS-1 Satellite, 
Apogee-Perigee Radial Component 
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Figure B-3. Altitude Graph for SSS-1 Satellite, Apogee-Perigee 

Component in Orbital Plane Normal to Radial Component 
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Figure B-4. Altitude Graph for SSS-1 Satellite, Apogee-Perigee 
Component Normal to Orbital Plane 
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MISSION AND TRAJECTORY ANALYSIS DIVISION, OODDARD SHOE FLIGHT CENTER RUN DATE 720912 

ORBITAL UNCERTAINTY ESTIMATE FOR SSS-1 < n»«i * 



Figure B-5. Altitude Graph for SSS-1 Satellite, 
Perigee-Apogee Radial Component 
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Figure B-6. Altitude Graph for SSS-1 Satellite, Perigee-Apogee 

Component in Orbital Plane Normal to Radial Component 
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MISSION AND TRAJECTORY ANALYSIS DIVISION, GODDARD SPACE FLIGHT CENTER RUN DATE 720912 

ORBITAL UNCERTAINTY ESTIMATE FOR SSS-1 < mei > 


Sm 9iS iEi mi ^51^5^5^51 

■■HU 



DATA DISTRIBUTION 


10 . 20 
RADIAL DISTANCE FROM CENTER OF EARTH * 1000 KILOMETERS 


NO DATA FOR THIS PERIOD 
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Figure B-7. Altitude Graph for SSS-1 Satellite, Perigee-Apogee 
Component Normal to Orbital Plane 
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MISSION AND TRAJECTORY ANALYSIS DIVISION, GODDARD SPACE FLlGVfT CENTER RUN DATE 720912 
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Figure B-8. 


Time Graph for SSS-1 Satellite, Radial Component 
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Figure B-9. Time Graph for SSS-1 Satellite, Component in 
Orbital Plane Normal to Radial Component 
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Figure B-10. Time Graph for SSS-1 Satellite, 

Component Normal to Orbital Plane 
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APPENDIX C - INTEGRATED GRAPHICS SOFTWARE (IGS) 

ERROR CODE 


TABLE 3-4 


IQS ERROR CODE 

THIS TABLE DESCRIBES THE MEAN I HO OP BACH POSSIBLE IOS ERROR 
MESSAOB. WHEN AN ERROR OCCURS, SUBROUTINE BRRZZ IS CALLED TO PRINT 
OUT AN ERROR NBSSAOB . THE ERROR NBSSAOB WILL RBAD AS FOLLOWS 
YOU HAVE COMNITTBD ERROR NO. ’NO* DU R I NQ THE PLOTTINQ OF FRANB NO. XX 
THE ERROR VALUE WAS VALUEll), VALUE(P), VALUE ( A ) . 

CONTROL IS RETURNED AFTER THE NBSSAOB IS PRINTED--THE JOB IS NOT 

TERN I NATBD . 


NO SUBROUTINE VALUE DESCRIPTION 


1 

OBTSNO 

NO 

ILLBOAL NODE SET NUMBER IN CALL. 

2 

LEON DO 

N 

ILLBOAL CHARACTER COUNT IN CALL. 

3 

L I NBSO 

NO 

ILLBOAL NUMBER IN CALL. 

4 

MODS SO 

1 TAPE 

ILLBOAL TAPB NO. IN CALL. » 

5 

NUNBRQ 

PMT 

ILLEGAL PORMAT IN CALL. 

6 

OBJCTO 

• - 

MAX X 0| Y LB MIN X OR Y IN CALL. 

7 

P AO BQ 

- 

ILLBOAL AROUMBNTS IN CALL. 

e 

POINTO 

N 

ILLBOAL NUMBER IM CALL. 

9 

SBQNTQ 

N 

ILLBOAL NUMBER IN CALL. 

10 

TAB SO 

N 

ILLBOAL NUMBER IN CALL. 

1 3 

LABBLO 

- 

ILLBOAL PORMAT IN CALL. 

1 2 

0 R I DO 

- 

OR ID TOO SMALL TO DRAW. 

1 1 

MLTPLO 

NL1NBS 

ILLBOAL NUMBER IN CALL. 

14 

TITLBO 

- 

ILLBOAL AROUMBNTS IN CALL. 

u 

SBTUPQ 

- 

ILLBOAL AROUMBNTS IN CALL. 

i| 

SUBJBO 

- 

MAX X OR ,Y BQ MIN X OR Y . 

i| 

LABBLO 

- 

ILLBOAL AROUMBNTS IN CALL. 

i! 

LABBLO 

- 

GRID TOO SMALL TO LABBL 

1 9 

LABBLO 

- 

ZERO SUBJECT SPAQfe 

20 

OR IDO 

- 

ILLBOAL AROUMBNTS IN CALL 

21 

SBTUPQ 

- 

NOT ENOUGH ROOM TO DRAW A ORID. 

22 

SBTUPO 

- 

DENSITY LB 0. 

2 3 

SBTSMO 

N 

ILLBOAL MODB SBT NO IN CALL. 

24 

SBTUPO 

- 

ORID WILL NOT PIT ON PAQB 

25 

TBXTO 

N 

ILLBOAL CHARACTER COUNT IN CALL. 

26 

LABBLO 

- 

LABELS WILL NOT PIT ON PAOB 

27 

PACKZZ 

- 

NO INITIALIZATION CALL TO MODBSO 

28 

ORAPHO 

N 

ILLBOAL AROUMBNT IN CALL. 

29 

SUBJBO 

- 

MINUS VALUB FOR LOO ORID 

30 

SBTUPO 

- 

TOO MANY CYCLES IN LOO ORID. 

31 

VBCTZZ 

- 

NO VECTOR C HA RAC TB R PONT INITIALIZED 

32 

SCALZZ 

X 

BAD X-COORDINATB . 

34 

VBCTZZ 

CHAR 

CHARACTER NOT IN PONT. 

35 

TITLBO 

- 

NOT BNOUQH ROOM TO TITLE ORID. 

38 

VBCTZZ 

CHAR 

RBQUBST FOR NON-EXISTENT VECTOR CHAR 

37 

NVBCZ 

HP 

BYTE NUMBER IS ZERO OR NBOAT1VB. 

100 

PSUBJO 


MAX. THETA EQUAL TO M I N THETA 
AND/OR MAX. RADIUS EQUAL TO MIN. 
RADIUS . 

33 

SCALZZ 

Y 

BAD Y -C OORD I MATE . 

101 

PSUBJO 

- 

NIM. RADIUS ORBATER THAN MAX RADIUS. 

102 

POBAPO 

N 

ILLBOAL AROUMBNT IN CALL. 

10 3 

VBCAZZ 

- 

MORE THAN 380. DEO. OP CHARACTERS. 

104 

VBCAZZ 

N 

ILLEGAL AROUMBNT IN CALL. 

105 

PLINBO 

N 

ILLBOAL AROUMBNT IN CALL. 

108 

POLPTO 

N 

ILLBOAL AROUMBNT IN CALL. 

107 

PS BO MQ 

H 

ILLBOAL AROUMBNT IN CALL. 

108 

PNLTLO 

N 

ILLBOAL AROUMBNT IN CALL. 

109 

POLBO 

N 

ILLBOAL AROUMBNT IM CALL. 

1 1 0 

PVBSSZ 

N 

ILLBOAL AROUMBNT IN CALL. 

1 1 1 

PLABLO 

- 

ILLBOAL PORMAT IN CALL. 

1 1 2 

PLABLO 

- 

ILLBOAL AXIS. 

1 1 3 

PLABLO 

- 

ZBRO SUBJBCT SPACB . 

1 14 

PVBAZZ 

N 

ILLBOAL AROUMBNT IN CALL. 


C-l 




TAILS J-4 


(COH*T) 


1 1 5 

6 0X0 

N 

1 1 e 

CLASBQ 

LBVBL 

202 

SUBJ 3D 

ZM 1 N 

20 1 

SUB J 3D 

YM I N 

200 

SUBJ 3D 

XM I N 

1 1 7 

PSBTQ 

MODS 

20 3 

PLOTS 3 

X 

204 

PLOTS 3 

X 

205 

PLOTS 3 

Y 

206 

PLOTS 3 

Y 

20 7 

PLOTS 3 

Z 

2 0 H 

PLOTS 3 

Z 


ILLBOAL ARGUMENT IN CALL 
ILLBQAL ARQUMBNT IN CALL 
MINIMUM Z = MAXIMUM Z. 
MINIMUM Y = MAXIMUM Y 
MINI MUM X = MAX IMUM X 
ILLBQAL ARQUMBNT IN CALL 
X MAX I MUM X 

X MINI MUM X . 

Y MAX IMUM Y . 

Y MINI MUM Y 

Z MAXIMUM Z 

Z MINIMUM Z 


02 



APPENDIX D - SD-4060 0CPLT4 SOURCE PROGRAM 
COMPILATIONS LISTING S 

This appendix presents a compilations listings of the SD-4060 OCPLT4 source 
program. The subroutines are listed as follows: 


Figure 

Title 

Page 

D-l 

JCL Used in Compilation of 
SD-4060 OCPLT4 Program 

D-2 

D-2 

MAIN Routine 

D-3 

D-3 

Subroutine DAT APT 

D-16 

D-4 

Subroutine TIMTCK 

D-22 

D-5 

Subroutine ALTCK 

D-23 

D-6 

Subroutine TITLES 

D-24 

D-7 

Subroutine TAPES 

D-27 

D-8 

Subroutine BSFTAP 

D-29 

D-9 

Subroutine T CON VO 

D-30 

D-10 

Subroutine A5READ 

D-33 

D-ll 

Subroutine B5READ 

D-34 


D-l 


// 7 *K ,>PP| T joa <GI 014 ie 41 E-.E.CCCC 30 t C 0 1 C 01 > .QQO . MSGLEVEL= U . 1 > 

// exec FORTRANH,PARM=« HAP , IC.OPT= 2 • ,REGION= 50 GK 

XXOEFAULT PROC FORTH AN= I EX AA 0 C • NBL K = 4G OOGOGIOC 

xxsauRce exec pgm=£fortran ,regich= 30 ok 00000200 

TEF&aa-i suBST-iTuTHCN aci- - pgm=x e k a a oc ♦REoicN=aaoK 

//SOURCE.SYSL IN DO CSN = FKS.CMT = <240 0 -S,,DEFER>.DI SP=(NEW .PASS) . 

// CCe=<RECFM=FB.LRECL=£C.eLKEI2E=32C0,DEN=3) . 

// L A B EL= ( 1 * BLP • .OCT ) #VCL = SEH = ScAN 

X/SY 5 LIN DC DSN=fc&aBJMaC . SPACE = < 3200 . (CN 3 LK# 1 0 ) .. .ROUND ) ,UNIT=DISK. 00000300 
IEF 652 I SUBSTITUTION JCL - D SN =£ 00 E JMO C . SP ACE= ( 32 C 0 . ( 40 . 1 0 > . . . R CUNC ) . UN IT = D ISK , 
XX BISP=<-MOG*PAS.S ) ,OCa=<-R£CFM=Ft(, LRECL^flO .BLKS I ZE = 320 0 ) OOOOC 4 CC 

//SCURCE . SYSPR INT CC SP AC E= ( C Y L , < S . 1 > > 

X/SYSPPINT CC SYSOCT=A . C C E= < R ECF M = VB A ,LRECL = 1 37, DLKSI ZE=7265 ) , 0CCCC5J0 

XX UN I T = < D I SK . £EP=SYSLI N ) 0 OCOOcOO 

XXSYSPUNCM DD SYSOCI=8,OCE=(HECFF=Fe.LRECU=80,BEKSIZE=723C) 00000700 

/✓SOURCE .SYSUOUMP DC SYSQCT=A 

X /SY S-U DU MP DO SYS0C t=A. SPACE =( TRK ,-SJ UOC0030C 

XXSYSUT1 CC SPACE = < TRK , <C .S )» . UNIT = (DI SK, SEP=( S YSLl N ,5 YSPR INT > ) O00C09OC 

XXSYS0T2 CC SPACE = < C YL . < 1 . 1) ) . OCOOIOOC 

XX UN IT = (D ISK ,SEP= (SYSLI N, SYSPRI NT. SYSUT1 ) > OOOOllCC 

//SOUKCE.SYSIN CC * 

// 



ALLOC 

. FCH ZONJEFtT 

SCL«CE 

I EF237 I 

0 C 3 

ALLOCATED 

TO 

SYSLIN 

I EF237 l 

3 3 1 

ALLOCATED 

TO 

SYSPR I T 

I EF 23 71 

332 

ALLOC ATEC 

TO 

SYSPUNCH 

I EF237 I 

333 

ALLOCATED 

TC 

SYS UD UVF 

I EF 2 371 

334 

ALLOCATED 

TO 

STSUT 1 

IEF2a7I 

232 

ALLOC ATEC 

TO 

S VSUT2 

I EF 237 I 

223 

ALLOC ATED 

TO 

S YSIN 


Figure D-l. JCL Used in Compilation of SD-4060 OCPLT4 Program 
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LEVCt 20^4 <AUG ?l-t OS/ 360 F O R T R AN H 

G O MP ILE R O P TI ONS — NAME- MAI N , GPT=0 2 .L I NECNT =50»S I2E-OOOOK , 

SOURCE. EECOICt NOLIST, NO DECK ,LOAO , MA P , NOEDIT # I D. XREF 


FSN 

0002 


10)4ERRl>R(50> , ERROR 1 ( 50 ) . ERR0R2< 50 ) • R ANG E7 ( 30 ) , I HOUR 2 ( 30 ) , ABSIC(JO) 

0008 




I .ASS IC1 < 30 ) i A SSI 02 < 30 ) » I I OAT » I IDAT 1 . IH 

0009 

ISN 

0003 


CQMMC.4 EfiRLJl , E RRLC2 . ERRL03 ♦ ERRP I 1* ERRH I 2 * ERP HI 3 


rstt 

000 4 


COMMON AMODE< 200) . "CON* MANY »LOG 


I SN 

0005 


COM WEN NiiS 1 .NSS2 ,N5SJ. NSS4.NSS5.NSS6.A5 . A7.B3.B6. IX IY 


r sr* 

0000 


OAT A INCRP* J5YOUT/3 ,0/ 


ISN 

0 0 C 7 


DATA NO ST OP/ 4H A b CD/ 


rsN 

0008 


DATA CCNN7THC/ 

0012 

ISN 

0009 


DATA S0RT3T/4H / 

0013 

rsN 

0010 


DATA TTMLY/4HTIME/ 

0017 

ISN 

rsn 

00 1 1 
0012 


DATA 3TRT 1/4HYYMM/ 

0021 


DIMENSION a 2 ( 120) i X 3 < 120) . X4{ t 2 C ) »X5( 120) t X 6 { 120) i X 7{ 120) *KO( 120) • 





1 X9(120).V2(120),Y5(I20).Y4(120).Y5(120).Y6(120).Y7(120).YO(120). 





2 Y9“t"120 ) 


I SN 

00 13 


DIMENSION XI { 1 20 ) . Y 1 l 1 20 ) 


ISN 

O0T4 


DIMENSION XMR31 1 500 .RANGE! 1500 ) .RANGE l < 1500) .PANGE21 1500 ) 


ISN 

0015 


DIMENSION SAVE1 < 120 ) »SAVt2( 120) *SAVE3< 120 ) »SAVE4( 120) 


rsn 

0010 


01 MENS I0M r^UTTM-ttr) 


ISN 

0017 


REAL48 TIMDUT .TIMtND 


ISN 

0018 


TNTEGER A5i A7. 03* B5 


ISN 

0019 


INTEGER UPPER. 1 . UPP E R2 . UPPER 5 . B I G 1 , B I G2 . B I G 3 . TOP 1 .T0P2.T0P3 


I5N 

0020 


WRITEC J5YOUT .2) 


I SN 

0021 

2 

FCRMATt 1H 1 ) 


I SN 

0022 


A5 *= ZD 


ISN 

0023 


A7 = 21 


tsn 

0024 


S3 = 22 


ISN 

0025 


05 = 23- 


rsN 

0026 


STD — o , 


ISN 

0027 


XRAN 1=0 • 


rsTj 

0028 


rx-rrsB 


ISN 

0029 


MTYRE=1 


rsis 

0030 


J INX=0 

0029 

ISN 

0031 


NOGO=0 

0030 

rsn 

00*32 


FSKIP=5 

0033 

ISN 

0033 


CALL MQUESGt AMODE, 0 ) 


I SN 

0030 


CALL SET 0MQ t AMODE i 19 1 T 360 * ♦ 


ISN 

0035 


CAUL 5ETSMG< AMOCE.20, 1 023.) 


rsN 

0036 


CAUL TITLES ( MTTPE . WSlCl P > 

0032 

ISN 

0037 


ITAPE= 1 

0033 

I SN 

0000 


I60P 1 = 1 

0034 

ISN 

0039 


JK=*0 

0035 

HIM 

0040 


PSKtP^-0 

0036 

I SN 

0041 


I ADD =0 

0037 

rsw 

0042 


I SET -0 

0030 

ISN 

0043 

1 

X SET 1=1 

0039 

ISN 

0044 


IT6ffH=0 

0040 

ISN 

0045 


I GQP=0 

0041 

i sn 

0046 


1*0 

0-0*42 

ISN 

0047 

res 

READ 102.TIMEY 

0044 

ISN 

004“E 

FORMAT fAAl 

0043 

ISN 

0049 


READ 3 0 50 » NSS 1 1 NSS2 * NS S3 *NSS4 *NSS5»NSS6 


raw 

0030 

0000 

PGIW4AT { 6 1 1 ) 


ISN 

0051 


IFCTIMEY.EU.SORTST. AND. NSS6.EG. 0 ) PRINT 9101 


I SN 

00 53 

9161 

PDRMATt OX » ' T IME AND ALT I T UDt * ) 


ISN 

0054 


lF(TIMEY.tU. TIMLY. AND. NSS6. E 0 .0 > PRINT 9102.TIMEY 


raw 

0000 

91 02 

fCRMATtOX. A-4-1 


ISN 

0 057 


IFtTIMEY.tuJ. JCJRTST. AND.NSS6.EC. 1 ) PRINT 9103 


ISN 

0059 

91 03 

FORMAT (bX , * ALT I 7 UDE * ) 


ISN 

C 060 


IFTTIMEY.EU.TIML Y. AND. NS36.EQ. 1 ) PRINT 9104 


ISN 

0002 

9104 

FCRMATtoX* ♦CONFLICTING USER I NPUT* A* A* CHECK DATA CARDS 8 AND 9'l 


ISN 

0063 


READU NCRU, 1 03) NCOMP* TFREU . A PO G EE * PER IGE 


i-sn 

0004 

ro-3 

F CM NAT (-ITfTXlF3»T,-iX »FDrO» 1 X*F8r0) 


ISN 

0065 


TFREC=TFF«eu/oO. 


i-sn 

0000 


XRANGt=< APcufcE-PEKTGE)/ 1000. 


ISN 

0 C 6 7 


RFREU=XRANGt/ 100 • 


ISN 

0008 


imiTEt J3YDUT • 10c) APOGEE ♦ PER I GE 


ISN 

0069 

106 

FORMAT { IH , • APCGEE= * ,F1 1 .0 . • PER I GEE= • *F 1 1 .0 ) 


rsN 

007 0 




ISN 

007 1 

1C5 

FORMAT < 1 PQ *oMNCUMP = * 1 2 » 4X . 6HTFREU= » F8 . 5 • 3HHR S « 4 X • 6H PFRE Q = , F 8. 4 • 



I UH KM/1C00) 

C *READ ANY OPTIGNAL INPUT TIME DIFFER*; NC E 
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D-3 



I~5N 

0072 


READ 3052. T3DIFF 


I SN 

0073 

3052 

FCRMAT ( F4 . 0 ) 


rsrj 

ctcta 


IX (T3DIFF.NE .0. > T"3DTFF =T 3CI FF/0C4 * 


ISN 

0076 


IF (T3DIFF.EQ.0. ) T 3D IFF = o0./ec4. 


tsk 

O'OTB 


33DIFF=( (T3DlFF*864.)/60.)+.l 


ISN 

0079 


IDIFF= IF IX(S JDIFF) 


ISN 

0X90 


PRINT 1311.T3DIFF 


ISN 

ooai 

101 1 

FORMAT ( In . • T 3 D I F F = • ,F 10 .5t 1 X . * CENT I DA Y 3 * > 


rsra 

0TTET2 


TFTNSob.EQ.UI PRINT 3. ID IFF 


ISN 

0084 

3 

FOR MAT ( 1 H ,* IDIFF=* . 13. IX, ‘MINUTES =T3DIFF*> 


rsN 

oce 5 


I BL A P= 3 

0046 

ISN 

0 0 do 


CALL T APES C I OL AH 1 

0047 

rsN 

Oca/ 


IF C7SET.NE.0I GC TO 25 

0049 



C * 

THt PURPOSE JF THIS SEGMENT IS TO DETERMINE END TIME AND 




c * 

STATIC TYPE 7AOJES 


ISN 

0089 

14 

CONTINUE 


rstf 

0090 

2 1X0 

READ ( B3* END = 3150) 


ISN 

0091 


GC TO 3100 


rsN 

0W2 

3T50 

BACKSPACE B5 


ISN 

0099 


6 AOK SPACE 03 


rsn 

3094 


CATTL CT5YEAD ( JTME , JTMtrl. ITTPEdt) 


I3N 

0C95 


REWIND 85 


rsw 

0096 


TMNS55.EQ.0 J PRINT 332 . JTME . JTME 1 


ISN 

0098 

332 

FORM AT ( 1 H ,‘JTME = • . I 6 • 2X, • JTME 1 = M6) 


rsN 

0 C 99 

Zb 

HEADCA5 . 1 04 . END“26l STRT2 


ISN 

0 100 

26 

IFCTIMLY.NE.STRT2) GO TO 25 



0 102 


reap rxa , ron rn.rr 2 .tro.iT 4 


I SN 

0103 

1 07 

FORMAT (40X. I 6. 1X.I4.22X. 16. 1 X . I 4 ) 


rsN 

0104 


T ImuUT — 0 • 


ISN 

0105 


IOmTlM(l) = I T 3/ 10 000 


rsn 

<rroo 


tCtinMt-21 = (IT3- - 1 00004 IOUT IM( 1 ))/ 100 


ISN 

0107 


I0UTIM13) = IT3 - ( 1 OOOOMOUT IM( 1 ) + 100*1 OUT IM ( 2) ) 


rsN 

010 6 


i cut t n m = m/i-co 


I SN 

0 1 09 


1 OUT IM ( 5 ) = I T 4 - 1 004I0UT IM( 4 ) 


rsN 

OTTO 


CA1_L "TC'CNVO tTtMOUTi JOUTIM* SEC) 


ISN 

0 111 


TIMEND = TIMOUT 


rsN 

01 1 2 


PRtTtT 1071 , T 30 IFF 


ISN 

0 113 


PRINT 108.IT1 .IT2. IT3. IT4 


ta* 

01U 

lOd 

FORM AT ( 1 rt . 1 5X . * 3T ART T I ME = ‘ .I6.1X.I4.11X.* END T I ME =*.I6»1X.I4) 


ISN 

0 115 . 


IF ( ISt T .EUi 1 ) GO TO 31 

0074 



C CHECK F"OR TIME PEXTCO COVET? E D ON VECTOR CtTWP AQITON TAPE 

7075 

ISN 

0 117 

23 

I =1 + 1 

0082 

rsN 

07 re 


IF CC3N.NE.CCNN ) GO TO 27 

0083 

ISN 

0120 

232 

CALL B5READ lITZMLd(X), ITIME9(I). ITYPE(I)) 


rsN 

17121 


IF (l.EU.l) GO TO 23 

0086 

I SN 

0123 


IF ( I.lT.7) GO TO 233 

0087 

raw 

0129 


EH 1-1 ) . AND • •£ T-J ME 9 < I ).E<i. I T IME9< I- 1->»4ND.-I TYP 

0088 




1E(I) .t J.ITYPEt 1-1 ) ) GO TO 232 

0089 

rsN 

0 127 


IF t IT I ML 0 t I ) .EC.ITIMEOt 1-2 ) . ANC. ITIME9< I l.EO.I TIME 9(1-2) . AND. I TYP 

0090 




1E( I ) .EO.ITYPEt 1-2) ) GO TO 232 

0091 

rsN 

0129 


IF ( IT I Mta ( I ) .ECJ. ITIMEdt 1-3) .AND. IT IME9( I ) . EQ . I T I ME 9{ 1-3) . AND. I T YP 

0092 




1E(IJ .= U . I TYPE ( I -3 ) ) GO TO 232 

0093 

rstr 

rm 


if ( iT-trtdt itve-c* it i Meat i-4).and« iTiMe^t ii.fcu. i t i-me«< 1 - 4 ) .and* i typ 

0094- 




lttl ) *bJ. ITYPEt 1-4) ) GO TC 232 

0095 

rsN 

07-33 


IF ( IT I MEU( I ) .EC.I TI ME8( 1-5 J .AMU. I T IM£9( 1 ) • EQ • I T I ME 9 < I -5 ) • AND* I TYP 

0096 




1E< I ) .EG. ITYPE ( 1-3) ) GO TO 232 

0097 

rsN 

Or35 


IF ( ITIMEdt I ) *E C » I T X ME 8 ( I-6).ANU. IT IM£9( I ) • EQ . 1 T I ME 9( 1-6 ) • AND* I TYP 

0098 




lEtll.EJ. I TYPE ( I -6 ) ) ^ O TO 232 

0099 

rsn 

0T37 

233 

IF rm^EBri7*ET.ITIME8I 1-1 ) ) GC TO 23 

oroo 

ISN 

0139 


IF ( ITIMEdt I ).GT. ITIMEdt 1-1 ) ) GC TC 24 

0101 

rsN 

01 4 I 


IF < IT I ME9 ( I ).LE.ITI ME9( I- 1 ) ) GC TO 23 

0102 

ISN 

0143 

24 

1=1-1 

0103 

I SN 

0144 


IF ( I.EiQ.O ) GO TO 27 

0104 

I SN 

0 1 46 


WRITE ( 83,202) NCSTOP, ITIMEdt I) . ITIME9( I) .ITYPEt I) 


rsn 

0 14/ 


PRIITT 202 , NU3TOP . IT7 PE B t I ) » I TI ME9T I ) . I TYPE t I ) 


ISN 

0148 


GO TO 24 

0 106 

rsrr 

0X49 

27 

CALL 85 READ (ITIVEB(l), ITIME9(1), ITYPEt 1 ) I 


ISN 

0 150 


IF ( ITI M£8( 1 ) .EQ. JTME) GO TO 17 

0 1 08 

rsH 

0152 


GTTTC 2 9 

0109 

ISN 

0 153 

17 

IF ( ITIME9C1 ).EG. JTME1 ) GO TO 28 

0110 

I SN 

0 155 

29 

TOTTt C33. Z 0 X 1 ITCirrOP.rrrREBt 1 ) ,TTIWE9t 11 .ITTPET 7) 


ISN 

0156 


GO TO 2/ 

0 112 

raw 

0 15/ 

ZB 

WRTTE CB3.20Z) SORT ST , IT I ME8 ( 1) , I T I ME9 ( 1 ) , X TYPE t l > 


I SN 

0159 


PRINT 202. SOFTST. ITIMEdt l ) , ITIME9( 1) , ITYPEt 1 ) 


rsrj 

0 159 

202 

FORMAT (IXi A4. 3X, Id, IX, 14, IX, 12) 


ISN 

0160 


END FILE 03 
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I SN 

0TB1 


REWIND B5 



I SN 

0162 


IF ( T I ML Y *EU. T I MEY ) GO TO 1004 



isn 

0IO4 

31 

CCNTINUE 



ISN 

0 165 


PRI.NT lOU’TJDlFF 



rsrj 

crt 66 

rtT4 

FORMAT T7X, A4) 



I SN 

0167 


IF ( XTGPH.EU. 1) GO TO 749 


0124 

rsw 

0109 


HT ( IScT , trtJ 4 Ht 80 TO 39 



I SN 

0171 


l START = | 


0126 

rsN 

C7TT2 


I-PTrNSSO.EO.O ) PRINT 211. I ST ART 



ISN 

0 174 

2000 

IF I TI MLY.EO.T IMEY) GO TO 25 


0127 

rsu 

0*^6 

3-9 

eeNT-t Ntl£ 



ISN 

0 1 77 

esy 

READ (A5.890) CHC 



I-SN 


t 

FORMAT 4-+*-»-A+4 



ISN 

0179 


IF (CHC .NE.STRT 1 ) GO TO 889 



rsN 

OT8-1 

40 

CENT I. TOE 



ISN 

0182 


L ARGE1 =0 



I-SK 

01 83 


EAR6E2~G 



ISN 

0184 


LARGE3=0 



ISN 

0185 


CALL A5REAO (ITIME5. ITIME6. RADI. IEXP1, RAD 2 . 

IEXP2. RADJ. 





1IEXP3. RANI. IEXP4, I3E0F) 



ISN 

o ias 


IF (-1 3£OF.£U»l ) GC TO 1 OCA 





C READ IN LINE OF VECTOR CCMPARISCN T 4 PE 


0134 

ISN 

0 188 


T IMDL7 = 0 . 



ISN 

0 189 


ICUTIM(l) s ITIME5/10000 



I SN 

0190 


IOUT 1 M ( £ ) = ( I T I ME E - 1 0 00 0 * I ULT I M { 1 ) )/ 1 0 0 



ISN 

01S1 


ICL T IV ( 5 ) = I TIMES - ( IOCCO*IOUTIM( 1 ) ♦ 100* IOUT I M( 2 ) ) 


FSN 

©1-92 


ICVI-1N44 -* = f TIME-6/4 00 



ISN 

0 193 


IOLTIMC5) = IUME6 1 00*1 OLTI M(4 ) 



ISN 

0 194 


CALL 7CONVO (TIMOUT, IOUTIM. SEC) 



ISN 

0 19 5 


TIMOUT = TIMDLT ♦ T3DIFF 



ISN 

0 196 


IF ( 7 I MDL T. GE. TI MENO) GO TO 1004 





C CHECK FCR END OF VECTOR COMPARISON TAPE 


C 1 37 

ISN 

04-98 

6 C 

UADI-RAD 1* U.C**IEXP1 


0 138 

ISN 

0 199 


RAD2=RAD2* 1 0. 0**IEXP2 


C 1 39 

ISN 

0200 


R AC3-=R AO J* 10 .0*» IE XP3 


0140 

ISN 

0 20 1 


RANI =( RANI* 10 .0* *IE XP4) /l COO. 0 


C 1 4 1 



C CCNVEFSGN OF NLM 6ERS TO REAL 


0142 

I SN 

0202 


IF( RAD 1 .LT. ERRLU1 ) RADI =ERRLO 1 



ISN 

0204 


I F ( RAD 2.LT • CRRL02 I RAD2=ERRL02 



ISN 

C2 C‘6 


I F (R AD3 *LT .SMRLC3 ) RAD3=ERRL03 



ISN 

0 2 C 8 


I F( RAD 1 .GT .£RRHI 1 ) LARGE! =1 



ISN 

02 10 


I F{ R AC 1 . GT .ERKHI 1) R AD l = RAD 1 / l 0 . 



ISN 

02 12 


IF( RAU2.GT. ERR HI 2) LARGE2=1 



I SN 

C 2 1 4 


IF( HAD2 .GT.ERRHI 2 ) R AD2 =R A D2 / 1 0 « 



ISN 

0216 


IF (RAC 2 .GT .ERRrtl 2) L ARGE3 =1 



ISN 

0218 


IF (R AD 3. GT .ERRHI 3 ) RAD3 =R AD3/1 0. 



ISN 

022 C 


IF ( IT IME5. LT. I IDAT) GO TO 40 


0 149 

ISN 

0222 


IF ( I IC AT .EQ . I T I ME 5 .A ND • 1 TI ME6. LT. I H ) GO TO 40 


0150 

ISN 

0224 


IF ( ITIMES.EU. I IDAT 1 ) GO TO 1C04 


0151 

ISN 

0226 


IF (IGCP1.NE.I) GO TO 2 0C8 


01 52 

ISN 

0228 


I T J=ITIME6 


0153 

I SN 

0229 


I TK=( I T I ME 6/ 1 00) *1 00 


0 154 

ISN 

0230 


IKL = I T J— I T K 


0 155 

I SN 

0231 


IF ( IKL . E Q • 0 ) GO TO 81 0 


0156 

ISN 

0233 


IF( IKL.LT.I CIFF) GO TO 310 



I SN 

0235 

C 2 

IF ( ISTART.EQ. 1 ) GO TO 75 


0 157 

I SN 

C237 


IF ( I START-. EC « 2 ) GO TO 70 


04 58 

ISN 

0239 


IF { ISTART .c 0 . 5 ) GC TC 80 


0159 

ISN 

0241 


IF ( ISTART.EQ .4 ) GO TG 90 


0160 



C CHECK TO SEE fchICH GRAPH IS TC BE PLOTTED 


016 1 

ISN 

C 24 3 

8 1C 

IF (RANI .GT. 1 30. C) GO TO fill 


0 162 

ISN 

0245 


JK= JK* 1 


0163 

ISN 

0246 


MHCURS < JK)=ITJ/1 CC 


0 164 

ISN 

0247 


RANCE5( JK)=RAN1 


0165 

ISN 

0248 


ERROR(JK) =R ADI 


0166 

ISN 

0249 


ERRCR1 ( JK)=RAD2 


0167 

ISN 

025C 


E RR CR 2 ( JK > = R A D 3 


0168 

ISN 

0251 


I F ( NSS5 .EQ. C ) PRINT 336 . JK .MHUU RS ( J K ) 



I Sry 

Q25 3 

336 

FCRMAT ( 1H . *MHGU«S ( *♦ 13-.* ) = •♦1-6) 



ISN 

0254 


CC TC €2 


0 172 

ISN 

0255 

2CCQ 

IF ( IT IME5.LT.IT IM88) GC TO 40 


0 173 

ISN 

0257 


IF ( IT IME5.GT. IT IM68) GO TO 2006 


0174 

ISN 

0 259 


IF < IT IME6.LT. IT IM99 ) GO TO 40 


0175 

ISN 

0 26 1 


I GOP 1 = 1 


C 176 
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D-5 



ISN 0262 
ISN 0264 
1-8 N G 265 
ISN 0266 
ISN 0267 
ISN 0269 
ISN 0269 
ISN 0270 
I-8N G27 1 

ISN 0272 
ISN 0273 

ISN 0274 
ISN G27-S 
ISN 0277 
ISN 0278 
ISN 0260 
ISN C 2 £ 1 
ISN 0263 
ISN S 2-64 
ISN 0265 
ISN 0286 
ISN 0287 
ISN 0283 
ISN 0289 
ISN 0290 
ISN 0282 
ISN 0293 
ISN 0294 
ISN 0295 
ISN C286 
l-SN 0297 

ISN 0288 
ISN 0300 
ISN 0 30 I 
ISN 0302 
l SN 0303 
ISN C3C4 
ISN 0305 
ISN 0 306 
ISN 0307 
ISN 0303 
ISN 0310 
ISN 0211 
ISN 0212 
ISN 0313 
ISN 0314 
ISN 0316 
ISM 0318 
ISN 0320 
ISN 0322 
ISN 0324 
ISN 0326 
ISN 0328 
ISN 0330 
ISN 0332 
l-SN 0 3 3 4 
ISN 0336 
ISN 0337 
ISN 0338 
ISN 0339 
ISN C 34 0 
ISN 0-34-1 
ISN 0342 
ISN C- 34 3 
ISN 0344 
ISN 0345 
ISN 0346 
I-SM 034? 
ISN 0349 
ISN 0 350 
ISN 035 1 
I SN 0352 
ISN C 353 


IF(NSS5.EQ.Q ) PRINT 335*1 TIME5# I T I ME6 
33S FORMMUH ,'ITIME5= *16 . 2X * • I T I ME 6 = '.I6> 

GO TO 62 
2 C C 6 PRINT 2CC7 

2CC7 FCfiMATC 1H1 .27HTI ME SPAN INCORRECT CN THIS VC REPORT) 

PRINT 2C08 . I T I N66 . I TI M99 

2CC9 FORMAT ( IHO * 3 4HEND TIME CF LAST CORRECT VC R EPORT * 1 X * 1 6 * l X , I 4 ) 
PRINT 20 10 , I T1 . I T2 * I T3 . I T4 

2-0-1 0 F-CRM A T ( 1 Hut 30H7 I Mfc SPAN CN THIS VC REPORT I S * l X , 16 • IX * I 4* 1 JX • I 6* 

1 1 X . 14 ) 

FR I N T 2C l 1 

2011 FORM AT ( 1H0 . 6 8H0CPL T4 WILL PROCEEO TC NEXT VC REPCRT TO SEARCH FOR 
1CCRRECT TIME SPAN) 

GO TO 127 

811 IF (UK.EO.Q) GO TO 61 
I T TK = I TK/1 DC 

IF ( I T TK «L T *M HOURS! JK ) ) I T TK = ITT K+24 
I HD IF = I T TK- NHOOR S { JK ) 

IF (IHCIF.GT.4) GC TO 61 
GO TC ez 

61 JK=JK4 1 

M FOURS ( J K ) = I T J/1 00 
R ANG E5( JK ) sRAN 1 
ERRCRt JK)=RAC1 
ERR CR 1 ( JK > =fi AC2 
ERR CR2( JK ) =RAD3 

IF ( NSSe (£6« 0 ) PRINT 336 . JK .M HO URS ( J K ) 

CU TO 62 

75 XP AN 1 =R AN 1 
TCP 1 =0 

T CP 2 =0 
TCP 3=0 

CALL AEREAD (ITIME5. ITIME6* KAOl , IEXP 1 , RAD2. IEXP2. RAD3, 
IIEXP3* RANI. IEXF4, I 3EQF ) 

[F ( I3E0F »E0 • 1 ) GC TO 10C4 
TIMOCT = 0. 

lOLTIM(l) = ITIME5/10000 

I0LTINC2) = ( I T I ME 5 - 1 0000* IOUT IM < 1 )>/ 100 

I OUT l M ( 2 ) = ITIME5 - < l 00C 0 * IOUT I M ( 1 ) 4 10 0 * I CUT I v< 2 ) ) 

IOUT IM (4 ) = ITIME6/1 00 

lOLTIM(E) = 1TIMEC - l 0 C* I OU T I M ( 4 ) 

CALL TCONVO (TIMDUT. I OUT i M , SEC) 

T I MOL T = TIMCLT ♦ T3DIFF 
IF ( T IMDUT.GE. TIMENO ) GO TO 1004 
RAC1=RAD 1* 10 .0** IEXP 1 
RAC2=RAC2*1J.0**IEXP2 
RAL)J = RAC 3*10.0** IEXP 3 
R AN1 = < RAN1*1 C .C** IEXP4)/10G0 .0 
IF(RAD1 •LT.ERRLCI ) RAD1=ERRL01 
IF(RAD£#LT.ERRL02) RA02=ERRL02 
IFIRAD2.LT. EHRLU2) RA03=ERRLQJ 
I F (RAO I .GT .ERRHI I ) TOPl=l 
IF(RAD1 *GT .ERRHI 1 ) R AD 1 =RAD 1 / 1 0 • 

IF<RAC2.GT.ERRHI 2 ) T0P2=1 

IF( R AC2 .GT.ERRHI2 ) R AD2=RAD2/ 1 0 . 

IF (RAD2 .GT. ERRHI 3) T0P3=i 
IFIRA02 .GT.ERRHI 2 ) RAD3=RAD3/l 0. 

IF ( XRAN 1 .GT .RANI ) GO TC 76 
I F ( XRnAN 1 .EG .RANI ) GO TO 76 
X 1 ( 1 ) = F A N 1 
V 1 ( 1 )=RAD1 
X 2 ( 1 >=RAN1 
V 2 ( 1 )=RAD 2 
X6( 1 )=RAN 1 
Y6U)=WAOJ 
I = 1 

I T ME 1 = 1 TIMES 
I TME 2= I T I ML 6 
I SET 1 = 0 
I ST A R T = 2 

I F (-N SS5.EQ.0) PRINT 211.ISTART 
211 FORM A T ( 1 H . « ISTART=* . I 3) 

PRINT 209 

2C 5 FCPMATUH .'PERIGEE TO APOGEE PASS') 

GC TC 40 

76 IF (IStTl.Nc.l) GC TO 77 


0 180 

0 1 ei 

0 183 
0165 


0 187 


0193 
0184 
C 1 8 5 
C 196 
0 197 
0 198 
0199 
C 2 00 
020 1 
C 2 02 

0203 

0204 
0 205 

0209 
02 1 0 


02 16 
02 17 
02 l 8 
02 19 


0226 

0227 

0228 

0229 

0230 

0232 

0233 
C 2 34 
0237 
02 38 
0 2 39 
0240 


0249 

025 l 
0252 
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D-6 



I-SN 

0 3-55 


PRINT 210 

0253 

ISN 

0 3 S6 


X2{ i ) =RAN l 

0254 

ISN 

0 357 


V 31 1 1 = RADl 

0255 

ISN 

0358 


X 4 ( 1 ) =R A N l 

0256 

I SN 

C359 


Y 4 ( 1 ) sRA02 

0257 

I SN 

0 3 6C 


XS ( 1 ) = R A N1 

0 258 

ISN 

C361 


Y5 1 1 J=RAD J 

0259 

ISN 

0362 


ITME 1= I T I ME 5 

0260 

ISN 

0363 


I TME2= I T IMF 6 

026 l 

ISN 

0 364 


I =1 

0263 

I SN 

0365 


ISTART = 4 

0265 

I SN 

0366 


I SE T 1 = 2 

0266 

4S-N 

C36? 


IF.(NSS£*EQ* 0) PRINT 211 •ISTART 


ISN 

C369 


GO TO 40 

0270 

ISN 

0370 

77 

X 7 ( 1 ) =RAN 1 

C27 1 

(SN 

037 1 


PRINT 210 

0276 

ISN 

0372 

21 C 

FCHMATUH • • APOGEE TO PERIGEE PASS*) 


I SN 

0373 

403 

Y 7 ( 1 ) = RAO 1 

0278 

l-SN 

037A 


X6< 1 ) =R A M 

0 279 

ISN 

0375 


Y6( 1 )-RAD2 

0280 

ISN 

0376 


X 9 ( 1 ) sRAN 1 

0281 

ISN 

0377 


Y9( 1 )=RA03 

C2e2 

ISN 

0375 


I TMEl =1 TIMES 

02 63 

ISN 

0379 


I TME 2 = I T I ME 6 

0264 

i-SN 

0380 


1 = 1 

0286 

ISN 

0 361 


I 2= 1 

0237 

ISN 

03£2 


I STAR T = 2 

0288 

ISN 

0383 


I SET 1 = 5 

0289 

(SN 

C 384 


I F ( N S S 6 • EG • 0 ) PRINT 211. ISTART 


ISN 

C 386 


GC TO 4 C 

0293 

I-SN 

0387 

70 

W? 4«ANl-*G.T*XiiAN44 GO JO 500 

0294 

ISN 

0389 


XR AN 1=R A N t 

0295 

ISN 

0390 


IF < ISET .EQ . 1 ) GC TO 71 

0296 

ISN 

0 392 


RCI F= ( X 7 ( I 1-RANI ) 

0297 

ISN 

0 39 3 


If ( R C I F • GE *RFRE 0 ) GO TO 71 


ISN 

0395 


GC TC 40 

0299 

ISN 

0396 

7 1 

1 = 14 1 

0300 

ISN 

0397 


X7( I )=RAN 1 

0301 

ISN 

0 398 


Y7( l 1 — R A 0 l 

0302 

ISN 

0399 


*6(1) =R AN 1 

0303 

ISN 

0400 


Yf( I ) =RA02 

0804 

I SN 

0401 


X 9 ( I ) =R A Nl 

• 0305 

ISN 

0402 


Y 9 ( I ) = RAD 3 

0306 

ISN 

C 4 C 3 


I SET =C 

0307 

ISN 

0404 


IF( LARGE l *50 . 1 1 TOPl=l 


ISN 

0406 


I F( L ARGE2. EQ. 1 ) T0P2=t 


ISN 

0408 


IF (LARGE2 .E C . 1 ) TOP 1 = 1 


ISN 

04 10 


GC TC 40 

0309 

ISN 

04 1 1 

60 

IF (RAN1.LT.XRAM) GC TC 600 

0310 

ISN 

04 13 


XR AN 1 = R AN 1 

C2 1 1 

ISN 

0444 


IF 4-lS-E T*£Q •-! ) GC TO 81 

0212 

ISN 

04 16 


RC I F = (RANI — XI (I)) 

0213 

ISN 

04 17 


IF (RDIF .GE .KFREC) GO TO ei 


ISN 

04 19 


GC TC 40 

03 15 

ISN 

0420 

61 

1 = 1-41 

03 16 

ISN 

042 1 


X 1 ( I ) =R AN 1 

02 1 7 

ISN 

04-52 


Y 1 ( I ) =R A 0 l 

C 2 1 8 

ISN 

0 4 23 


X 2 ( I ) = R A N 1 

03 19 

ISN 

0 4 24 


Y 2 ( I ) =R A D2 

C3 2C 

I SN 

0425 


X6 ( I ) =R A N 1 

C 2 2 1 

ISN 

0426 


Y 6 ( I ) = RAD 3 

0 J22 

ISN 

0427 


I SET = 0 

0323 

ISN 

0428 


1 F( LARGE 1 .tQ • 11 T CP 1=1 


ISN 

04 30 


l F( LARGE2 .E J . 1 > TOP2 = l 


ISN 

0* 32 


IF( LARGER . EO. 1 ) TGP3=1 


ISN 

0434 


GC TC 40 

0224 

I SN 

C 4 35 

90 

IF (RAM.GT.XRAN1) GO TO 700 

0225 

I SN 

0437 


XfiAM =RA Nl 

0 2 26 

I-SN 

C4 3S 


IF (ISLT.EU.l) GO TO 91 

0327 

ISN 

C 4 4 0 


R C IF = ( >3 ( I ) - R A N l ) 

0328 

ISN 

0 44 1 


IF (RDIF .GE .RFREQ) GO TO 91 


ISN 

044 3 


GC TC AC 

0330 

ISN 

C444 

9 l 

1 = 14 1 

0221 

I SN 

C 4 4 5 


X 3 ( I ) = R A M 

0332 

ISN 

C446 


Y 3 ( I > = RAD I 

0333 

ISN 

C447 


X4 ( I ) = RAN 1 

0334 
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ISN 0A48 
ISN 0449 
ISN 0450 
ISN C451 
I SN 04 52 
ISN 0454 
ISN 0456 
ISN 0458 
ISN 0 4 59 
ISN 0460 
I-SN 0461 

ISN 0462 
ISN C 463 
ISN 0464 
ISN 0466 
ISN C467 
ISN 0468 
ISN 04 7-0 
ISN 0471 
ISN 0472 
ISN 0473 
ISN 0474 
ISN 0475 
ISN 9477 
ISN 0479 
ISN C 480 
ISN C 4 8 1 
ISN 0482 
ISN 04e3 
IS*. 0 484 

ISN 0485 
ISN 0487 
ISN 0488 
ISN C489 
ISN 0490 
ISN 049 1 
ISN 0493 
ISN 0495 
ISN 0497 
ISN 0499 
ISN C501 
ISN 0 SC 3 
ISN C 50 5 
ISN 0507 
ISN C 5 09 
I SN 0510 
ISN 0511 
ISN 0512 
ISN 0513 
ISN 0514 
ISN 0515 
ISN 0516 
ISN 0517 
ISN 0513 
ISN 0520 
ISN 0522 
ISN 0524 
ISN 0525 
ISN 0526 
ISN 0527 
ISN 0528 
ISN 0529 
ISN OSJO 
ISN 0531 
I SN 0532 
ISN 0533 
ISN 0534 
ISN 0535 

ISN 0536 
ISN 0538 
ISN 0539 
I-SN 0540 
ISN 0541 


Y4( I ) =RA02 
X £ ( I > ANl 
Y 5 ( I » = R AD3 
I SE T =0 

IF < LARGE WE 0. I ) TGPt=l 
I F ( L AR GE 2 .E J • 1 ) TCP2=l 
I F( LAPCEJ .EQ . 1 > T CP 3= 1 
GO TO 40 

1CC4 I EN 0 = 1 T4 /ICO 
PRINT 1 C 1 4 

10 14 FCRUAT (IH* 7-HAT ICG4) 

C IF IRST= 1 

I 1= I 
I =1 B 

IF(N$S6. cQ • 1 ) GO TO 2004 
L STAR T = C 


ITGPH= 1 

IF C T IME V .EQ .T INLY ) GC TO 31 
CALL BSFTAP ( AS 1 
GC TC 31 
749 PRIM 4 C 9 9 
4099 FORMAT ( 6H ITGPH) 

9689 READ (A £ * 39 J ) CHC 

IF ( CFC .NE • STR T 1 ) GO TO 9889 
I F (IT AP b tN£ • 1 ) GO TO 752 


K J= C 

7 £ C CCNTINLE 
UPPER 1 = 0 
UPPER2=0 


UPPER 2 = 0 

CALL AS R CAD (ITIML5. ITIME6* WADI. 
1IEXP3. RANI. IEXP4, 1 3E OF ) 

IF ( 1 3LCF . E C • 1 ) GO TO 20C4 
RADI=RAU1*1C.C**IEXP1 
RAD 2 = R AD 2* l C . C**IEXP2 


RAD 3 =R AD 3* 
RAN l = ( P AN l 
I F( RAC 1 .LT 
I F ( RAD2.LT 
IF ( fiADj.LT 
IF (RAD 1 • GT 
IF (RAD 1 .gT 
I F (RAD 2 .GT 
l F ( R AC 2 .GT 
IF( RAC3. GT 
I F ( RAD2.GT 


13. J**IEXP3 

• 10. 0* * I E X P4 1/10GC.0 
•ERRLOl) RAU1=ERRLC1 

• c R RL02 ) RAD2 = EWRLU2 
•£RRL03> RAC3=£RRL03 
. £ RRH I 1 ) UPPER 1=1 
•ERRHI1) RADI =RAD1 /IC . 

• ERRH12) UPPER2 =1 
•ERRHI2) RAD2=RAD2/1 C. 
. E RR h I 5 ) UPPER 3= 1 
•ERRHlJ) R AD3=RAD3/ 1 j • 


I EXP 1 . R AD2 • 


IEXP2. RAD 3. 


1 = l 

T IMOUT = :. 

10UTIM(l) = I T I ME 5/1 0 GO 0 

IOL T IM ( 2 1 = ( I T I ME £ - 1 0 00 C * I CUT I M ( 1 )) / 1 00 

1CUTIM2I = I TIMES - ( 1 0000* ICUT 1 M ( 1 ) 1 00 ♦ 1 OU T I M ( 2 ) ) 

I CL T 1M ( 4 I = ITIME6/100 

ICL T IM ( £ ) = I T I ME 6 - 1 00*1 OL T I M ( 4 ) 

CALL TCCNVO (TIMDUT, ICUTIM, SEC) 

1IMCLT = TIMOLT ♦ T 3D IFF 
IF < 1 IMCLT.GE .T IMEND ) GO TU 2004 
IF ( IT IMES.LT. I ICAT ) GO TC 750 
IF ( ITIMcE.EC.I ICATI ) GO TO 2005 
X T IME = FLCAT ( I T IME6 ) 

XT IME=M IME / 100. C 
7 £ I XFRS(1)=XTIME 
ICA 1E = IT I ME 5 
RANGE ( 1 > =M AC 1 
RANGE 1 ( 1 ) =RAC2 
R ANGE2 ( 1 ) =R AD 3 
7 £ 2 CCNTINLE 
8 TG1=0 
D I G 2 = 0 


0 IG3=0 

CALL A5READ (ITIME5. ITlMb6. KADI. IEXPI. HAC2. IEXP2. RAD3. 
I IEXP2. RANI. IEXP4. I 3 EOF ) 

IF ( 1 2 ECF .EC . 1 ) GC TO 2CC4 
RAC1=RAC1*10 .O** IE XP 1 
RAD2=RAC2* 1 C .0** IEXP2 
R AD 3 -U AO 3 » 1 C-.Q** IE XP 3 
RAN 1=( RANI* 1 0. 0**1 EXP4 ) /I OuO.O 


0335 

0336 

0337 
0 3 38 


0339 
C 34 0 


034 I 
0 34 2 
0343 

0 246 

0347 

0348 

0250 


0 355 

0356 

0357 

0358 


0267 

0268 
0269 
0370 
C271 
C 3 73 
C 3 74 
0375 
0 276 


0379 

0380 
0281 
0282 
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D-8 



ISN 

0542 


IF< RAD 1 .LT .ERRLC1 ) RAD1=ERRLC1 




ISN 

0544 


IF(RAD2.LT.=RRLC2 ) RAD2 =E RRLQ2 




ISN 

0546 


IF<RAC2.LT.ERRL03) RAD 3=£ RRL03 




ISN 

05 4d 


I F( RAC1 • GT.ERRHI 1 1 rfIGl=l 




ISN 

0S5Q 


IF < R AD 1 • G T • g-RRH I 1 4- R AO 1 =R AD 1 ✓ 1 0 . 




I SN 

0552 


IFIRAD2 .GT.ERRHI 2) 3IG2=1 




ISN 

0554 


IF (RAD c .GT.ERRHI 2) R AD2 =R AD2 / 1 0 • 




ISN 

0556 


I F (RAC 2 .GT .ERRHI 2 ) 8IG3 = 1 




ISN 

0556 


I F ( R AC 3 « GT .tRRHl 3 ) RAD3— RAD3/1 0. 




ISN 

0560 


T IMCLT=C • 




ISN 

456- 1 


ICCTIM<44 -= -4 T IME5/1 0000 




I SN 

0562 


I CL T I M ( 2 ) = ( I TIM E5 - 1 0000 * I0UT I V ( l M / 1 3 0 




ISN 

0563 


iciiiim;) = iriM£5 - < ioooo*ioutim i ) ♦ 100 * tout im< 2 ) 1 



I SN 

0564 


I CL T IM ( 4 > = I T I ME 6/100 




ISN 

0565 


ICLT IM ( £ ) s ITIME6 - 1 03*1 OUT I M(4 ) 




ISN 

0566 


CALL T CUN V0 ( T I MDUT • I CUT I M . SEC) 




-ISN 

<1567 


TIMOLT = TIMOLT * T 3D IFF 




ISN 

0568 


IF < T l MDL T. Gt • TI MEND. AND. ITAPE.EO.O ) GC TC 

20C4 



ISN 

0570 


IF ( IT IME6.EC. 0) GC TO 781 



03e7 

I SN 

0572 


ITJ=( IT IME6/1 30) *IC0 



0368 

ISN 

0573 


IKL - IT IME6-IT J 



0389 

ISN 

0574 


IF (IKL.EU.O) GO TO 7o2 



C 390 

ISN 

0576 

7 6 1 

X CH6CK = FLU A T ( ITIME6) 



C391 

ISN 

0577 


XCHE CK =X CHECK / 1 00 • 0 



0392 

ISN 

057b 


IF ( ITIME8.EQ.10ATE) GO TO 755 



0393 

I SN 

0580 


c\ 

n 

n 



0394 

ISN 

058 I 

755 

F IKL=FLCA T( I T J/1CC) 




ISN 

05e2 


XCF ECK = FLUA T < I KL 1 / 69* 




LSN 

0563 


XC FECK = XCHEC X*HI KL 




ISN 

0564 


Y CF ECK = XCHEC K— XHRS ( I ) 




ISN 

0565 


IF (YCHECK.LT. TFREO) GO TO 752 




ISN 

0587 


1 = 1*1 



0397 

I SN 

0588 


XFRS ( I > = XCHECK 



0398 

I SN 

C 5 69 


RANGE ( I ) =RAD 1 



0 399 

I SN 

0590 


RANGE 1 ( I ) =R AD 2 



0400 

ISN 

0591 


R ANCE2 < I )=R AD 3 



C 4 0 1 

ISN 

C 59 2 


IFieiGl.EQ.il UPPER 1 = 1 




I SN 

0594 


IF (8 IG2.EQ. 1 1 UPPER2 =1 




I SN 

0596 


IF (8 1G3.EQ. 1 ) UPPER3 = 1 




ISN 

C 59 0 


GC TO 7*2 



C 402 

ISN 

0 5S9 

762 

If(KJ • EQ • 0 ) GO TO 7 « 5 




ISN 

060 I 


IFC I 1 IME6/ 100.LE . IHOUR2 (K J) ) GO TC 781 




ISN 

0603 

7 € 5 

K J=K J* 1 




I SN 

0604 


I FOUR 2 ( K J 1 = I TIME 6/1 JO 



0404 

ISN 

06C5 


RANGE7 ( K J 1 =RAN1 



0405 

I SN 

0 6C6 


AES IC < K J )=RAD 1 



0406 

ISN 

060-7 


A6S IC 1 <K J 1=«AD2 



0407 

ISN 

0608 


AESIC2IK J )=RAD3 

* 


C 4 C 8 

ISN 

0609 


GC TC 781 



C4C9 

ISN 

06 10 

7*3 

IF { I .EG . 1 ) GO TO 750 



04 10 

I SN 

06 12 


IF (XChECK.EC.Q. 01 GO TO 808 



04 1 1 

I SN 

0614 


1 = 1-1 



C4 12 

ISN 

06 15 


I S T = C 



C413 

ISN 

06 16 


I ST 1 =0 



04 14 

ISN 

06 17 


*SK-l F-=-l 



0415 

I SN 

06 18 


MTYPE= 1 




ISN 

06 19 


CALL TITLES < MT YPE . MSK I P ) 



C4 19 

ISN 

0620 


X IX = XNCRMZ( AFCDE . C. ) 




ISN 

0621 


IF ( XFiRSI 1 1 .NE .0 .0 .AND. IF IRST.EQ. 1 ) IST = XHRS(l> 


C421 

ISN 

0623 


PRINT 337, l CA TC. I T IMEO 




ISN 

06 24 

33 7 

F-GWM AT ( 1 H « ■ IC A T-E- = - * . I 6 • 2 X , • I T I ME 6 -=• ' • I 6 1 




I SN 

0625 


CALL LEGNDG < AMODE ,642. . 121 • .12 , I 2HHCLRS CF 

CAY ) 



I SN 

0626 


IF< UPPER 1 .EQ. 1 1 







1 C ALL LtC ND G ( A MODE ,237. ,612. ,94, 94 H * * * * CNE 

OR M CPE POINTS 

FAVE EXC 





2EE0EC LPPER Y LIMT ANU ARL DIVIDED 6Y 10 BEFORE PLOTTING 

** ** ) 


ISN 

0623 


CALL S E T SM G ( A MODE ,14,0. ) 




ISN 

C629 


C4LL L INESGI AMOOE, 1*1. XHRS .RANGE) 




I SN 

0630 


CALL SE T SMG < A MODE ,14,3. 1 




ISN 

0631 


J NO TE = C 



0430 

I SN 

06 32 


CALL ALTCK ( C v ,JNDTE) 



0431 

ISN 

0633 


I ND T E = 4 



0432 

ISN 

06 34 


CALL CATAPT ( ICATE, I ST , I T I VE 5 . I ST 1 , X I X , INCTE) 



I-SN 

06-35 


IF ( NCCMR.EQ. 1 1 GO TO 960 




ISN 

0637 


F TYPE =2 




ISN 

C 6 38 


CALL TITLES ( M T YFE , MSK IP) 



04 35 

ISN 

0639 


CAt_L LEGNDG (AMODE .642. .121. .12. 12HHCLRS CF 

DAY ) 
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ISN C640 IF < UPPER2 .EG . 1 ) 


1C/»LL LECNDG( A MODE ,2 37. , 8 12. ,94 , 9AH* * <M> CNE OP MOPE POINTS HAVE EXC 
2EEDEC UPPER V LIMIT ANu AkE DIVIDED eY 10 BEFORE PLOTTING **++ ) 


ISN 

0642 


CAlL SET SMG ( AMQD6 , 1 4 ,0 • > 



ISN 

064 3 


CALL L I N t£ SG ( A M OD E » I + 1 .XHRS.RANGtl ) 



ISN 

0644 


CALL SETSMG ( AVODE ,1 4,3. ) 



ISN 

Of 46 


JNDTE= 1 


0443 

ISN 

0646 


CALL ALTCK ( K J , JNDTE ) 


C444 

ISN 

0647 


INC IE =3 


C 4 45 

ISN 

0643 


CALL CAT AP T( IDATE.IST ,ITIME5,IST1 , X 1 X , [NOTE) 



ISN 

0649 


V TYPE = 3 



I SN 

C65C 


CALL TITLES < m TYPE , M SK I P ) 


0448 

ISN 

0651 


CALL L EG NDG { A MODE ,642. • 1 2 l • • 12.12HHCURS CF CAY) 



ISN 

0652 


I F ( LPPER 2 .LG . 1 ) 






10 ALL LE6NDG ( AMCDE ,237. , e 12. , 94* S4H* *** CNE OR MORE POINTS 

HAVE EXC 





2EFDEC UPPER Y LIMIT AND ARE DIVIDED EY 10 dEFGRE PLOTTING 

W* ¥* ) 


ISN 

0654 


CALL StTSMGI AMODE.14,0. ) 



ISN 

0655 


CALL L I NESG ( A MODE « I + 1 , XHRS ,RANGE2) 



I SN 

0656 


CALL SFTSMG( AMODE.14,3. ) 



ISN 

0657 


JNC TE = 2 


C 459 

ISN 

0653 


CALL ALTCK (K J , JNOTE ) 


C460 

ISN 

0 6 59 


I NC T E =3 


C 4 6 l 

ISN 

0660 


CALL CATAPTI 1 0 ATE , I ST . I TX ME5 • I ST 1 , X I X , I NOTE ) 



ISN 

066 l 

SCO 

CONTINUE 



ISN 

066-2 


K J — 0 


0463 

I SN 

0663 


IF I R £ T = 0 


C456 

ISN 

0664 


M SK IF ~C 


C457 

ISN 

0665 


I TAPE=1 


C458 

ISN 

0666 


IF (NOGC.EQ.l) GO TO 2004 


C 464 

ISN 

0669 


IF (ITIME5.EC.IT2) GO TO 757 


0 465 

ISN 

0670 


GC TC 750 


0466 

ISN 

0671 

757 

IF ( IEND .EG . I T I ME C ) GO TO 2004 


0467 

ISN 

0673 


IF I I T IME6 .EQ.IT4) GO TO 2004 


0466 

ISN 

0675 


PRINT l C 7, I TI ME6 , IEND , IT4 



ISN 

0676 


I TAPE— G 


0469 

ISN 

0677 


ISAVE=1 


0470 

ISN 

0673 


KK= 1 


0471 

ISN 

0679 


GO TC 750 


C472 

ISN 

0680 

9£3 

I TGPH=1 


C482 

ISN 

068 1 


IF ( XO-ECK.EQ.O.O > GO TO 80 8 


0483 

ESN 

06 ©3 


I SAVE-C 


0464 

I SN 

0684 


I TAPE= 1 


0465 

ISN 

0685 


L START = l 


04 86 

ISN 

0686 


GO TO 753 


0467 

ISN 

0687 

8 C 8 

XCHECK=24 .0 


C 488 

ISN 

0 6 83 


1 = 1 + 1 


C489 

I-S-N 

06-80 


XF-RS-I I )=XCHECK 


04 SO 

ISN 

0690 


RANGE! I ) =R AC l 


C 4 9 1 

ISN 

0691 


RANGE 1(1) =R AD2 


C 492 

I SN 

0692 


R ANGE 2 ( I ) = RAD 3 


C 49 J 

ISN 

C69J 


IF (b IG l .EQ. 1 ) UPPER1 =1 



ISN 

0695 


I F( g IG2 .EO • I ) UP PE R 2 = 1 



I-SN 

0697 


I F<8 1G3.EG . 1 ) OP PE R 3= 1 



ISN 

0699 


IF (ISAVL.EC.l) GO TO 953 


C 4 9 4 

I SN 

070 I 


GC TC 753 


0 49 5 

I SN 

C7C2 

2 C C 4 

I T I M 6 6 = I T 2 


0496 

1 SN 

0703 


IT IM99 = I T 4 


C497 

ISN 

C 7 04 


IF ( IT IMtS.EQ. I I0AT1 ) GO TO 2C95 


0498 

ISN 

0 70 6 


I GOP 1 =C 


C499 

ISN 

0707 


IF ( 7 IME Y. EG. T IMLY ) GO TC 1007 



ISN 

0709 


NC GC =0 


C 502 

ISN 

C 71 0 


IF ( T IME Y ,FC • TI ML Y. AND. I TAPE .EQ. 1 ) GC T3 1037 


0503 

ISN 

0712 


ie= i 


C5C4 

ISN 

07 13 


i = i i 


C5C5 

ESN 

07 14 


I TGPH = C 


CSC 6 

ISN 

07 16 


IF{ I «L T . 2 ) GO TO 40C7 


0507 

ISN 

0717 


IF ( ISTAKT.EQ.2) GO TO 1005 


0508 

ISN 

0719 


IF ( I3TART.E0.3) GO TO 10C6 


0509 

ISN 

072 1 


M SK IP = C 


C5 1 1 

ISN 

0722 


M T YP E = 1 



ISN 

0723 


CALL TITLES(MTYPE.MSKIP) 


CSi 3 

ISN 

07 2** 


CALL SET SMG( AMOOE . 1 4 ,0 . ) 



ISN 

0725 


CALL L INESG! AMOOE , I . X3 , Y3 ) 



ISN 

C 726 


CALL SETSMG(AMCDE.14,3. ) 



ISN 

0727 


DC 1 19 J= 1 , I 


C 5 20 

ISN 

0728 


SAVE 1 ( J )=Y3( J) 


C 5 2 1 
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D-10 



ISN 

G 72 9 


SAVE2I J )=Y4< J) 

C 522 

ISN 

0730 


S A V E 3 ( J )=Y5 I J ) 

0523 

ISN 

0731 


S A V E 4 ( J > = X3 ( J > 

C £ 24 

ISN 

0732 

1 19 

CONTINUE 


I SN 

0733 


X 3< 1 ) = X 3 ( I ) 

0526 

I SN 

0734 


Y3( 1 > = Y3( I ) 

0527 

F&N 

0735 


X 4 I 1 ) — X4 ( I > 

0 5 28 

ISN 

C736 


X5( l >=X5 ( I ) 

0529 

ISN 

0 737 


Y 4 ( 1 >= Y4( I > 

0530 

ISN 

0730 


Y £ { l ) = Y5 ( I ) 

053 1 

ISN 

0739 


12=1-1 

0532 

I SN 

0740 


I GCP= 1 

05 33 

ISN 

0741 


GO TO 1007 

C 534 

ISN 

0742 

1CC5 

CONTINUE 


ISN 

07-4-3 


MSK IP^=-C 

C E36 

ISN 

0744 


*TYPE= 1 


I SN 

C 7 45 


CALL TITLES I MTY PE • MSK l P) 

C 5 38 

I SN 

0746 


CALL SETSMGI AMODE • 1 4 , 0. ) 


ISN 

G747 


CALL LlNESGIANCDc.I *X7.Y7I 


ISN 

0746 


CALL SETSMG IAMODE • 1 4.3. ) 


ISN 

0-749 


CO 1-22 3=1*1 

0545 

ISN 

C 750 


SAVEl C J) =Y7( J) 

C £ 46 

I SN 

0751 


SAVE 2IJ)=Y3(J) 

0547 

I SN 

0752 


S A vE 3 ( J ) =Y 9 ( J ) 

0548 

I SN 

C 753 


SAVE 4 ( J)=X7( J) 

0549 

ISN 

0754 

122 

CONT INUE 


ISN 

07-55 


X7( l ) = >7( I 1 

C 5 5 1 

ISN 

0756 


Y7< 1 J = Y 7 ( [ ) 

0552 

ISN 

0757 


xe ( i ) =x£ ( i ) 

0553 

1 SN 

0758 


X9 < 1 ) =X 9( 1 ) 

0554 

I SN 

0759 


Y8< 1 ) = Y8(I) 

0555 

ISN 

0 76 3 


>91 1 )=Y9 (I ) 

C 5 56 

ISN 

0 76-1 


12=1-1 

0557 

ISN 

0762 


I GOP = 1 

C £ 5 9 

ISN 

0763 


GO TC 1 C 07 

0559 

I SN 

0764 

10 C 6 

CCNT INGE 


I SN 

C765 


*SK I F=2 

056 1 

ISN 

C766 


M T YPE = 1 


ISN 

0767 


CALL TITLES IMTYPe . MSKI PJ 

C 563 

ISN 

0 76e 


CALL SETSMGI AMODE . 14 .0 . ) 


ISN 

0769 


C Al_L L INESGI AMODE • I .XI ,Y1 ) 


ISN 

C 770 


CALL SE 1 SMGI AMODE * I 4 .3 . 1 


I SN 

0771 


OC 12 E J = l ♦ l 

0570 

I SN 

0772 


SAVE 11 J) = Y 1 ( J) 

C 5 7 I 

ISN 

0773 


SAVE 2 ( J J =Y2( J) 

0572 

ISN 

0 774 


SAVE’ I J ) =Y6 ( J ) 

0573 

ISN 

C 7 7 5 


S AV64I J ) = X 1 ( J ) 

0574 

ISN 

0776 

125 

CCNT INLE 


I SN 

0777 


XI ( 1 ) =X 1 ( I ) 

0576 

I SN 

C 778 


Y 1 ( 1 )=Y 1 ( I ) 

0577 

I SN 

0779 


X 2 I 1 ) = X2(1 ) 

0 5 78 

ISN 

0780 


X 6 ( 1 ) = X 6 I I ) 

0579 

ISN 

0761 


Y2< 1 1 = Y2( II 

0580 

ISN 

0782 


YC( 1 ) = Y 6 ( I ) 

0591 

ISN 

0783 


12=1-1 

C5e2 

I SN 

0784 


I GCP = 1 

058 3 

ISN 

0785 

1007 

JINX=1 

0584 

ISN 

C 7 86 


IF (TIME Y.EO.TI MLY) GO TO 4008 


ISN 

0 788 

4C ( 7 

I = C 

cses 

ISN 

0789 


IF( IGOP.EQ. I ) 1=1 


I SN 

0791 

4CC8 

I SET =1 


I SN 

C792 


I BL A P=C 

0587 

I-SN 

0-7-93 


CALL TAPE8418LAP4 

0588 

I SN 

0794 

1 27 

CGNT INUE 


ISN 

0795 


ix i v = e 


ISN 

0796 


GC TO 25 

C £ 9 5 

ISN 

0797 

SCO 

I 1ME2 = I TI ME 5 

0596 

ISN 

C 79 Q 


ITM£4=ITI ME £ 

0597 

ISN 

0799 


1 = 1- 1 


ISN 

0800 


IF( I *LE« 1 • AND* I GOP • NE • 0 ) GC TC 597 


ISN 

0802 




ISN 

0803 


IF( I .LT.C .AND. IGOP.EO. 0) GO TO 800 

0599 

ISN 

0805 


1F( J INX .£Q. I ) GO TO 555 

0600 

ISN 

0807 


MSK IP=0 

C 6 0 1 

ISN 

0808 


MTYPE= 1 


ISN 

0 809 


CALL T I TLES (MTYPE . MSKI P) 

0603 
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I SN 
ISN 

ISN 
ISN 
I SN 
I SN 
ISN 
ISN 
ISN 
ISN 
ISN 
I SN 


ISN 
ISN 
ISN 
ISN 
I SN 
ISN 
ISN 
ISN 
ISN 
I SN 
I SN 
I SN 
ISN 
ISN 
ISN 
ISN 
I SN 
I SN 
ISN 
ISN 


I SN 
I SN 
ISN 
ISN 
ISN 
ISN 
ISN 
ISN 
ISN 
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I-SN 
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ISN 
ISN 
I SN 
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ISN 
ISN 
ISN 
ISN 
I SN 
I SN 
ISN 
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0810 
06 11 
oa 1 2 

0813 
06 14 
C 8 1 5 
0 81 ^ 
08 17 
08 18 
ce 1 9 
0821 
0922 
C823 
<jB2 5 
0 8 26 
0827 
C 628 


C 8 30 
0621 
C632 
0833 
0835 
C636 
0837 
C 6 26 
0 8 40 
084 1 
0642 
0843 
084 4 
0e45 

0847 

0848 

0849 

0850 

0651 

0652 


G fi 54 
0865 
0856 
0657 
0e58 

0859 

0860 
0862 
0663 
08-84 

0865 

0866 
C 8 6 7 
0869 
0 6 70 

0871 

0872 

0873 
0674 


0870 
C 67 7 

0678 

0679 
0660 
0881 
C 882 
0 88 3 
0884 
C 8 86 
0887 
0889 
C«9l 
0 8 92 


5 X I X = XNORM 2* 4 MOOt ,X7(1 ) ) 

X I I X =XNQRM Z ( A MODE ,X7 <1 >> 

X I IXsXI I X-6S. 

1 = 1- 1 

CALL SETSMG ( AMCDE • 1 4 , 0. ) 

CALL LINE SG ( A MODE . I + I » X7 » Y7 ) 

C.ALL SfcTSMGtAMOOE ,1 4,3« J 

597 JNC T £ = C 

CALL 7IMTCK ( JK * JN DTE ) 

IF (IGCP.NE.C) GC TO 502 
GO TO 504 

EC 2 XIX=XNORMZ( AMCDE .SAVE4 ( 1 )) 

IF < I *GT • 1 IGCTCSG4 
X l I X = XNOK MZ( A MODE .SAVE 4 ( I 2+ 1 ) ) 

X I i X = X 1 1 X— 65 . 

5 C A l NOTE =C 

IF (TCP 1 .EG, 1 .AND. I .GT. 1 > 

1 CALL LEGNDG ( A MODE *237. .612# .94. 94 H* *** ONE OR MORE POINTS HAVE EXC 
2EEJED UPPER Y LI MI T AND ARE Cl VICED GY 10 BEFORE PLOTTING «***) 
Call CATAPT ( I TME 1 , I T WE2, I TME3 . I TME4 , XI X, I NOTE) 

CALL LEGNDG! AMODE .XIX, 170. .8 . 6H * APCGEE) 

CALL L£GN0G( AMODE. XI I X , 170. .9 ,9HP£R IGcE *) 

IF< KCCMP.EC. I ) GO TO 515 
MSKIP=C 
M TYPE =2 

CALL TITLES ( V T YPE , MSKI P) 

I F( I .LE • 1 .AND . IGOP.NE.O > GO TC 598 
CALL SfcT SMG( AMODE , 1 4 ,0 . ) 

CALL L IN5SG( AMODE. I *1 . XO , Y8 ) 

CALL SETSMG( AMODE . 1 4,3. ) 

598 JNO TE = 1 

CALL T IM TC K (JK.JNDTE) 

IF (IGOP.NE.O) GU TO 505 
GC TC £ C 7 

£ C £ CALL SETSMG( AMODE, 14 ,0. ) 

CALL LINE SGI AMODE. 12 + 1 , SA VE4 .SAVE 2) 

CALL SE T SMG ( AMOOE ,14.3. ) 

5C7 IND TE = 1 

IF< TOP 2 .£0 . 1 . AND. I . GT. 1 ) 

I C ALL L EGNDG( AMODE . 2 37. . 8 l 2. *94 , 94H* *** CNE OR MORE POINTS hAVE EXC 
2EEOEO UPPER Y LIMIT AND AH t DIVIDED BY 10 QEFCRE PLOTTING *♦** ) 
CALL DA T A PT ( I TME 1.1 TME2.I TMEJ.ITME4 .XlX.INDTE) 

CALL LEGM>G( AMOOE, XIX,170..8,8H* APOOEt ) 

CALL LEGNDG( AMOOE ,XI I X,1 70. ,9,9HP£R IGEE *> 

M SK IP = 0 
M T YP E = 3 

CALL TITLES ( M T YPE , M SK I P > 

IF< I .LE. 1. AND. IGOP.NE.O) GO TC 599 
CALL SETSMG (AMODE .14,0. ) 

CALL LINE SG ( A MODE , I 1 , X 9 » Y9 ) 

CALL SE T-SMGX AMOOE.L4..3*-! 

£99 J NO T £ =2 

CALL TIM7CK (JK.JNCTE) 

IF (IGOP.NE.O GC TO 508 
GO TC £ 10 

5C8 CALL SET SMG ( A N ODE ,14,0.) 

CALL L INESGtAMOOE 2*1 .SAVE* , SAVE3) 

CALL SET SMG( AMODE .14 ,J. ) 

£10 INJTE=1 

IF(TCP2.E0, 1 • AND . I • GT. I > 

1 C ALL LEGKOG( AMCDE.2J7. , d 1 2 . , 94 , 94 H* * ** CNE OR MCRE POINTS hAVE EXC 
2EEDED UPPER Y LIMIT AND ARE Cl VIDEO BY 10 BEFORE PLOTTING **¥*) 
CALL CATAPT ( I T ME l.IT Mfc 2.1 TME 3. 1 TME4 -.XI X.INDTEJ 
CALL L EGNDG ( AMODE .XIX, 170. ,8, BH* APCGEE ) 

CALL LEGNDG ( AMODE ,X II X, 17C. .9.9HPER I GEE *> 

5 1 5 CCNT INUE 

GC TC 60 C 

6 CO I TME 2= I T I ME 5 

ITM£4= IT I ME 6 
1 = 1-1 


0606 


C 6 1 2 

ce 1 3 
oe 14 

06 15 


cei7 


0621 

C623 


C 6 30 
C 6 3 1 

0632 

0633 


C 64 1 


0645 

0647 


0654 

0655 

0656 

0657 


C 665 


C 6 70 
067 1 
0672 


IF( I .LE. 1 .AND 
1=1+1 

I F( I .LT. t .AND 
IF ( J I NX .EO. 1 > 
M SK IP=2 
M T YPE= 1 


I GOP • NE . 0 ) 

GC 

TC 

697 


IGCP. EQ. 0) 

GU 

TO 

30 0 

0674 

GU TO 666 




0675 


0676 


Figure D-2. MAIN Routine (10 of 13) 


D-12 



C 6 70 


ISN 0893 
I SN 0 8 s. A 
ISN 0855 
ISN 0096 
ISN 0097 
ISN 0853 
ISN C 999 
ISN 0900 
ISN 0901 
ISN 0902 
ISN 0903 
ISN 09C5 
ISN C 50 6 
ISN 0907 
ISN 09C9 
ISN 0510 
ISN 09 i I 
ISN 0912 


ISN 091** 
ISN 0915 
I-S-N 0 516 
ISN 0917 
ISN 09 19 
ISN 0920 
ISN 0921 
ISN 0922 
ISN C 924 
ISN 0925 
ISN 0926 
ISN 0927 
ISN 0923 
ISN 0929 
ISN 0931 
ISN 0532 
ISN 0533 
ISN 0924 
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ISN 0936 


xsn ossa 
ISN 0 5 29 
ISN 0940 
ISN 0941 
ISN 0542 
ISN 0543 
ISN C544 
ISN 0546 
ISN 0947 
ISN 0948 
ISN 0949 
ISN 0950 
ISN 095 1 
ISN 0 95 3 
ISN 0954 
ISN 0955 
ISN 0956 
ISN 0957 
ISN 0958 


ISN 0960 
ISN 0961 
ISN 0562 
ISN 0563 
ISN 0564 
ISN C 9 65 
ISN 0966 
ISN C 56 7 
ISN 0560 
ISN C 570 
ISN 0971 
ISN 0573 
ISN 0575 


CALL TITLES ( M T YPE • MSKX P) 

666 XlX=XNORMZ( AMOOE .XI ( 1 ) ) 

* I IX=XNCRMZ (AMODE. XI ( I ) ) 

X I IX=X I I X-65. 

1=1- 1 

CALL SETSMG (a WOOL .1 4 , 0. ) 

CALL L INESG ( A MODE .I*1.X1.Y1) 

CALL SET SMG ( A MODE .14.3. ) 

657 JNOTE = C 

CALL TIMTCK (JK.JNDTE) 

IF ( IG CP • NE • 0 ) GU TO 602 
CO TO 60 4 

6C2 * IX=XNORM Z( AMOOE ,SAVfc4( 1 ) ) 

IFd.GT.DCiC TO604 

X I IX=XNCR*Z< Ay ODE, SAVt4< 12*11) 

X I I X = X I I X — 6 5 • 

604 INDTE=0 

IF(T0P1.EQ. 1 .AND. I .GT.l ) 

ICALL LEGND G ( A MODE *237* • S 1 2 • *94* 94 H* * * * CNE OR MORE POINTS HAVE EXC 
2£ EDED UPPER Y LIMIT AND ARE DIVIDED BY 10 BEFORE PLOTTING **♦*) 
CALL CA TAPT ( IT ME 1*1 TME2. I TME3 .1 TME4 ,XIX, INDTE) 

CALL LEGNDG { AMODE ,X I I X* 1 7C. *9 .9HAP0GEE *> 

C AL-L LEGNOGI AM CDC-, XIX , 4-7-0 P-ER4G8E) 

IFtNCOMP .EO . 1 > GO TO 615 
M SK I P-2 
V T YPE = 2 

CALL TITLES ( M TYPE . MSK I P ) 

I F C I «LE • I • A NO • I G OP • N.E . 0 ) GO TO 698 
CALL SETSMG (A MODE tl 4 , 0* ) 

CALL L INESG (AMODE .1*1 • X2 • Y2 ) 

CALL S6TSMG4AMODE , 14.3* i 
65 8 J ND T E = 1 

CALL TIMTCK (JK.JNDTE) 

IF (IGCP.NE.C) GO TO 605 
GO TC 607 

6C5 CALL SETSMG (AMODE, 1 4, 0. ) 

CALL L INE-SG4 AMODE »-I 2*4 . SA VE4 , SA VE2-) 

CALL SET SMG( AMODE . 14 , 3. ) 

6 C 7 I ND TE = 1 

IF ( TCP < .EG. I.AND.I.GT.l ) 

ICALL LE GND G ( A MODE .237. *812. .94. 94 H* 4*4 CNE OR MCRE POINTS HAVE EXC 
2EEDED UPPER Y LIMIT AND ARE Cl VIDEO BY 10 BEFORE PLOTTING *«**) 
CALL, DA TAPT (ITME1.IT ME 2 • I TME3.I TMEA »X I X . INDTE ) 

C ALL LEGNDC( AMODE .XI I X , 170. .9.9HAPCGEE *) 

CALL LEGNDG ( AMODE .X I X . I 70 . , 9 *9H* PERIGEE) 

MSK I P=2 
M TYPE =3 

CALL TITLES ( MT Y PE , M SK I P) 

IFt I *LE • l.ANO. IGCF.NE. 0) GU TO 699 
CALL SETSMG (AMODE . I 4 , 0. ) 

CALL LlNcSG( AMODE .1 *1 ,Xo .Y6) 

CALL SET SMG( AMODE . 1 4 ,3. ) 

699 J NOTE = 2 

CALL TIMTCK (JK.JNDTE) 

IF (1GUP.NE.0) GC TQ 6u8 
CO TO 610 

6(8 CALL SETSMG( AMODE .1 4,0. ) 

CALL LINESGt AMODE. 12*1 .SA VE4 , SA VEJ > 

CALL SET 5MG( AMODE • I 4 . 3. ) 

6 1C I NDTE= l 

IFlTGPc.cU. 1 .AND. I .GT.l ) 

‘CALL LEGNDG < A MODE .237. ,812. , 94 .94H**** CNE OR MORE PCINTS HAVE EXC 
2EECEC UPPER Y LIMIT AND ARE DIVIDED EY 10 BEFORE PLOTTING ***♦) 

CALL CA7APT(I TMt 1 , I TME2 . I TME3 .IT ME 4 ,X I X , I NOTE) 

CALL LEGNDG( AMCDE .XI X, I 70. , 9 ,5H* PERIGEE) 

CALL LEG ND G ( A MODE .XIIX.17C..9. 9H APO CEE *> 

6 15 CC N T I N LE 
GC TC 800 
7CC ITME2=ITIME5 
I TME < = MI ME 6 
1 = 1-1 

I M I .LE • 1 • AND • IGOP. NE. 0) GO TC 797 
1 = 1 ♦ 1 

IF( I .LT.6.AND. IGOR.EQ. 0) GO TO 300 
IF( J INX.EG. I) GO TC 777 
M SK IP = C 
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I SN 
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I SN 
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l-SN 
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1 041 
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ISN 

1 044 

ISN 
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10 46 
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1C47 

I SN 
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I SN 

l 049 

ISN 
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ISN 

1© 64 

ISN 

1052 

ISN 

1053 

I SN 

1054 


mtype= i 

CALL TITLES ( M TYPE , M SK I P ) 

777 X LX=XNGR MZ(AM,CDE .X3I t) ) 

XI IX=XNCRMZ( AMOOE*X3{ I ) ) 

X I I X = X I I X- 6 £ • 

1 = 1-1 

CALL SET SM G I A MODE .14*0. ) 

CALL L INE SG ( A MODE * I + l , X 3 • Y3 ) 

CALL SET SMG (AMCDE.14,3. ) 

797 JNOTE=C 

CALL T-IMTCK (JK.JNDTE) 

IF (IGOP.NE.O GO TO 702 
GC TO 7C4 

7C2 X I X = XNOR MZ ( A MCDE • S AVE4 ( 1 ) ) 

IF ( I .G T • DGCTC7C4 
> I IX = XNO«VZ (A MCDE , SAVE4I I 2* 1 ) ) 

X I I X=X 11 X-6 5 • 

7C4 INCTE=C 

IF( TOPI .EQ. 1 • AND* I* GT • I ) 

1C ALL LEGNDG I AMODE, 237. ,812. .94. 94H**** CNE CR MOPE POINTS HAVE EXC 
2EED60 UPPER Y LIMIT AND ARE DIVIDED BY 10 BEFORE PLOTTING *♦**) 
CALL CATAPT< I T NC l , I T ME 2 . I TME3.I TME4 • X I X , I NDT E ) 

CALL LEGNDG < A MODE .XI X. 1 70* .3* 8H * APCGEE) 

CALL LEONDGI AMODE ,XI IX . 1 70. , 9, 9HPER IGEE *) 

I F ( NCOMP * £0 • 1 ) GO TO 7 15 
MSKIPsC 
M TYPE = 2 

CALL TITLES ( MT Y PE . M SK I P) 

IF < I *LE • 1 • A NO * IGCP* NE • 0) GC TO 79d 
CALL SETSMG (AMODE , 1 4 ,0 • > 

CALL L INE SGI A MODE .14-1 .X4.Y4) 

CALL SET SMG( AMODE , 1 4 ,3. ) 

79d JNDTE=l 

CALL TIMTCK ( JK . JND TE ) 

IF (IGOP.NE.O) GC TU 7 05 
CO TO 707 

7 C 5 CALL SETSMG( AMODE .1 4,0. ) 

CALL LINE SGI AMUDE . 12 «■! . SA VE4 . SA VE 2 ) 

CALL SET SMGI AMODE . 1 4 . J. ) 

7C7 I ND TE= 1 

IFITOPi.EU. l.AND.I.GT.l) 

1 C ALL LEGNDG (AMODE .237* . 812* .94 ,94H**** CNE OR MORE POINTS HAVE EXC 
2EECEC OPPER Y LIMIT AND ARE DIVIOCD EY 10 BEFORE PLOTTING *♦**> 
CALL DA TAPTI I TME 1 . I TME2 . I TME3 .IT ME 4 . X I X . 1 NO T E ) 

CALL LEGNDG! AMODE. XIX. 170* .3 .6H* APOGEE) 

CALL LEG NDG ( A MODE .XIIX.170. ,9.9 HPER I C E E *) 

MSK IP=C 
MT YPE= 2 

CALL TITLES ( MT YPE ♦ MSKI P) 

IF( I *LE. 1 . AND. IGOP .NE.C ) GC TC 799 
CALL SET SMG ( A MODE .14.0.) 

CALL LINE SG I A MODE ,1 *1 .X5.Y5) 

CALL SET SMG ( AMODE .14.3.) 

799 JND TE- 2 

CALL T1MTCK (JK.JNDTE) 

IF (IGOP.NE.O) GO TO 70S 
GO TC 710 

7C8 CALL SETSMGI AMODE. 14,0. ) 

CALL L I Nt SGI A MCDE .12+1. SA VE4 ♦ SA VE3 ) 

CALL SETS MG I A M OOE .14,3.) 

710 IND TE = 1 

1 F ( TCP2.EO. 1 .AND. I .GT. 1 ) 

1 C ALL LEGNDG! AMOOE. 237. .812. .94 ,94H**** CNE OR MORE POINTS FAVE EXC 
2EEDED UPPER Y LIMIT AND ARE DIVIDED BY 10 BEFORE PLOTTING *♦**) 
CALL DA TAP T < I T 1 .IT ME 2.1 TME3.I TME4 .XlX. I NOTE ) 

CALL LEGNDG I AMODE »X I X • l 70. « 3. 8H * APOC-EE) 

CFLL LEGNdGI AMODE *X 1 I X . 1 7C . . 9 , 9 HPER I GEE *> 

715 CCNTINLE 
C 800 BACKSPACE 5 
8C0 BACKSPACE AS 
J IN X = 0 
JK = 0 
I-GOP =C 
1 = 0 

REWIND B 3 
ISTART=1 


C 756 


C 762 
0763 
0 7C4 
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C 7 73 
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END 
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LEVEi 20*1 X-AUG JXA 


OS/3O0 FORTRAN M 


COMPILER OfiTJOMS - NAME= M AI N. OPT = 0 2 • L I NECNT =58 • S I ZE =0 OOOK . 

SOURCE. EBCDIC *N0L 1ST. NO DE CK , LO AD . M A P, NOEDIT • ID. XREF 


W*AJ 

9003 


-SUBROUTINE DATA&T U74C-U.1Im£3*LT:mE3.IXmi--a. * Tk f I NOTE? 

ISN 

0003 


COMMON I T I ME 8(99 99 ) * I T I ME9( 9999 > • I T YPE { 9999 > . MHOURS! 100) .RANGES! 10 
1-0) .ERROR! 60 ) .ERROR 1UO). ERR0R2! 50 > . R ANGE7 £ 301 » LHntittgj f t AfftyyctTO) 

1.ABSICK30) . ABS I C2( JO) . I I DAT . IIDATl ,IH 

1 -SN 

0004 


C-OMNON E-RRLOI . E RRL02 *E WRLOJ. ERR HI 1 * ERRH I 2 ♦ ERRHI 3 

ISN 

0005 


COMMON AMODE ! 200 ) .CON. MANY.LOG 

FSN 

0006 


COMMON NSSi . N5SLSNS S3 rNSS4rNSSS • NSS6. A8. A7 .8 3*8 Sr i X-4V 

ISN 

0007 


DIMENSION M3NTH( 12) .NJNTH( 12) 

I SN 

0008 


DIMENSION XYl 2000 ) • DCOS ( 2000 )» RP ! 4000 ) • BOXX < 5 ) « BOXY < 5 ) 

ISN 

0009 


OATA S0RTST/4H / 

f-SN 

0010 


DATA 58XX/166.. 1182* .1 182. • 166* . 1 66 • / , BOX Y/ 92* • 92* . A. . 4 92^/ 

ISN 

00 1 1 


INTEGER A 5 • A7. 03. B5 

I-SN 

00-13 


TF4NSS4~*CO*-l-) S3 TO 114 

ISN 

0014 


MONTH! 1 1= 1 J1 

ISN 

0015 


MONTH! 21=226 

ISN 

0016 


MONTH! 3 ) = 331 

ISN 

0017 


MONTH! 4 >=430 

ISN 

0018 


MONTH! S ) =55 1 

I-SN 

001-9 


MONTH! 6 1*680 

ISN 

0020 


MONTH! 7 )=7Jl 

ISN 

0021 


MONTH! dl“831 

I SN 

0022 


MCNTH! 9 1=930 

ISN 

0023 


MON 7Ht 101=1031 

ISN 

0024 


MONTH! 1 1 >=1130 

rsN 

0025 


MONTH! 1-21=1251 

ISN 

0026 


MONTH! 11=101 

ISN 

0027 


MONTH! 2 1=201 

ISN 

0028 


NONTH! J)=301 

ISN 

0029 


NONTH! 4 1 = 401 

ISN 

0030 


NQNTH! 51=501 

T5N 

OO-H 


NONTH! o 1—60 1 

I SN 

0032 


MONTH! 7 ) = 70 1 

TSN 

0055 


NONTH! 8 >=801 

I SN 

0034 


NONTH! 9 )=90 1 

TSN 

0035 


MONTH! 1 0>=100l 

ISN 

0036 


NONTH! 1 1)=1 101 

TSN 

0057 


NONTHf 1-2 1*4-20 1 

ISN 

0038 


YEAR=F LOAT ! I TME 1 ) 

t-SN 

0039 


I-YEAR* YE AR/l 0000.0 

ISN 

0040 


I YE AR= I YEAR* 100 00 

ISN 

0041 


IF ! MOD! I YEAR.41.EQ.0) MO NTH( 2 ) =MONTH ! 2 ) ♦ 1 

ISN 

0043 


DO 459 LL=1 . 12 

I-5N 

0044 


MONTH! LL»*MONTH!LL ) *1 YEAR 

ISN 

0045 


NONTH! LL ) =NGNTH < LL ) ♦ I YEAR 

I-SN 

0046 

459 

CONTINUE 

ISN 

0047 


MONTH! 1 ) = NONTH! 1)110000 

I-SN 

0048 


CALL LfcGNDG! AMODE. 197. .105.. 17. 17HDATA DISTRIBUTION) 

ISN 

0049 

511 

ISTOP=0 

I-SN 

0060 


M=© 

ISN 

0051 


KD I F =0 

I-SN 

0052 


MDIF=0 

ISN 

0053 


LD IF =0 

ISN 

0054 


IF 1INDTE.EJ.1) GO TO 522 

ISN 

0056 


IF ! I NDT E . E 0. 3 ) GO TO 522 


0009 


0838 

0841 


0847 

Q flAB 

0849 

0850 

0851 

0852 

0853 

0854 

0855 
0855 

0857 

0858 

0859 
0880 
0861 
0883 
0863 
0354 
0865 
03 6 5 
0867 
0368 
0869 
Q870 

0871 

0872 

0873 

0874 
0 875 
0376 

0877 

0878 

0879 

0881 

0 882 

0883 

0884 
0 885 
Q886 
0887 
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D-16 



ISN 

0058 


J J=0 

I SN 

0059 


I XY — 0 

ISN 

0060 


XRRsO 

I SN 

006 1 


LM = 0 

ISN 

0062 


L = 0 

ISN 

006 J 


J = 0 

ISN 

0064 


K=0 

I SN 

0065 

522 

PRINT 121. ITME 1 • ITME2, ITME3, ITME4 

ISN 

0066 

1 21 

FORMAT ( IHO. 10X • 14HPERI00 PLOTTED • 1 OX • I 6 • IX • I 4* 5X . 1 6 • 1 X « IAJ 

ISN 

0067 


DAT A 9= FLO AT ( ITME 1 ) 

I SN 

0066 


WR I TE( 6*2) D AT A9.ITME1 

ISN 

0069 


DATA l = FLOAT ( I TME4) 

I SN 

0070 


DAT A 1* DAT Al/ 100.0 

ISN 

0071 


I D AT Al = OAT A 1 

ISN 

0072 


DATA lsFt.OAT< IOATA1 > 

ISN 

0073 


I F ( INDTE.EQ.O .OF. INDTE.EO. 1 ) DAT A 1 =DATA 1 ♦ 1 . 0 

ISN 

0075 


0ATA2sFL0AT< ITME2) 

I SN 

0 0 76 


DAT A2= D AT A2/ 100 -0 

ISN 

0077 


I DAT A2 = D A T A2 

I SN 

0078 


DATA2=FLOAT( IDATA2) 

I SN 

00 79 


KD tF= I T ME 3- 1 T ME 1 

I SN 

0080 


DO 567 I = 1 , 1 1 

I SN 

0081 


IF ( ITME1 .EQ. MONTH! I ) . AND. I TME3.E0 .NON T H ( I*| ) > KDIF-l 

ISN 

0083 


IF ( ITMEi.EG. MONTH! I ) • AND. ( ITME3-N0NTH! 1+1 IJ.EQ.l ) KDIF=2 

ISN 

0085 


IF ( IT ME 3 .E U. NONTH (14-1). AND.! MONTH! I ) — I TME 11* Ed • 1 ) KDIF=2 

ISN 

0087 


! F( I TME 3. EG* NINTH ( 14-1). AND* ( MONTH ( I 1-ITME 1 ) .EQ. 2) KDIF = 3 

ISN 

0069 


IF< ITME3.EQ. ( NO NTH ( I +1 ) + 1 > « AND. ( MONT H( - I 1-ITME 1 ) .EQ.1) KDIF=3 

ISN 

0091 


I F ( I TME 3* EG. ( NO NTH ( I 4- 1 1 4- 2 ) • AN D* ITME 1 .EG • MONTH ( I )) KDIFs3 

ISN 

0093 


IF! I TME 3. EG. NONTH ( 1*1 ) • AND* (MONTH! Il-!TMEl)*EQ*3> KOI F»4 

ISN 

0095 


IF( ITME 3. ED* < NONTH ( I 4- l ) 4-1 1 .AND* ( MONTH ( I )-I TME 1 ) • EO. 2 > KDlF = 4 

ISN 

0097 


IF ( I TME 3* EG. ( NONTH ( I 4- 1 ) 4-2 ) • AND* ( MONTH< I 1-ITMEl ) *EO* 1 1 K0IF*4 

ISN 

0099 


IF ( I TME 3* EG. ( NO NTH ( I 4- l ) 4-3 ) • AN D* I T ME 1 *E Q .MONTH ( I )) KDIF=4 

ISN 

O101 

56 7 

CONTINUE 

ISN 

0102 


I F ( ITME 1 .EG. MONTH! 12) . AND. I TM E3 .EG .NONT H ( 1 ) ) KDIF=1 

ISN 

0104 


IF! ITME 1 .EG. MONTH! 12). AND.! ITmE3-N0N t H( 1 )).EQ.l ) KDlFs2 

ISN 

0106 


IF! ITME J. EG. NONTH! 1 ) .AND. (MONTH! 1 2 ) - I TM E 1 ) . EG . 1 1 KD IF=2 

ISN 

0108 


IF! I TME 3. EG. NONTH! 1 ) .AND. ! MONTH! 121-1 TM El). EG. 2) KDIF=3 

ISN 

0110 


IF! ITME J. EG. ( NINTH! 1 )41 ) .AND. (MONTH! 1 2 ) - I TME 1 ) . EG . 1 ) KDIF=3 

ISN 

0112 


IF (ITME 3. EG. (NONTH! 1 ) 4-21 .AND. ITMEI.EQ .MONTH-4 12) ) K0IF=3 

ISN 

0114 


IF! I TME 3. EG. NONTH! 1 ) • AND. ( MONTH! 1 2 ) - I TM E 1 ) . EG . 3 ) KDIF = 4 

I SN 

0116 


I F ( I THE 3 . EG. ( NINTH! 1)4-1) .AND. ( MONTH! 121-ITME11.EQ.2) KDIF=4 

I SN 

0 118 


IF( I TME 3. EG. ( NONTH! 1)42) .AMD. ( MONTH! 12)-TTMEl).EQ.l) K0IFs4 

ISN 

0120 


I F ( I TME 3. EG. ( NONTH! 1)4-3) .AND. ITME 1 .EG. MONTH! 12) ) KDIF=4 

I SN 

0122 


DATA 1= D AT Al 4-24. *FLO AT ( K D I F ) 

I SN 

0123 


CALL 08 JCTG! AMODE .X I A. 45.. 1 15 1 • .90. ) 

I SN 

0124 


IF( NSS 5 .E G. 1 ) 30 TO 405 

ISN 

0126 


PRINT 123, XIX 

ISN 

0127 

123 

F CRM AT ( 1H • F 5 . 0 ) 

ISN 

0128 

405 

CALL SETS Mo ( AMODE *24,0. ) 

ISN 

0129 


CALL SUB JEG! A 1 JDt, DA’ r A2 , 0. .DATA i , 1 . ) 

ISN 

0130 


CALL LE oNDG ! AMODE , l 70. * 30. , 5. 5HHOU3S) 

ISN 

01 J 1 


CALL LL GNDG 1 A MODE • l 70. • 9. , 4 , 4HD ATE ) 

ISN 

01 32 


CALL LEGNDG! AMOCE * 1 70. . 82., 2, 2HXY ) 

ISN 

0133 


CALL LEGNDG! AMODE *170.,o7.,2, 2HL M ) 

ISN 

0134 


CALL LEGNDG! AMODE, 1 70 • • 52 • , 2, 2HRR ) 

ISN 

01 35 


CALL L I Nt 33 ( AMOCE, 5, BOX X « BOXY ) 
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0888 

0889 

aaoo 

0891 

0892 

0893 

0894 

0898 

0899 

0900 

0901 

0902 

0903 
0984 

0905 

0906 

0907 

0908 

0909 

0910 
091 1 
8912 


0913 


D-17 



ISN 

0136 


IF (NSS5.EU. 1 > GC TQ 600 


ISN 

01 36 


PRINT 123. X I X 


tSN 

0+39 


600 0A+A4»DATA2 

0949 

ISN 

0 1 40 


XDIF=DATA 1 -DATA 2 

0950 

ISN 

0-1-43 


I F { NSS 5.EQ. 1 ) GO TO 406 


ISN 

0143 


PRINT 123. XIX 




C 

SUBTRACT END AND START HPS FOR TOTAL 

0954 

ISN 

0 144 


406 JDIFF=XOIF 

0955 



C 

INTEGER DIFFERENCE CF HRS 

0956 

ISN 

0145 


OAT A5=DAT A2 

0957 

ISN 

0146 


PPINT 222, JDIFF 

0958 

I SN 

0147 


222 FORMAT (1H0.I5) 

0959 

ISN 

0146 


I OAT 2— DAT A2 


ISN 

0149 


CALL NUMBRGl A MODE »XNORMZ( AMODE, DAT A5 ) » 3 0. , 2 . T DA T2 ) 




c 

POINT TO COUNT FROM FOR NEXT TIME LABELS 

0965 

I SN 

0160 


JUlFF= JDIFF+1 

0966 

ISN 

0 15 1 


IF ( MANY . E Q. 1 ) JDIFF=JDIFF/4 


I SN 

0 153 


DO 700 1=1, JDIFF 

0967 

ISN 

0154 


DATA4=0A r A4F 1.0 

0969 

ISN 

0155 


DAT A5=D AT A5+ 1 . 0 

0970 

ISN 

0156 


I F ( MANY , E Q, 1 1 D AT A4 = DAT A4+3 .0 


ISN 

0 158 


I F ( MANY #EQ« 1 ) DAT A5 = DAT A5+J.0 


ISN 

0160 


IF ( OAT A 4 »G£ .24,0) GO TO 224 

0971 

ISN 

0162 


IF (DATA5.GE.DATA1 ) GO TO 594 

I SN 

0164 


227 I DAT 4=L)A T A4 


ISN 

0165 


CALL NUMB RG ( ANODE • X NOP M2 ( AMODE, DAT AS ) , 3 0. .2, IDAT4) 


I SN 

Olotj 


GC TO 700 

Q 974 

ISN 

0167 


224 DAT A4 = DAT A4 — 24 . 0 

0975 

ISN 

0166 


DAT 1 0=D AT A9F 1 .0 

0976 

ISN 

0 1 1>9 


WRITE! o»2) U AT 1 0 , I T ME 1 

ISN 

0170 


DC 456 LL=l , 12 

0977 

ISN 

0171 


LTME 4= DAT 1 0 

0978 

t SN 

0172 


LTMt£=L TME 4-MQNT-H <LL 1 

0979 

ISN 

0173 


IF < LTnt.tQ. 1 • AND.LL.EO. 12 ) D AT 1 0 = N0NTH ( 1 ) 

0980 

ISN 

0175 


IF (LTME.EU. 1 .AN0.LL.NE.12) 0 AT l 0=NONTH< LL ♦ 1 ) 

0 981 

ISN 

0177 


458 CONTINUE 

0982 

ISN 

01 78 


OATA9=OAT 10 

0983 

ISN 

0179 


XPOS=XNORMZ( AMODE .DAT A5 ) 


ISN 

0180 


XPOS=XHC3-13. 


I SN 

0181 


I DAT 1 0 =D AT 1 0 


ISN 

0182 


CALL 3BJCTG( A M(10fc,X IX, 45.. 11 82., 90. > 


ISN 

0 183 


CALL NUMBRG ( AMODE , X POS • 9 • • 6 . IDATio ) 


I SN 

0184 


CALL Q 6 JCTG! AMODE , X I X , 45 • , 1 1 5 1 . , 90 . ) 


I SN 

0185 


WPITE(o,2) DAT10.ITME1 


I«N 

01 86 


2 FORMAT ( 1H , F 9. 2 , I 6 > 


ISN 

0167 


IF ( CATA5.EU.0ATA 1 ) GO TO 594 

0987 

I SN 

0169 


GC TO 227 

0988 

ISN 

0 1 90 


7 CO CONTINUE 

0989 

I SN 

0191 


594 ID AT 4=D AT A4 


ISN 

0192 


XHRtND = XNCRMZ (AMODE , D AT A5 ) 


ISN 

0193 


CALL OBJCTG( AMOOE.XIX.45.. 1182. ,90. 1 


I SN 

0194 


CALL NUMBRG ( AMQUE , X HREND . 30.. 2. ID AT 4 > 


ISN 

0195 


CALL Ot» JC T G ( AMODE, XIX .45..1151..90.) 


ISN 

0196 


IF (IN0TE.EU.1) GO T 0 305 

0991 

ISN 

0196 


IF (1ST OP • E 0. 1 ) GC TQ 1 

0992 

ISN 

0200 


IF (INDTE.EO.J) GO TO 305 

0993 
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D-18 



ISN 

0202 

595 

K=K +■ 1 

0995 

ISN 

0203 


IF (K.GE.1600) GO TO 450 


ISM 

0205 

596 

«fcAO (43,302) ENDTPE.ITIME6(K),ITIME9(K.).ITYPE.IK) 


I SN 

02 0 

302 

FORMAT (IX, A4 « 3X , 1 1> , IX, 14, IX, 12) 


I SN 

0207 

305 

IF (MDIF.EQ.O) GO TO 612 

0999 

ISN 

0209 


IF (ENDTPE.EQ.SORTST ) GO TO 911 

1000 

ISN 

021 1 


IF! ITIME8(K> .LT.ITMEl ) GO TO 596 

IOOI 

ISN 

0213 


IF < ITIME8!K ) .EQ.ITME1 ) GO TO 597 

1002 

ISN 

0215 


Gti TO 202 


ISN 

0216 

597 

IF ( ITIME9!K) .LT.ITME2) GO TO 596 


I SN 

021 8 

202 

IF (K*EQ, 1) GO TO 598 

1045 

ISN 

0220 


IF ( ITIME8!K).E0.ITIME8(K-1 ) > GO TO 203 

1 006 

I SN 

0222 


LM=0 

100-7 

ISN 

0223 


IR«=0 

1008 

ISN 

0224 


I XV— 0 

1409 

ISN 

0225 


GO TO 598 

1010 

ISN 

0226 

203 

IF < ITIME9(K) .EU.ITIME9!K-1 ) ) GO TO 110 

1XU1 

I SN 

0220 


L M=0 

1012 

ISN 

0229 


I«R=0 

1013 

ISN 

0230 


I X Y = 0 

1014 

ISN 

0231 


GO TO 598 

1015 

ISN 

0232 

l 10 

IF ( ITYPE! X) .EO.TTYPEC K- 1 ) ) GO TO 595 

1016 

ISN 

0234 


IF C IT YPE < K ) , EQ * 1 ) GO TO 1 1 4 

1017 

ISN 

0236 


IF ( ITYPE (K) .EQ.9) GO TO 114 

1018 

ISN 

0238 


IF < ITYPE(K) • EQ • 2 ) GO To 115 

1019 

ISN 

0240 


IF ( ITYPE(K).EQ.3) GO TO lib 

1020 

ISN 

0242 


IP (IXY.EQ.l) GO TO 595 

HHH 

ISN 

0244 


GO TO 598 

1022 

ISN 

0245 

1 14 

IF (IRR.E0.1) GO TO 595 

1023 

I SN 

0247 


GO TC 598 

1024 

ISN 

0246 

115 

IF (LM.eO.l) GO TO 595 

1025 

ISN 

0250 

596 

IF ( ITIMEB(K) .LT.ITME3) GO TO 600 

1026 

ISN 

0252 


IF { IT1ME600 •GT.ITME3) GO Tq 611 

1047 

ISN 

0254 


If* ( IT IME8( <> .EO. ITME3) GO Tri 599 

1028 

ISN 

0256 


GC TO 600 

1029 

ISN 

0257 

599 

IF 1 IT I ME9( K) .GT.ITME4 ) GO To 811 

1030 

I SN 

0259 

600 

XTI=0. 5 

1031 

ISN 

0260 


LDIF=IT iMEti(K)— ITME1 

1032 

ISN 

0 261 


00 5o8 111=1,11 

1033 

ISN 

0262 


IF ( ITME1 .EO. MGNTH( I II ) . AND • I T ! ME6 ( K ) .EQ.NPNTH! I IT* 1 ) ) LDIF=l 

1 034 

ISN 

0264 


IF ( IT ME 1 .EQ. MONTH! Ill) • AND • ( I T IME6 ( K > -N )NTH( I I 1 + 1 ) ).EQ. 1 ) L0lF=2 

1035 

ISN 

0266 


IF ( (MONTH! X 1 1 ) -I T MEl ) , EO# 1 .AND, IT IME8! K ) • E 0, NO NTH < I I I + 1 ) ) LDIF = 2 

1036 

ISN 

0266 


IF! ITIME6(K). EQ.NONTH! I I I ♦ 1 ).ANO,< MONTH! Ill ) — IT ME 1 ) ,EO,2) L0IF = 3 


ISN 

0270 


I F ( I T I ME8 ( K ) .EQ. (NONTH! I T I ♦ 1 ) ♦ 1 ) • AND.! MONTH! III)-ITMEl).EQ.l) 




1 

- CO I F= 3 


ISN 

0272 


I F ( ITIME8 ( K) .LQ. (MONTH! IXX«-1)+2).AND.ITMEI .EQ. MONTH! Ill)) LDIF=3 


ISN 

0274 


I F ( ITIME8! X) .cQ.NONTH! I I I ♦ l ) .AND. ( MONTH! III)-ITME1).EQ.3) LDlF=4 


ISN 

0276 


IF(ITIMt6(K).6Q. (MONTH! I I I *1 > ♦ 1 ) .AND.! MONTH! I I I ) - I T ME 1 ) .EQ • 2 > 




1 

LDIF=4 


ISN 

0276 


IF( ITIME8(K) • E Q . ! NONT H( I 1 1 ♦ 1 ) *2 ) . AND • ( MONT H ( IIt)-lTMEl).EQ.l) 




1 

|_0IF=4 


ISN 

0280 


IF! IT I ME6 ( X) . EQ. (MONTH! I ! I + 1 ) +3 ) • AND. I TM E 1 • EQ . MONTH (III ) ) LDIF=4 


ISN 

0262 

566 

CONTINUE 

1037 

ISN 

0263 


IF ( ITME 1 . E Q .M'JN^H ! 12). AND. I T I ME 8 ( K ) .E 0. NONT H ( 1 ) ) LD IF=1 


ISN 

0285 


IF< ITME 1 .EJ. MONTH! 12). AND. ( IT I ME8 ( K ) -NONTH ( 1 ) ).EQ. 1 ) L0IF=2 


ISN 

0287 


IF! I TIME* (ft) .EQ.NOMTH! 1 ) .AND. (MONTH! 12)- ITME l ).EQ.l ) LDIF = 2 
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D-19 



ISU 

0239 



IF ( ITi-MEd ( K.J ,EQ,NQNTH( 1 ) .AND. (MONTH ( 12I-1TME1 ) . EQ.2 ) LDIF*3 


ISN 

0291 



IF( ITIME8(K) .EQ. (NONTH( 1 )*l ) .AND. ( MONTH ( 12I-ITME1 I.EQ.l ) L0IF = 3 



0203 



liLt-LTI ME a ( fCl . EQ >3 NUMXH L ± AND . I T_ME_1 .E Q.MQNTH (_1 2) 1 LD IF = 3 


ISN 

02 9b 



I F ( I TI MEQ ( K ) . EO.NONT H( 1 ) . AND. (MONTH ( 12)-ITME1 ) .EQ. 3 1 LDIF=4 


ISN 

0297 



4F4-ITIM&6< KUEQ* (NONTH( 1 ) ♦ 14 • AND* ( MQNTH(-12 1 — XTME-LJ .EQ.2 ) LOIE=4 


ISN 

0299 



IF( I TIMES ( K) .EQ. (NONTHC l M-2 ). AND. < MONTH ( 12)-I TMEll.EQ.il LDIF=4 


ISN 

0201 



IF( ITIME6( <) .EQ. (NONTHl 1 ) * 3 ) * AND. I TME1 . EQ. MONTH ( 121 1 LD IF =4 




C 

TO 

GET TOTAL NUMbER OF DAYS 

1036 

tSN 

0303 



PTXME9 4-K.I-1-T.IME0U. 1 *-24 0 041-0 IF 


ISN 

0304 



IHRS=ITI ME9 ( K ) / I 00 


ISN 

030b 



Ml NUTS* IT I ME 9 < K ) — I HRS* 10 0 


ISN 

030b 



YTI = FLOAT ( I HRS ) TFLP AT ( MINUTSI/60. 


ISN 

030-7 


669 

IF -4 I7-YP£(-K )*bO« l ) GO TO- 601 

1045 

ISN 

0309 



IF < ITYPEIK) .EQ.9) GO TO 601 

1046 

f-SN 

931 1 



GO TO 992 

1047 

ISN 

0312 


601 

JJ*J J+ l 

1046 

I SN 

0 313 



RR< J J ) = YT I 

1049 

ISN 

0314 



I RR= 1 

lO^O 

ISN 

0315 


6 02 

ITIME9< K) = ITIME9(K |-240 0*LDIF 


ISN 

031b 



GO TO 5y5 

1 0 53 

ISN 

0317 


612 

CONTINUE 


ISN 

0318 



MDIF=5 

1055 



c 

SHIFT FIRST DATE TO THE LEFT ON ALTITUDE PLOT TO ELIMINATE DATE OVERRUN 


ISN 

0319 



XPOSS=XNORMZ< AMODE. C AT A2 I 


ISN 

0320 



IF (NSSo .EQ. 1 ) XPGS5=XPOSS-38. 


ISN 

0322 



CALL NUMO RG ( AMO DE . XPOSS , 9. • 6 . ITMF.l ) 


I SN 

0323 



IF <INDT€.EQ.3> GO TO 811 

1 asft 

ISN 

0325 



IF (INDTE.EQ.1) GO TO 811 

1 059 

ISN 

0327 



IF (DATA1.EQ.OATA2I GO TO 811 

1060 

I SN 

0329 



Gr TO 596 

1061 

ISN 

0330 


99 2 

IF ( ITYPE(K) .EQ.2) GO TO 993 

4062 

I SN 

0332 



IF ( ITYPt (K1.EQ.3) GO TO 993 

1063 

13N 

0 334 



J — J ♦ 1 s, 

4964 

ISN 

0335 



XY(J)=YTI 

1065 

ISN 

033b 



I XV* 1 

1066 

I SN 

03 37 



GO TO oO 2 

1067 

I SN 

0336 


993 

L=Lfl 

1 066 

ISN 

0339 



DCOS (L 1 - YTI 

1069 

I SN 

0340 



LM=-1 

iQ70 

ISN 

0341 



GO TO o02 

1071 

ISN 

0342 


61 1 

IF (K.EQ.l) GO TO 812 

1072 

ISN 

0344 



PRINT J3J.JJ.L.J 


ISN 

0345 


333 

FORMAT ( 1H • 1 OX • 3 1 HR ANGE/R ANGfc KATE PASSES PLOTTED • 2X • I 6 ■ 5X • 24HNI N 

1074 





UNPACK PASSES PLOTTED. 2X. I6.5X, 17HXY PASSES PLOTTED .2X * I 6 > 

1075 

ISN 

934b 



GO TO 1 

1026 

ISN 

0347 


612 

CALL L E GNUG ( AM JOE .490. .48. .23.2 3HN0 DATA FOR THIS PEPIODI 


ISN 

0346 



PRINT 855 

1076 

ISN 

0349 


655 

FORMAT < 1H . 1 0 X . 25HN0 PASSES FOP THIS PERIOD) 

1 079 

ISN 

0350 



I ST 0P= 1 

1 080 

ISN 

0351 


615 

IF ( IS TOP • ME . 1 1 GO TO 595 

1061 

i-SN 

0 353 


1 

-IF (3.EQ.0I GO TO 100 

1062 

ISN 

0355 



XT I = 0. 75 

1083 

ISN 

035b 


1 1 

M-M* 1 

1064 

ISN 

0357 



CALL SETS MG ( AMOOE . 14.0.) 


ISN 

0356 



CALL NUMdRGt AMOOE.XYC M) .XT I .-l. 1H*| 


I SN 

0359 



CALL SETS MG ( AMO CE .14.3.) 



Figure D-3. Subroutine DATAPT (5 of 6) 
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ISN 

0360 


IF ( M« EQ • J i GJ TO 100 

1086 

ISN 

0362 


GO TO 11 

1087 


0363 

4 00 

M=0 

toaa 

ISN 

0364 


IF (L.EO.O) GO T 3 13 

1089 

I SN 

0366 


XT 1 = 0* 5 

1090 

ISN 

0367 

12 

M = M+1 

1091 

ISN 

0 368 


CALL SETS MG ( AMOCE • 1 4 • 0 • ) 


ISN 

0369 


CALL NUMB RG { A MODE • DCOS <M)#XTI»— 1» 1 H* ) 


ISN- 

0370 


CALL SET SMG ( A MODE * 1 4 • 3 • ) 


I SN 

0371 


IF { M.EG.L) GO TO 13 

1093 

ISN 

0373 


GO TO 12 

1096 

ISN 

0374 

13 

M— 0 

1095 

ISN 

0375 


IF (JJ.EQ.O) GO TO 111 

1096 

TSN 

0377 


XT I — 0 • 25 

1097 

ISN 

0378 

4 4 

M=M4 l 

1 098 

ISN 

0379 


CALL SETS MG ( A MODE • 1 4, 0< ) 


ISN 

0380 


CALL NUMBRG(AM00E#KR(M)*XTI.-1* 1H*> 


ISN 

0381 


CALL SETS MG ( A MODE • 1 4 t 3* ) 


I SN 

0382 


IF (M.EO.JJ) GO TO 111 

1 190 

ISN 

0384 


GO TO 14 

1 101 



C * BEGIN OCPCT3 SEQUENCE 


I SN 

0385 

911 

PPINT 3400 


ISN 

0386 

3400 

FORMAT (oOH REQUESTED T I ME SPAN TO BE PLOTTED EXCEEDS OBSERVATION T 




1 I ME ) 


ISN 

0 3 8 7 


REWIND B3 


ISN 

0388 


CALL EXITG(AMODE) 


ISN 

0 389 


STOP 


I SN 

0390 

450 

IT I ME8 { 50 ) = I TI ME 8 ( K — 1 ) 

1173 

ISN 

0391 


I T IME9( 50 )=1 TIME9I K-l ) 

1-174 

ISN 

0392 


ITYPCI 50 )=I TYPE (K-l > 

1 175 

I SN 

0393 


K =5 0 

1 1 76 

I SN 

0394 


GO TO d95 

1177 

ISN 

0395 

i 1 1 

REWIND S3 


ISN 

0396 


RETURN 

t 1 T9 

I SN 

0397 


ENO 

1180 


Figure D-3. Subroutine DATAPT (6 of 6) 
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LEVEL 20 . L l AUG 711 


HS/360 FORTRAN H 


COMP-lLER -OPTIONS - NAfctE= V A I N . 0PT=02 • L I NECNT = 58 • S I ZE = 0 OOOK , 

S3URCE.EBCDIC.NIDLTST, NODECK ,L0 AD ( HAP t NOEO IT , ID , XREF 


l&u 

0092 


SUBRQUTl-NE T I MT €K 1 JNO T E-4 

11-02 

1SN 

0003 


COMMCN I T I ME d ( 9 999 ) , I T f ME9( 9999 > , I TYPE ( 99991 ,MHCURS( 100 ) »PANGE5( 1 0 





1-01 TERROR <-50 ItERROR 1 < 50) , EPRCR2< 50) , R ANGE 713 Q)tT HOUR 2 ( 30 1.A0SI ClOO) 





1 . ASS I C 1(30), ASS t C2 ( 30) , I I DAT, I I DAT 1 , IH 

0009 

ISN 

0004 


COMMON 6RRL01 *ERRLP2,ERRL03 .EOPHI l , ERR HI 2 , ERRHl 3 


ISN 

0005 


COMMCN A MODE (200), CON, MANY ,LOG 


ESN 

OOOt. 


COMMON NS Si r-NSS 2*NSS3 » NS 54 . NSS5 , NS S6. A5 t-A7tB3tBSt 1 X-l-Y 


ISN 

0007 


INTEGER A5, A f » 83, 6b 


ISN 

0008 


IF( NSSs.EQ. 1 ) GG TO 407 


I SN 

0010 


PRINT 200, JK 

1187 

ISN 

001-1 

200 

FORMAT ( 1H , 1 IHDO I NOEX IS* ID 

LI 88 

ISN 

00 12 

407 

CALL LEGNDG( AMOOE ,237. , 900. , 70, 70HNUMBERS ON CURVE INDICATE UNI VER 





1SAL TIME IN HOURS ALONG THE TRAJECTORY) 


ISN 

0013 


DO |0 1 = 1, JK 

1191 

ISN 

0014 


I F< JNOTE .EQ.O.AND.LOG.EO.O) YlV=YNORMZ( A MODE • ERRORLI > ) 


ISN 

0016 


I F ( J ND TE ,£u,0 • ANC.LOG.EO. 1 ) Y I Y=YNOR MZ ( A MODE , ALOG 1 0 ( ERR0R( I ) ) ) 


ISN 

0018 


IF( JNDTE.EQ. 1 • AKD. LQ G. E Q. 0 ) Y I Y =YNORMZ ( AMQDE , ERROR 1 < I)j 


ISN 

0020 


IF ( JNDTE.EQ. l.ANO.LOG.EO.l) Y!Y=YNORMZ( AMODE , ALOG l 0 ( ERROR 1 < I ) ) ) 


I SN 

0022 


IF< JNOTE.EO. 2.ANO.LOG.EO.O) Y ?Y^V NO RNZ< 4 MODE tERROR 24--14) 


ISN 

0024 


IF( JNDTE .EQ. 2. ANC.LOG.EO. 1 ) YIY=YNOPMZ( A MODE , AL0G1 0 ( ERR OR2 ( I ) ) ) 


ISN 

0026 


X I X = XNOR MZ ( AMOOE, »ANGE5( I ) ) 


ISN 

0027 


IF(NSSb.£Q. 1 ) GC TO 4C8 


I SN 

0029 


PRINT 100,XIX»YIY 


I SN 

0030 

1 00 

FORMAT ( 1H » 2F o . 0 ) 


rsw 


408 

CONTINUE 


ISN 

0032 


XIX=XIX-J. 


ISN 

0033 


YIY=YI Y+15. 


ISN 

0034 


I F ( NSS 5 o EQ. 1 ) GO TO 4 09 


ISN 

0036 


PRINT 10 1 ,MH0URS( I ) 


ISN 

0037 

101 

FORMAT ( 1H ,13) 


ISN 

0030 

409 

CALL NUMBRG< AMOOE * X I X . Y IY* 2. MHO UR S ( I)) 


ISN 

0039 

10 

CONT INUE 

1209 

I SN 

0040 


RETURN 

1210 

ISN 

0041 


END 

1211 


Figure D-4. Subroutine TIMTCK 



LEVEL 20* 1 (AUG 71 i 


US/360 FORTRAN H 


CCBtPlLER OPTIONS - NAME= M Al N • OPT=02. L I NECNT = 58 • S I ZE = 0 0 OOK • 


I-&N 9002 
ISN 0003 


ISN 0004 
ISN 0005 
I 6N- 0006 
ISN 0007 
ISN 0 008 

ISN 0009 
ISN 0010 
ISN 0012 
ISN 0014 
ISN 0016 
ISN 0016 
ISN 0020 
ISN 0022 
ISN 0023 
ISN 0024 
I3N 0026 
ISN 0027 
ISN 0026 
ISN 002V 
ISN 0030 
ISN 0031 
ISN 0033 
ISN 0034 
ISN 0036 
ISN 0036 
ISN 0037 
ISN 0036 
ISN 003V 


SOURCE. EBCDIC. NOL 1ST. NO DECK ,L0 AD • P , NOED I T • I D • XREF 
SUBROUTINE ALT€K(KJ. JNOTE ) 

COMMON IT I ME 6 ( 9999 ) , I TI ME9! 9999 > • I TYPE ( 9999 > • MHCURS ( 100 ) • RANGES ( 1 0 
10 ) .ERROR! 50 ) . ERROR 1 ( 50 ) • ERRCR2! 50 ) • R ANGt 7 ( 30 ) • I HOUR 2 ( 30 I . ABSIC! 30) 
1 . AbSIC 1(30). A0S I C2 ( 301.1 I DAT . 1 1 DAT 1 • I H 
COMMON ERRLOl . ERRL02 . ERRLQ3 . ERRHJ 1 . ERRHI2.ERRHI 3 
COMMON AMODE ( 200) . CON. MANY .LOG 

COMMON N >S 1 • NSS2 • NSS3 • NSS4 . NSS5 « NSSb . A5 . A 7_. B3 »U5« IXlY 
INTEGER A5. A 7. 83. d5 

CALL LEGNDG! AMODE. 237*. 900..91, 91HNUMBERS ON CURVE INDICATE RADIAL 
1 DISTANCE FROM CENTER OF EARTH ALCNG THE TRAJECTORY KM* 1000) 

DO 10 1=1, KJ 

I F ( JNDTE .EQ.O .ANC.LOG.EO.O) Y I Y = YN ORMZ ( AMODE . ABSICt I ) ) 

IE (-JNDT-E-.EQ.0 .AND.LOG.EQ. 1 ) Y I Y= YNORMZi A MODE . ALQG 10 t ABSI.CMJJJ 
IF( JNDTE .EG. 1 . AND. LOG. EQ. 0 ) Y X Y = YNO RMZ ( AMODE, ABSIC 1 ( I) ) 

IF< JNDTE .ED. 1 . ANO. LOG. E Q. 1 ) Y I Y sYNDRMZ ( AMO DE , ALOG 1 0 ( ABS I C 1 ( I ) ) ) 

IF( JNDTE. ED. 2. ANC.LOG.EO.O) Y I Y =YND RMZ ( AMODE, ABS I C2 ( I ) ) 

IF< JNDTE .EQ. 2. ANO. LOG. EQ. t ) Y IY=YN0 RMZ( AMODE. ALOG l 0< ABSXC2( I ) ) ) 
HOUR2=FLOAT( IH0UR2( I ) ) 

XIX=XN0RMZ( AMOOE.HOUR2 ) 

IF (NSS5.EQ. 1 ) GC TO 410 
PRINT 100.xlX.YlY 

100 FORMAT! 1H ,Fo.0,F6.0) 

410 CONTINUE 
XIX=XI X-12. 

Y I Y«YI Y+15. 

IF (NSF5.EQ. 1 ) GO TO 411 
PRINT 101 . RANGE 7 ( T ) 

101 FORMAT ( 1H .F10.4) 

411 IRAN 7= R ANGE7 ( I ) 

CALL NUMbRG! AMO DE , X I X , Y I Y , 3 , IPAN7) 

10 CONTINUE 
RETURN 
END 


Figure D-5. Subroutine ALTCK 


U 4 3 

0009 

1219 

1223 

1234 

1235 

1237 

1238 

1239 
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LEVEL 20. I ( AUG 7U 


OS/ 360 FORTRAN H 


COMPILER OPTIONS - NAME = MAI N. QPT = 02 #L I NECNT=58 • SI ZE=OOOOK ♦ 

SOURCE# EBCDIC* NOL 1ST. NODECK #L0 AD » MAP# NOEDIT • ID# XREF 
I-6N M M S U » RO UT-I-N£ TITLES < MT-YPE ♦ MSJClO 

ISN 0003 COMMCN IT X ME 0 ( 9 999 ) • I T I ME9( 9999 ) • I TYPE ( 9999 ) • MHQURS (100) • RANGES I I 0 

10 ) .ERROR! 50 l #ERR0R1 ( SOI .ERROR 2 ( SO > , R AN 3E 7 ( 30 1 #J HQUR2 ( 30 ) • ABSICOO) 
I # ABSXC 1(30)# AdS I C2 ( 30 ) • I I D AT # I I DAT 1 # I H 


ISN 

0004 


CGMMGN ERPLQ1 *ERRL02#ERRLO3. ERR Hi 1 ♦ ERRH I 2 . ERRHI 3 

ISN 

0005 


COMMON A MODE (200) • CON# MANY#LOG 

HrN 

OOOto- 


C-GM-MGN N5S1 #NSS2»NSS3# N SS4 ^NS S5.NSS6 # Ab ♦ A7 *B3*S5 #-XX£Y 

ISN 

0007 


OATA I NCRD • JSY JU T /b • 6/ 

ISN 

0008 


DATA I YDARK/ 1/ , YLAbEL/1 •/ 

I SN 

0009 


INTEGER Ab# A7. 83. 8b 

ISN 

0010 


REAL*6 SNAME 

I SN 

001 1 


CALL 06 JCTG( AMODE# 204. . 1 45. # 1 1 6 1 • , 95 3 • ) 

ESN 

0012 


CALL PAGEGt AHOOE.O • 1 « 1 1 

ISN 

00 13 


IF ( MSK IP#EU#6 ) GO TO 2001 

ISN 

0015 

1 

IF ( MSK I P«E J»b) GO TO 50 

I SN 

0017 


GO TO ( 35»4b# 5b ) . MTYPE 

I SN 

0016 

35 

Y0*Y 0 1 

ISN 

0019 


YT=YT1 

ISN 

0020 


YGRI D=* YGR I D 1 

ISN 

0021 


IF(LOG#EC#0) I YCARK=YL A81 /YGR IO « 

I SN 

0023 


I F ( LOG • E C • 0 ) YLABEL=YLA81 

I SN 

0025 


F MT Y =F MT Y 1 

ISN 

0026 


GO TO bd 

I SN 

0027 

45 

Y8=Y82 

ISN 

0028 


YT=rY2 

I SN 

0029 


YGR IO= YGRI 02 

ISN 

0030 


I F ( LOG • E 0 • 0 ) I Y0ARK = YL A82/YGP ID2 

I SN 

0 0 32 


IF(LOG.EG.O) YLABEL=YLA62 

ISN 

00J4 


PMTY=FMTY2 

ISN 

003b 


GO TC 68 

t-sn 

0036 

53 

YB=Y83 

ISN 

0037 


YT= YTJ 

ISN 

0038 


YGRID= YGRIDJ 

ISN 

0 0 39 


I F ( LOG • E Q • 0 ) I YDARK = YL AB3/YGR ID3 

I 9N 

0041 


IF(LOG.EQ.O) YL ABEL = YL AB3 

ISN 

0043 


F MT Y =F MT Y 3 

ISN 

0044 

58 

CONTINUE 

I SN 

0045 


IF (MSKIP.EU.O) GO TO 60 

ISN 

0047 


IF (MSKIP. EU.2) GO TO 80 

ISN 

0049 


CALL SU8JEG( AMODE . 0 •* YB . 24 # ,YT ) 

ISN 

0050 


IF ( LOG • E 0# 1 ) CALL SETSMGt AMODE • 24 • 1 • ) 

I SN 

0 0 52 


CALL SETSMG( AMODE. 14. 0# ) 

ISN 

0053 


CALL GRIDG< AMOUE # . 5 . YGR I D# 2 . I YDARK ) 

ISN 

0 0 54 


CALL SLTSMG( AMODE. 1 4,3. ) 

ISN 

0055 


CALL L A6E LG ( AMODE #0# 1 . • 0 • 2 ) 

ISN 

005b 


GO TO VO 

I SN 

0057 

60 

CALL SUB J EG ( AMODE #XR #Yti#XL#YT 1 

I SN 

0058 


IF(LOG.EQ.l) CALL SETSMG ( AMODE • 24 , 1 . ) 

ISN 

00b© 


CALL SETS*G( AMODE# 14,0. ) 

ISN 

0061 


CALL G R 1 0 G ( AMODE #-XGRID#YGRIO# I XDARK # I YDARK ) 

ISN 

0062 


CALL 3tTSMG( AMOOE. 14,3. ) 

ISN 

006 3 


CALL LABELG( AMO CE • 0 • -XL ABEL » 0 • F MTX ) 

ISN 

0064 


GO TO 90 

ISN 

0065 

80 

CALL SUB J EG ( AMODE #XL»YB#XR#YT) 


1241 

0009 


12 50 
1251 


1255 
12 56 


1259 


1263 


Figure D-6. Subroutine TITLES (1 of 3) 
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ISN 

0066 


I F ( LOG « E 0 « 1 ) CALL SETSMG! AMODE . 24 , 1 . ) 

ISN 

0068 


CALL SETS MG ( AMODE. 14, 0. ) 

ISN 

9069 


CALL GRlOG! AMOOE.XGRID.YGPIC. IX DANK. I YU ARK 1 

ISN 

0070 


CALL SETSMG! AMODE, 1 4, 3. ) 

ISN 

00 7 1 


CALL L AbELG ( A MODE • 0 • XL AH EL ♦ 0. FMTX 1 

ISN 

0072 

90 

CALL LABEL G( AMO DE , l . Y LABEL • 0 • FMTY > 

ISN 

0073 


GO TO 200 

ISN 

0074 

50 

READ(INdND.lOO) CON.SNAME. I SAT. T PUN » LOG • MANY 

I SN 

0076 

1-00 

FORMAT (Al, IX, A8. IX. 15 .lX.l6.lX.il. IX. 11) 

ISN 

0076 


WRITE! JSYOUT.9100) CON.SNAME. I SAT, I RUN 

ISN 

0077 

9100 

FORMAT (1M .6X.A1.2X, A8 * IX, I5.2X. Io) 

ISN 

0078 


RE AD ( I NCRD • 10 1) XL.XR.YEl .YB2.YB3.YT1.YT2.YT3 

ISN 

0079 

101 

FORMAT ( 8F 10.0 ) 

ISN 

0080 


WRITE! JSYOUT , 30 0) XL.XP.YB1 .YB2.YBJ.YT1 .YT2.YT3 

ISN 

0081 

300 

FORMAT! IH , 8 F 1 0 • J) 

ISN 

0082 


READ! I NCRD , 1021 XGR I 0 . X L ABEL » FMTX 

ISN 

0083 

1 02 

FORMAT 12F10.0.F3.1) 

ISN 

0084 


I F ( X GR ID* NE . 0 « ) TXDARK = XLABEL/XGR!D 

ISN 

0086 


READ! I NCRD. 104) YGR ID1 , Y GP l D2 • Y GR !D3i YLAbl * YLAB2.YL- 



I 

l FMTY1 .FMTY2 .FMTY3 

ISN 

0087 

1 04 

FORMAT IOF 10.0,3(F3.1 *1X) ) 

ISN 

0088 


READ (I NCRD. 10 s) 5RRL01 . ER p LC'2 .ERRL03.ERRHI 1 , ERR HI 2 , 1 

ISN 

0 0 89 

1 C 5 

FORMAT t 6F10.0) 

I SN 

0090 


READ 1 0 3. 1 I DAT . I IDA- 1 

I SN 

0091 


POINT 9103. I IDAT.I IDATl 

ISN 

0092 

9103 

FORMAT 16X.215X.I6)) 

ISN 

0093 


READ 133. IH 

ISN 

0094 

133 

FORMAT! 14 ) 

ISN 

0095 

1 03 

FORMAT ( I 6. 1 X • I 6 ) 

I SN 

0096 


CALL SETSMG! AMOOE, 14,3. 1 

ISN 

0097 


CALL SETS MG! AMOOE .100,3.) 

ISN 

0098 


CALL SETSMG! AMOCE, 1 04, -.75) 

I SN 

0099 

2001 

CONTINUE 

ISN 

0100 


CALL SETSMG! AMOOE. 45, 1 .5) 

ISN 

0 101 


CALL LEGNOG! AMOOE. 562. , 800., 14, 14HPE»I0D COVERED) 

ISN 

0102 


CALL SETSMG! AMODE, 45. • 75) 

ISN 

0103 


CALL NUMBRG! AMOCE , 6 22. , 700. .o. I ID AT > 

ISN 

0 104 


CALL NUMBRG! A MODE , o92. .700..6.IIDAT1) 

ISN 

0105 

200 

CALL NUMBRG! AMOOE. 9 3 7. • 100O..6. I^UN ) 

ISN 

0 1 06 


CALL NUMBRG! A MODE, 913.. 980. .5. IS AT) 

ISN 

0107 


IF 1 MEKIP.EQ.5.0R.MSKIP. EQ.6) GO TO 30 

I SN 

0109 


IF(MTYPE.EU.2) GC TO 20 

ISN 

0 111 


I F ( MTYPE .ED. 3 ) GO TO 25 

ISN 

0 113 


CALL SETSMG! AMODE. 50. 270. 1 

ISN 

0 1 14 


CALL LEGNDG! AMODE • 142. . 900. • 1 J. 19HR ADI AL COMPONENT 1 

I SN 

0 115 


CALL SETSMG 1 AMODE . 50 . 0 • > 

I SN 

0116 


GO TC 30 

I SN 

0117 

20 

CALL SETSMG ( AMUDE » 50 .2 70. ) 

I SN 

0118 


CALL LEGNDG! AM3CE. 142. . 900. , 40. 




1 40HCOMP IN ORB. PLANE NORMAL TO PAD COMP KM) 

ISN 

01 1 9 


CALL SETSMG! AM3DE. 50. 0. ) 

ISN 

0120 


GO TO 30 

ISN 

0121 

25 

CALL SET SMu ( A MODE, 50, 2 70. ) 

ISN 

0122 


CALL LEGNDG! AMODE. 142. ,900. .29. 29HC0MP NORMAL TO ORI 

I SN 

0123 


CALL SETSMG! AMODE, 50, 0. ) 

ISN 

0124 

30 

CALL SETSMG! AMOCE. 45. 1 .5) 


1330 


1331 
1 332 
1333 


1364 


1366 


1371 


Figure D-6. Subroutine TITLES (2 of 3) 
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ISM 0125 

I-SN 04-26 
XSN 0127 

ISN 0128 
ISN 0129 
ISN 0130 
ISN 0^34 
ISN 0132 
ISN 0134 

ISN 0*35 
ISN 0136 


CALL LEGNUG ( AMQDb » 24 <t«, 980*. 32. 32H0 RBI T AL UNCERTAINTY ESTIMATE FOP 
1 I 

CALL SCTSMGt AMUOE, 75 ) 


CALL LEGNDG ( AM JOE • 26*: • . 1000., 83. 8 3HM T SS I ON 
1 D4 VISION-, GOOOARD SPACE FLIGHT CENTER 
CALL LEGNUGt AMJOE , 849. , 980 . , lb, 16H ( 

CALL SLTSMG( AM00E.45, l .5) 

CALL LEGNOG< AMTOE, 777. , 980. , 8, S NAME ) 

CALL SETSMGI AMQOE, 4b,. 75) 

IF( HSKIP.EG. l.HC. MS KIP.EQ.5.0P. MSK IP.EQ.6) 
CALL LEGNOGI AM00E.412. . 121. ,54, 

1 54 HR A D I AL DISTANCE FROM CENTER OF EARTH * 
40 RETURN 
ENO 


AND TRAJECTORY ANALYSIS 
RUN UATE4 
) ) 


GO to 40 

1000 KILOMETERS) 


1 399 

1400 


Figure D-6. Subroutine TITLES (3 of 3) 
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LEVEL 20.1 (AUG ?1J 


OS/ 360 FORTRAN H 


COMPILER OPTIONS - NAME= M AI N • UPT =02 *L I NtCNT = 56 • S I ZL = 00 00 K , 

SOURCE , EBCDIC* NHL 1ST » NODECK • LO AD » M A P * NO ED I T • ID. XREF 


ISN 

00 OP 

SUBROUTINE T APES <IBL API 

I AD 7 

I SN 

0003 

COMMON I T I ML 6 ( 9999 I • I T T ME9< 99 99 > • 7 T YPE ( 9999 ) .MH'JURS< 100) . RANGES ( 10 




10 > rtRMOR( 50 ) • ERROR 1(50). ERROR 2 ( 50 ) . R ANGE7< 30 ) • I HOUR 2 ( 30 > . AOSIC430) 




l * AQSXC 1(30). ABS I C2( 30 ) . I I D AT . 1 1 DA T 1 , I H 

0009 

ISN 

0004 

CCMMCN ERRL J1 .ERRL02 .ERRLOJ.6PRHI I , fcRRHl 2 . ERRHI 3 


ISN 

000b 

COMMON A MODE ( 200 ) * CON* MANY .LOG 




€ * BEGIN OCPtT 3 SEQUENCE 


ISN 

000b 

COMMON NS S 1 .NSS2.NSS3. NSS4 • NSSb • NSSb • A5 • A7 » B3 • B5 • IX IY 




C * END 0CPLT3 SEQUENCE 


ISN 

0007 

REAL * 6 TAPE.SCRTST 


I SN 

oooe 

DIMENSION TAPE! 25) . IFILE( 25) 

1406 

I SN 

0009 

DATA S3RTST/8H / 




C » BEGIN OCPLT3 SEQUENCE 


ISN 

0010 

INTEGER A5. A7. B3. B5 




C * END OCPLT3 SEQUENCE 


ISN 

00 1 1 

IF < IULAP.EO.O > GO TO bO 

1410 

ISN 

0013 

IF (IBLAP.EQ.l) GO TO 60 

141 L 

ISN 

00 1 5 

DO 10 1-1.25 

1412 

ISN 

001b 

READ 190. TAPE< I ) . IFILE< I ) 

1433 

ISN 

0017 

PRINT 9 190 .TAPE (I).IFILE(I) 


ISN 

0016 

9190 FORMAT ( 7X * Ab • 5X .12) 


ISN 

0019 

190 FORMAT <Ab.lX,l2> 


ISN 

0020 

IF (TAPc< I ).EQ.SO«TST) GO TC 20 

1415 

ISN 

0022 

10 CCNTINUE 

1416 

ISN 

0023 

20 DO 30 J=! .25 

1-4- 1-7 

I SN 

0024 

READ 190* TAPE (J).XFILE(J.) 

1416 

ISN 

002b 

PRINT 9190 *TAPE( J ) . IFILE( 3) 


ISN 

002b 

IF (TAPE(J).EQ. S. n RT ST ) GO TO 40 

1419 

ISN 

0026 

30 CONTINUE 

1420 

ISN 

002 9 

40 J=J-1 

1421 

ISN 

00 JO 

ICOP=l 

1422 

I SN 

0031 

I DC 2= 1 

1423 

ISN 

0032 

I VC= 1 

1424 

ISN 

00 33 

I DC= I 

1425 

ISN 

00 34 

GO TO 200 

1426 

ISN 

0035 

bO IC3P=I COP4-1 

1427 

ISN 

00 3b 

I F ( 1 CO P • GT • I F I L E ( I V C ) ) GO TO bOO 

1426 

I SN 

00 36 

GO TO 200 

1429 

I 9N 

0039 

bOO IVC=IVC+1 

1430 

ISN 

0040 

IF ( IVC.EO.I ) GO T C1 199 

1 431 

ISN 

0042 

IC0P=1 

1 A32 



C * BEGIN OCPLTJ SEQUENCE 




C * ADVANCE A 5 tq £OF 


ISN 

004 3 

BACKSPACE Ao 


ISN 

0044 

3540 READ ( A5* 104 , END=3550 > 


I SN 

004b 

1 04 FORMAT ( 7X , A4 ) 


ISN 

004O 

GO TO 3540 


ISN 

0047 

3b50 CCNTINUE 




C ♦ END OCPLTJ SEQUENCE 


ISN 

0046 

PRINT 250 ,TAPE< IVC) 

1433 

ISN 

0 049 

250 FORMAT ( 1H l , 26H0PERAT0R PLEASE MOUNT ' r AP E.2X.Ab.2X. 2 OHON A-5 AND HI 

1434 



IT START) 

1 435 



C PAUSE 


I SN 

0050 

GO TC 200 

1437 


Figure D-7. Subroutine TAPES (1 of 2) 
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FSN 

0051 

60 

IOC2=IDC2F 1 


1436 

ISN 

0052 


IF ( I DC 2 • GT • IF I LE ( I DC ) ) GO TO 70 


1439 

1-SN 

0064 


GO TO 200 


1 44Q 

I SN 

0055 

70 

I DC~ IDC+ 1 


1441 

ISN 

0056 


I F ( I DC • GT • J > GO TO 199 


1442 

ISN 

0058 


I DC2= 1 


1443 

ISN 

0050 


PRINT 251 • T APE ( IOC ) 


1444 

ISN 

0060 

251 

FORMAT ( 1H1 , 26H0PERAT0R PLEASE MOUNT TAPE . 2X . A6. 2X . 20HON 

B— 5 AND HI 

1445 




IT START) 


1446 



C 

PAUSE 



I SN 

0061 


GO TG 200 


1 446 

ISN 

0062 

190 

PRINT 201 


1449 

ISN 

0063 

201 

FORMAT { 1H1 .o2HALL REQUESTED TAPES HAVE BEEN PROCESSED - 

EXECUTION 

1-450 




1 TERMINATED ) 


1451 

ISN 

0064 


CAL.-L EXITG1AMODE) 



I SN 

0065 


GO TO 300 


1453 

ISN 

0066 

200 

RETURN 


1454 

ISN 

0067 

300 

STOP 


1455 

ISN 

0068 


END 


14.56 


Figure D-7. Subroutine TAPES (2 of 2) 
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LEVEL- 2 0.1 (AUG 71J 


JS/360 FORTRAN H 


COMPILER OPTIONS - N AME= MAIN* OPT =02 .L I NEC NT = 5d» S I ZE=OOOOK * 

SOURCE. EBCDIC. NOLI ST. NODECK ,LO AO . MAP . NOEDIT , 10. XREF 


X 644-0002 


SUBRQUT-TARi- »GFTAP(NF-) 



C 

*0CPLT3 SUBROUTINE 



e 

•»TH£~ P-O RPOSfc OF This SUBROUTINE IS TO BACKSPACE TO THE BEGINNING 



c 

* OF A COMPARE TAPE SEGMENT 

ISN 

0003 


COMMON IT IMtd< 9999) . IT I ME 91 9999 ) . I TYPE ( 9999) .MHOURSI 100 ) *RANG€&t 10 
10 ) .ERROR (50) • ERROR! < SO ). ERROR 2 ( 50) .P ANGE7I 30 ) . I HOUR 2( 30 ) • ABSI CC 30 1 
F**USIC M 3<H . AHSI€2( 30 ) . I I OAT, I I OAT 1 *-lH 

ISN 

0004 


COMMON ERRL01 . ERRL02 • ERRLQ3* EPRFl 1 • EPPHI 2 • ERR HI j 

ISN 

0005 


CeMMCN A MODE ( 200 ) • CON. MANY .LOG 

ISN 

0006 


COMMON NSS1 .NS32.NSS3. NS S4 . NS. e 5 .NS56.A5 .A7 .B3.B5. IXIY 

fSN 

0007 


OAT-A CHCK |/-l H6/ 

ISN 

ooce 


CATA CHCK2/IH7/ 

ISN 

0009 


INTEGER AS. A 7* 83. B5 

ISN 

00 10 

3040 

DO JOoO I = 1. 34 

ISN 

00 11 

305 0 

BACKSPACE NF 

ISN 

0012 


PEA 0 < NF. 3030) A3C0MP 

I SN 

0013 

3080 

FORMAT <1X. Al) 

I SN 

0014 


PRINT 4000. A3CCMP 

ISN 

0015 

4000 

FORMAT <8H A 3C0 MP= . A 1 ) 

ISN 

0016 


IF ( A3C JMP.EQ*CHCK 1 .OR. A 3COMP *E0* CHCK 2 ) GO TO 3040 

I SN 

ooid 


oo 3140 I = i • e 

ISN 

00 19 

3140 

BACKSPACE NF 

I SN 

0020 


READ (NF, 30B0) A3C0MP 

ISN 

0021 


IF (A3C0MP.EU.CHCK1 .OR. A3C0MP.EQ.CHCK2) GO TO 3040 

I SN 

0023 


DC 3150 1=1. 15 

ISN 

0024 

(j 

in 

o 

BACKSPACE NF 

ISN 

0025 


RETURN 

I SN 

0026 


END 


Figure D-8. Subroutine BSFTAP 


0 0 0 9 
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LEVEL 20.1 ( AUG 71 ) 


OS/ 360 FORTRAN H 


COMPILER OPTIONS - NAME= M Al N • QPT=02 . L I Nt CNT =58 . S I ZE=0 OOOK . 






30JRCE , EBCDI C,NQL 1ST, NODECK .LOAD ,MAP, NOEOIT . ID, XREF 


I SN 

0002 


SUBROUTINE T CON VO ( T I MDUT • lOUTIM.SEO 

QQQ3Q 



c 



00030 



c 



00040 



c 



00050 



c 

PUPPOSE 


00 060 



c 

THIS MODULE IS DESIGNED TO CONVERT CALENDAR TIME TO INTERNAL 

00070 



c 

000 

UNITS (CENTIOAY) AND VICE VERSA 

QQQflQ 



c 



00090 



c 



00100 



c 



00110 



c 

JYEAR 

ALWAYS EQUAL TO 1957 

00120 



c 

MCNTH 

ALWAYS EQUAL TQ 9 

00130 



c 

OAV 

ALWAYS EQUAL TO 18 

QQ-1AO 



c 

NBRD AY 

ARRAY CONTAINING T HE NUMBER OF DAYS PREVIOUS TO THE ITH 

00150 



c 


MONTH 

00160 



c 

T I MDUT 

NUMBER OF DUT*S fOOF. R/1B7S7 TO THE CALENDAR TIME 

00170 



c 

I OUT I M 

ARRAY CONTAINING T Ht YEAR , MONTH , DAY , HOUR AND MINUTE OF 

00180 



c 


CALENDAR TIME 

00190 



c 

J 

USED FOR LEAP YEARS 

003QC 



c 

M 

CONTAINS THE LAST THO DIGITS OF THE YEAR 

00210 



c 

SEC 

SECONDS OF CALENDAR TIME (LESS THAN A MINUTE) 

00220 



c 

K 

NUMBER OF DAYS FROM 9/18/57 TO JAN 1 OF THE CALENDAR YEAR 

00230 



c 

IDREF 

NUMBER OF DAYS FROM 9/18/57 TP THE CALENDAR DAY 

00240 



c 

T I MSEC 

NUMBER OF DAYS FROM 9/18/57 Tp CALENDAR TIME 

00250 



c 

cstc 

NUMdER OF SECONDS IN THE CALENDAR DAY 

nn^fA/y 



c 

L 

SET TT 0 IF NOT LEAP YEAR SET TO 1 IF LEAP YEAR 

00270 



c 



00280 



c 



00290 

I SN 

0003 


REAL 

T I MDUT , C SEC , T I M S EC 

00300 

ISN 

0004 


DIMENSION I OUT I V ( 5 ) • NOR DAY (12) 

00310 

I SN 

0005 


OATA NSRUAY / 0 . 3 l > B9 . 90. 1 dO, 1 b 1 . 1 8 1 , 2 1 2 , 2-7 3. 30*. 33A2 

00320 




4 , JYEAR/57/ , MO NTH/ 9/ » JD A Y/ IB/ , JDREF/O/ 

00330 

ISN 

0006 


IFITIMDUt .GT. 0.0) GO TO 10’ 

00 340 



c 



00350 



c 



00360 



c 


«* «««****».«*-**« « ****** • « «* * ** ** *** **** ***« *** «*****w********Q0370 



c 



00380 



c 

* CCMPUTES 

NUMBER OF CENTIDAYS BETWEEN THE REFERENCE DATE AND A 

00390 



c 

* REGUESTED 

DATE 

00400 



c 



00410 



c 



K400420 






00430 



c 



00440 

I SN 

OOOtt 


M = MOD ( I'JUTIMI 1 ) , 1900) - 1 

00450 

ISN 

0009 


I SUM = 

M* 100 00 TOUT 1M( 2 ) MOO* I OUT IM< 3) 

00460 

ISN 

00 10 


IF ( 1 5JM • GE • 5o 0 9 1 d ) GO TO 444 

00470 

ISN 

0012 


T I MDU T = 

- 1DJ 

00480 

ISN 

0013 


RETURN 


00490 

I SN 

0014 

444 CONTINUt 

OOSOO 

I SN 

001 5 


K- 1 0 4 


00510 

I SN 

0016 


iFlM.t-J 

• So) K = — 2 0 1 

00520 



C 

COK^OTLS NUMBER 1 F OATS FROM REFLRENCE DAT£ TO BEGINNING OF TEAR 

00530 

ISN 

ooia 


lF(M.t() 

• So • OR* M*EQ*57) GO TO 507 

00540 

I SN 

0020 


DC 1 I = 

SB, M 

00550 


Figure D-9. Subroutine TCONVO (1 of 3) 
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I SN 00 21 K = KH36;> 00560 

ISN 0022 IF <MGD(I.4) .E0. 01 K = K*1 00570 

i-sn 0024 i continue onsao 

C DETERMINES NUMUE R OF DAY 5 FROM THE BEGINNING OF REFERENCE YEAR TO THE 00590 
C BEGINNING OF YEAR FCF DATE REQUESTED 00600 

ISN 0025 567 J=0 00610 

ISN 0026 IF < IMOUI IUUTIMI 1 1 .4) .EO. 01 .AND. (I0UTIMI2) .GT. 2)1 0 = 1 00620 

C ADDS ANOTHER DAY TO COUNT IF THE REQUESTED DATE IS A LE'P YEAR ANO MONTH00630 
C IS 6REA-TEP THAN FEBRUARY 00640 

ISN 00 2b I = I0UTIMI2) 00650 

ISN 0029 IDPEF=K + NBR0AY( 11*1 OUT INI 3 1 ♦ J 00660 

C COMPUTES TOTAL NUMBER OF DAYS FCRM REFERENCE DATE TO REQUESTED DATE 00670 

ISN 0036 T IMSEC= < I DREF- JCKEF 1 *86400 ♦ I OUT t M < 4 1 * 360 0 ♦ I OUT l M( 5 1* 6 0 00680 

ISN 0031 T I M SEC = T I MSE C * SEC 00690 

C CC-MFUTES TOTAL NUMBER OF SECONDS FROM REFERENCE DATE TO REQUESTED. TIME 00.700 
t SN 0032 T IMDUT = TI MSEC/d£4. 0 60710 

C DETERMINES NUMBER OF CENTIoayS BETWEEN THE TWO DATES C0720 

ISN 0033 RETURN 00730 

c 00740 

r 00750 


C 444*4*4 *************** **************** ***************** ** ** ***44 ****** *4*0 0760 


c . 00770 

C * COMPUTES THE CALENDER U AT E GIVEN THE NUMBER OF CENTIOAYS FROM THE 00780 

C * REFERENCE DATE 00790 

C « 00800 

c ...».....*.«»**.*.*.****«******•*•*** *400 8 10 

c 66829 

C 00830 

ISN 0034 10 CONTINUE 00840 

ISN 0035 CSEC * ( DMJD ( TI *OUT *0.o0-7, 102 1 )*8o4D0 00850 

C DETERMINES NUMBER OF SECONDS LESS THAN A DAY 00860 

ISN 0036 SEC=DM0U<C5EC.601 1 00870 

ISN 0037 1 0U T IM< 5 )*OMOD< CSEC • 36 02 l/eO.O 00880 

ISN 0038 IOUTIMI 4 )=C3EC / 3600.0 00890 

C THE ABOVE THREE STATEMENTS DETEF M IHE . RE SPECT I VELY .THE NUMBER OF SECONDS. 00900 
C MINUTES AND HOURS CF THE REQUESTED DATE 00910 

ISN 0039 timoT=TIMdUT*,S787D-7 00920 

ISN 0040 K=NBMDA Y < MJNTH) ♦ JO AY ♦ I F I* ( T I MDT 7100.) 00930 

C DETERMINES NUMBER UF DAYS FROM THE BEGINNING OF THE YEAR OF THE REFERENCE00940 

C CATE 00950 

ISN 0041 IUUTIMI 1 ) =3 YEAR 00960 

ISN 0042 1 1 L= 0 00970 

ISN 0043 IFIMCDI IQOTIMIl 1,4) .EQ. 0) L = 1 00980 

ISN 0045 IF ( K . LE . ( 365 * LI) GO TO 12 00990 

ISN 3047 l OUT I M I 1 1 = 1 OUT I M ( l 1 Tl OlOOO 

ISN 0048 K = K — 3o 5 — L 01010 

ISN 0049 GO to 11 01020 

C THE ABOVE SEGMENT CALCULATES THE NUMBER OF VEAFS F R CM THE PREVIOUSLY 01030 

C CALCULATED NUMBER OF DAYS 01040 

ISN 0050 12 J=0 01050 

ISN 0051 DO 13 1*2.12 01060 

ISN 0052 IF (I . GE . 3 1 3=1 01070 

ISN 0054 IF I K .LE. INBRCAYIII * J*L1) GO TO 14 01080 

ISN 0 05o 13 CONTINUE 01090 

ISN 0057 1=13 01100 

ISN 0058 14 I OUT IM I 2 1 = 1-1 01110 


Figure D-9. Subroutine TCONVO (2 of 3) 
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C DETERMINES the month within that year 0X120 

isn 0059 if (i.ea.j) j=o 01130 

ISN 0004 IOUTIM<3> = K- NURD AY <1-1 > -<J«E> n i 1 An 

C DETERMINES THE NUMBER OF DAYS WITHIN THAT MONTH 01 ISO 

ISN 0002 RETURN 01 ISO 

ISN 0003 END 01170 


Figure D-9. Subroutine TCONVO (3 of 3) 
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L£V€L 20.1 CAUG 71) 


OS/ 36 0 FORTRAN H 


COMPILER OPTIONS - NAME - M AI N , OPT=02 ,L I NECNT.= 58 • S I ZE = 0 Q.00K • 

SOURCE .EBCDIC. NOLIST .NODECK, LOAD, MAP. NOEDIT .ID. XREF 


| Q QQ2 


A- b R E AO tT IM € 6 . 

IEXP3. RANI, IEXP4, I3EOF) 


p AP i , I EX P 1 . P AP 2. I E kp? ( 


€ 

C THE PURPOSE OF ThIS SUBROUTINE IS TO PASS VALID A5 VALUE b FROM 
C THE VECTOR COMPARISON TAPE TO THE CALLING PROGRAM 

C 



8908 


10 ) .ERROR (50) . ERROR 1(50). Ef?ROR2( 50 ) ♦ R ANGE 7 ( 3 0 ) . I HOUR 2 ( 30 ) • ABSI C( 30 ) 
l ,-AbSIG 1(30) r ABSIC2( 30) • I I OAT. 1 I DAT 1 • IH 

ISN 

0004 


COMMCN ERRLOl • E C RL02 • ERRL 03 . ERRH I 1 , ERR HI 2. ERR HI 3 

J-SN 

ooes 


COMMON AMOO€( 200)-. CON. MANY, LOG 

I SN 

0006 


COMMON NS SI . NSS2 .NS S3. NSS4 .NSS5.NSS6.A5.A7.B3.B5.IXIY 

ISN 

0GG7 


DATA CHCK l/lrto/ 

ISN 

0008 


DATA CHCK2/1H7/ 

I SN 

0 009 


OATA STT/4HYYMM/ 

I SN 

0010 


INTEGER A5 . A7. 03. U5 

I SN 

001 1 


I JEOF=0 

I SN 

0012 

3000 

CONTINUE 

ISN 

00 13 


IXI Y*I XI Y+l 

ISN 

0014 


IF (IXIY.GT.33) GO TO 2001 

ISN 

0016 

3001 

READ (AS. 30 I 0 ■ E RR=300 I ■ 

l ENO=3120) CHECK, I T I ME 5 , ITIME6. RANI, IEXP4. 

2RA01, IEXP1, RA 02, I EXP 2 , RA03, IEXP3 

I SN 

0017 

3010 

FORMAT (IX. At. Ib • IX, 14, 7X , F6.0* IX. 13. 15X. FB.6. IX. 13. 

12 X“* Fb.G, IX, 13. 2X, F 8 • 6 . IX, 13) 

ISN 

0018 


IF (CHECK. NE.CHCK1 ) GO TO 3100 

ISN 

0020 


ITIME5 = o*l 0000041 TIMES 

I SN 

0021 


GO TO 3130 

ISN 

0022 

3100 

IF (CHECK. NE.CHCK2) GO TO 3000 

ISN 

0024 


ITIME5 = 7* 1000004ITIME5 

ISN 

0026 


GO TC 3130 

I SN 

0026 

2001 

I X I Y =d o 

ISN 

0027 

889 

HEAD(A3,890,END=3129) CHC 

I SN 

0 0 2d 

890 

FORMAT ( IX, A4) 

ISN 

0029 


IF (CHC.NE.STT) GO TO 889 

I SN 

0031 


GO tq JOOO 

ISN 

0032 

2F20 

PRINT 3125 

ISN 

0033 

3125 

FORMAT (40H END CIF FILE ENCOUNTERED ON COMPARE TAPE) 

ISN 

0034 


I3EDF= 1 

ISN 

0035 


RETURN 

I SN 

0036 

21 29 

PRINT 3124 

ISN 

0037 

3124 

FORM AT ( 2X , ■ END OF FILE ENCOUNTERED UY READ 389 IN A5READM 

ISN 

093© 


I 3EOF= 1 

ISN 

0039 

3130 

CCNT INUE 

ISN 

0040 


RETURN 

ISN 

0041 


END 


Figure D-10. Subroutine A5READ 


0009 
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LEVEL -20*1 (AUG 71 J 


0S/3O0 FORTRAN H 


COMPILER OPTIONS - NAME = M AI N , OPT^O 2 ,L I NE CNT =58, S I2E=00 OOK • 


I-SN 0002 


SN 

00 

03 

SN 

0004 

SN 

0005 

SN 

0006 

SN 

OO 

07 

SN 

00 

OB 

SN 

00 

09 

SN 

00 

10 

SN 

00 

1 1 

SN 

00 

12 

SN 

00 

13 

SN 

00 

14 

SN 

00 

1 5 

SN 

00 

16 

SN 

00 

17 


SOURCE. EBCDIC. NOL 1ST. NODECK .LOAD, MAP. NOEDI T . I D , XREF 
SUBROUTINE BSREAO (I3YMO.I3HM, J-3T-VP ) 

C 

C THE PURPOSE OF THIS SUBROUTINE IS TO SUPPLY PROPER CALENDAR 
C DATE INFORMATION FROM THE WORKING FILE TO THE MAIN PROGRAM 

C 

COMMON I T I MEB (9999 ) , I T I ME9 ( 9999 > , ITYPE( 9999) ,MHOURS< 100),RANGE5( 10 
10) , ERROR (501. ERROR 1 ( 50 ) . ERROR 2 ( 50 ), RANGE 7(30) f IHOUR2UOJ ^ARSIXlTQ) 
1.ABSICK30) * A BS I C2( 30) , I I DAT. I I DAT 1 ,IH 
COMMON ERRL01,ERPL02,ERRL03,ERRHI 1 . ERRHI2.ERRHI 3 
COMMON A MODE (200) .CON. MANY .LOG 

COMMON NSS 1 . NSS 2 , NSS 3 , NSS 4 . NSS5 .NSS6.Ab.A7 .B3.B5. IX IY 
DIMENSION I 0U T I M ( 5 ) 

REAL *6 T I MOOT 

INTtGER*2 I3TYPE 

INTEGER A5, A7, B3 . B5 

READ (B5) A, IJTYPB, 8. C, TIMDUT 

I3TYP = I3TYPB 

CALL TCONVO (TIMDUT, IQUTIM, SEC) 

I3YMD = 1 00 * ( 10 0*1 OUT I M ( 1 ) ♦ I OUT IM(2))+I OUT I M< 3 ) 

I3HM = 100* IUUTIM< 4 ) ♦ IOU T IM(5 ) 

RETURN 

END 


0009 


Figure D-ll. 


Subroutine B5READ 
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